Magnetic Exchange Interactions in BaMn2As2: A Case Study of the Ji-J2-Jc 

Heisenberg Model 



D. C. Johnston/' 2 R. J. McQueeney/'^ B. Lake,3'4 A. Honecker,^ M. E. Zhitomirsky,^. ^ 

R. Nath/'* Y. Furukawa,-^'^ V. P. Antropov/ and Yogesh Singh^'''' 

^ Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA 
^Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011, USA 
Hahn-Meitner-Institut, Glienicker Str. 100, D-14109 Berlin, Germany 
Institut fiir Festkorperphysik, Technische Universitdt Berlin, Hardenbergstr. 36, D-10623 Berlin, Germany 
^Institut fiir Theoretische Physik, Universitdt Gottingen, D-37077 Gottingen, Germany 
''Service de Physique Statistique, Magnetisme et Supraconductivite, 
UMR-E9001 GEA-INAC/UJF, 17 rue des Martyrs, 38054 Grenoble Cedex 9, France 
'' Max-Planck-Institut fiir Physik Komplexer Systeme, 
Nothnitzer Strafle 38, D-01187 Dresden, Germany 
(Dated: January 26, 2013) 

BaMn2As2 is unique among Bar2As2 compounds crystallizing in the body-centered-tetragonal 
ThCr2Si2 structure, which contain stacked square lattices of 3d transition metal T atoms, since it 
has an insulating large-moment (3.9 /iB/Mn) G-type (checkerboard) antiferromagnetic AF ground 
state. We report measurements of the anisotropic magnetic susceptibility x versus temperature 
T from 300 to 1000 K of single crystals of BaMn2As2, and magnetic inelastic neutron scattering 
measurements at 8 K and ^^As NMR measurements from 4 to 300 K of polycrystalline samples. The 
Neel temperature determined from the x{T) measurements is Tn = 618(3) K. The measurements 
are analyzed using the Ji-J2-Jc Heisenberg model for the stacked square lattice, where Ji and J2 are 
respectively the nearest-neighbor (NN) and next-nearest-neighbor intraplane exchange interactions 
and Jc is the NN interplane interaction. Linear spin wave theory for G-type AF ordering and 
classical and quantum Monte Carlo simulations and molecular field theory calculations of x{T) 
and of the magnetic heat capacity Cmag(T') are presented versus Ji, J2 and Jc- We also obtain 
band theoretical estimates of the exchange couplings in BaMn2As2. From analyses of our x{T)i 
NMR, neutron scattering, and previously published heat capacity data for BaMn2As2 on the basis 
of the above theories for the Ji-J2-Jc Heisenberg model and our band-theoretical results, our best 
estimates of the exchange constants in BaMn2As2 are Ji « 13 meV, J2/J1 ~ 0.3 and Jc/Ji ~ 0.1, 
which are all antiferromagnetic. From our classical Monte Carlo simulations of the G-type AF 
ordering transition, these exchange parameters predict Tn ~ 640 K for spin S = 5/2, in close 
agreement with experiment. Using spin wave theory, we also utilize these exchange constants to 
estimate the suppression of the ordered moment due to quantum fluctuations for comparison with 
the observed value and again obtain S — 5/2 for the Mn spin. 



PACS numbers: 75.30.-m, 75.40.Cx, 75.50.Ee, 76.60.Es 
I. INTRODUCTION 

The observations of superconductivity up to 56 K 
in several classes of Fe-based superconductors in 2008 
(Refs. 1-4) have reinvigorated the high-Tc field follow- 
ing the discovery of high-Tc superconductivity in the lay- 
ered cuprates 25 years ago.^'^ Interestingly, the Fe atoms 
have the same layered square lattice structure as the Cu 
atoms do. Even though the maximum Tc of the Fe- 
based materials is far below the maximum Tc of 164 K 
for the cuprates,^ the Fe-based materials have generated 
much interest because the superconductivity appears to 
be caused by a magnetic mechanism^ as also appears to 
be the case in the cuprates. One of the many motivations 
for carrying out detailed measurements on the Fe-based 
materials is to see if these studies can clarify the super- 
conducting mechanism in the high-Tc cuprates for which 
a clear consensus has not yet been reached despite 25 
years of intensive research. 

Many studies of the magnetic properties of the Fe- 



based superconductors have been carried out."*'^ For the 
FeAs-based materials such as Bai_a;K3;Fe2As2 with the 
body-centered-tetragonal ThCr2Si2 structure, the mag- 
netic susceptibility x increases approximately linearly 
with increasing temperature above Tc or above the Neel 
temperature Tn of the nonsuperconducting parent com- 
pounds up to at least 700 K.^'^° In a model of local 
magnetic moments on a square lattice with strong an- 
tiferromagnetic (AF) Heisenberg interactions, this type 
of behavior is explained as being due to the measure- 
ment temperature (T) range being on the low-T side of a 
broad maximum in x{T) at higher temperatures.^ On the 
other hand, many magnetic measurements of the FeAs- 
based superconductors have been explained in terms of 
itinerant magnetism models, and indeed the consensus 
is pointing in this direction, although this view is not 
universal.^ 

In this context it is very useful to have a benchmark 
compound with the same ThCr2Si2-type structure and 
similar composition as many of the Fe-based supercon- 
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FIG. 1: (Color online) Crystal and magnetic structures 
of BaMn2As2.^^ The crystal structure is body-centered- 
tetragonal ThCr2Si2-type in which the Mn atoms form a 
square lattice within the a6-plane, the axes of which are ro- 
tated by 45° with respect to the a and b unit cell axes. The 
Mn atoms in adjacent layers are directly above or below each 
other along the c-axis. The magnetic structure is a G-type an- 
tiferromagnetic structure in which the magnetic moments of 
nearest-neighbor spins are antiparallel both in the afe-plane 
and along the c-axis. The ordered moment at T = is 
3.9(1) /^b/Mu." 



ductors, but for which a local moment model must be 
used to explain the magnetic properties. Such a com- 
pound is BaMn2As2 because it has an insulating ground 
state. ^^'^^ The crystal structure of BaMn2As2 is shown in 
Fig. 1.^'^ It is a small-band-gap semiconductor^^'^^ with 
an activation energy of 30 meV.^^ The electronic struc- 
ture calculations of An et al.^^ for the predicted con- 
ventional G-type (checkerboard) AF state give a small 
band gap of 0.1-0.2 eV, qualitatively consistent with the 
experimental value of the activation energy^^ that is ex- 
pected to be a lower limit to half the band gap. The 
anisotropic x of single crystals was previously measured 
at T < 400 K.^^ These data indicate that the com- 
pound is in a coUinear AF state at these temperatures, 
with the ordered moment direction along the c-axis, and 
with a Tn significantly above 400 K. From subsequent 
magnetic neutron powder diffraction measurements, the 
Neel temperature was determined to be Tn = 625(1) K 
and the AF structure was found to be a conventional 
G-typc (checkerboard) structure in all three directions 
as shown in Fig. 1, with an ordered moment direction 
along the c-axis in agreement with the x(T) data, and 
with an ordered moment ^ = 3.9(1) /ie/Mn at 10 K, 
where /zb is the Bohr magneton.^'' These characteristics 
are radically different from those of the similar FeAs- 
based metallic parent compound BaFe2As2 with the same 



room temperature crystal structure. BaFe2As2 has a 
much smaller ordered moment fi « 0.9 ^e/Fe and much 
smaller Tn = 137 K than BaMn2As2, the structure of 
BaFe2As2 distorts to orthorhombic symmetry below a 
temperature Tg « Tn instead of remaining tetragonal as 
in BaMn2As2, the ordered moment direction is in the 
afc-plane instead of along the c-axis, and the in-plane 
AF structure is a stripe structure (see the bottom panel 
of Fig. 2 below) instead of G-type.^ These large differ- 
ences between the magnetic properties of BaMn2As2 and 
BaFe2As2 evidently arise because BaMn2As2 is a local 
moment antiferromagnet whereas BaFe2As2 is an itiner- 
ant antiferromagnet. 

An intriguing aspect of the in-plane electrical resistiv- 
ity p{T) data for BaMn2As2 single crystals is that above 
~ 100 K the slope of the resistivity versus temperature 
changes from negative (semiconductor-like) to positive 
(metal- like). ^^'^^ The p{T) of a material can be written 
in an effective single carrier model as 
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where e is the magnitude of the electron charge, and 
n{T) and p{T) are respectively the effective conduction 
carrier density and the effective carrier mobility, respec- 
tively. Thus, a positive temperature coefficient of resis- 
tivity can be obtained for a band semiconductor if the 
increase in carrier concentration with increasing temper- 
ature is slower than the decrease in mobility with in- 
creasing temperature. Our ^^As NMR measurements in 
Sec. XI were in fact initially motivated in order to ad- 
dress this issue. As stated in that section, we found 
no evidence for a Korringa contribution to the ^^As nu- 
clear spin-lattice relaxation rate that would have indi- 
cated metallic behavior, and indeed we could interpret 
the data from 50 to 300 K in terms of a local moment 
insulator model. Furthermore, there is no evidence from 
the previously published neutron diffraction,^^ resistivity 
or heat capacity^^'^^ measurements for any phase transi- 
tion from a band insulator at low temperatures to a metal 
at high temperatures. Thus in the absence of experimen- 
tal data to the contrary, our present interpretation of 
the positive temperature coefficient of resistivity above 
~ 100 K is as discussed below Eq. (1) above. 

A related Mn-based compound is Sr2Mn3As202 which 
consists of Mn2As2 layers that are the same as in 
BaMn2As2, alternating along the c-axis with Mn02 lay- 
ers with the same structure as the Cu02 layers in the lay- 
ered cuprate superconductor parent compounds. Due 
to geometric frustration effects, the Mn spins in the 
Mn02 layers do not show any obvious long-range mag- 
netic ordering for T > 4 K, but the Mn spins in the 
Mn2As2 layers show long-range G-type AF ordering be- 
low Tn — 340 K with a low-temperature ordered mo- 
ment fj. = 3.50(4) ^B/Mn."'i'5 Thus, in both BaMn2As2 
and Sr2Mn3As202, the Mn spins in the Mn2As2 lay- 
ers exhibit the same G-type AF structure and signifi- 
cant reductions in the ordered moments from the value 
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/i = gSfiB = 5 /ie/Mn that would be expected for the 
high-spin S = 5/2 Mn+-^ ion with spectroscopic spht- 
ting factor g = 2. 

The main goal of the present work was to deter- 
mine the magnitudes of the exchange interactions in the 
fiducial compound BaMn2As2 and their signs, i.e., AF 
or ferromagnetic (FM). Experimentally, we extend the 
single-crystal anisotropic x{T) measurements from 300 
to 1000 K, significantly above Tn- We also report inelas- 
tic magnetic neutron scattering measurements at 8 K and 
''^As NMR measurements from 4 to 300 K on polycrys- 
tallinc samples. We analyze these data using the J1-J2- 
Jc Heisenberg stacked square spin lattice model for which 
we develop extensive theory. This model has also been 
investigated recently by other groups. "'^^"^^ We calculate 
the spin wave dispersion relations for this model. We re- 
port classical and quantum Monte Carlo simulations and 
molecular field theory calculations of x(T) and the mag- 
netic heat capacity C'mag(T). Wc extract the values of Ji, 
J2 and Jc by fitting our experimental data for BaMn2As2 
by these theoretical predictions for the Ji-J^-Jc niodel. 
From our classical Monte Carlo simulations of the heat 
capacity of coupled layers, we derive; a formula for Tn 
versus the exchange parameters which yields a Tn very 
close to experiment from the independently-derived ex- 
change constants, which indicates that the spin on the 
Mn ions is 5/2. We also utilize these exchange constants 
to determine from spin wave theory the suppression of the 
ordered moment due to quantum fluctuations for compar- 
ison with the observed value, and again arrive at the es- 
timate of 5* = 5/2 for the spin of the Mn+^ ions when the 
additional expected suppression of the ordered moment 
due to hybridization and/or to charge and/or magnetic 
moment amplitude fluctuations, which arise from both 
on-site and intersite interactions, are taken into account. 
Finally, we report band- theoretical calculations of J\ , J2 
and Jc for BaMn2As2. 

The applicability of the local moment Heisenberg 
model to a specific compound depends on the degree of 
variation of atomic magnetic moments and interatomic 
exchange parameters found from electronic structure cal- 
culations for the relevant magnetic ordering configura- 
tions. Such variations are usually found to be small in 
magnetic insulators. In the case of BaMn2As2, our band 
theory analysis in Sec. XII indicates that insulating char- 
acter is conserved for both the Ncel and stripe antiferro- 
magnctic structures, as observed, with a tiny metallicity 
appearing in the ferromagnetic case. As noted above, 
An et al. previously estimated that the band gap is 0.1- 
0.2 eV from electronic structure calculations for G-type 
AF order in BaMn2As2.^^ Moreover, our direct calcula- 
tions of the atomic magnetic moment and exchange cou- 
plings for different spin configurations demonstrate that 
the ordered moment variations do not exceed 10-12%, 
while the exchange coupling variation is only about 5%. 
The largest change appears for the ferromagnetic state 
which due to its high energy is expected to contribute lit- 
tle to thermodynamic properties. Finally, our determina- 



tion of a self-consistent set of antiferromagnetic exchange 
coupling parameters in BaMn2As2 from both static and 
dynamic experiments confirm the validity of our analyses 
in terms of the local moment Heisenberg model. 

The remainder of the paper is organized as follows. 
The experimental details for the sample preparation and 
characterization of BaMn2As2 and for the various mea- 
surements are given in Sec. II. The Ji- J2- Jc Heisenberg 
model is introduced and defined in Sec. HI. The inelas- 
tic neutron scattering measurements of polycrystalline 
BaMn2As2 and the analysis of these data are presented 
in Sec. IV. This includes the presentation of spin wave 
theory for the Ji-J^-Jc model of the G-type antiferromag- 
net in Sec. IV A that is used to fit the neutron data and 
to later obtain an estimate of the spin wave contribution 
to the heat capacity at low temperatures in Sec. VIII B 
and to analyze the ''^As nuclear spin-lattice relaxation 
NMR data below Tn in Sec. XI B. The high-temperature 
anisotropic magnetic susceptibility measurements of sin- 
gle crystals of BaMn2As2 are presented in Sec. V. The 
predictions of molecular field theory (MFT) and related 
topics for the J\-J2-Jc Heisenberg model are given in 
Sees. VI, VII and the Appendices. Comparisons of the 
MFT predictions with our experimental x{T)i Cp{T) and 
ordered moment /i(T) data for BaMn2As2 are given in 
Sec. VIII. In this section we also calculate the spin wave 
contribution to the heat capacity at low temperatures 
assuming a negligible anisotropy gap in the spin wave 
spectrum and compare this contribution with the experi- 
mental heat capacity data at low temperatures. Classical 
and quantum Monte Carlo simulations of Tn, x(T) and 
CmagiT) in the J\-J2-Jc Heisenberg model are presented 
in Sec. IX and comparisons with the experimental Tn 
and x(T > Tn) data arc carried out in Sec. X. The 
NMR measurements and analysis are given in Sec. XI, 
and our band-theoretical estimates of the exchange cou- 
plings in BaMn2As2 are presented in Sec. XII. Our spin 
wave theory results for the suppression of the ordered 
moment are given in Sec. XIII. From a comparison with 
the experimental ordered moment, we infer that the spin 
on the Mn ions is 5/2. A summary of our results and of 
our most reliable values of the Ji, J2 and Jc exchange 
constants and of the spin value derived for the Mn ions 
in BaMn2As2 is given in Sec. XIV. 



II. EXPERIMENTAL DETAILS 

A 25-g polycrystalline sample of BaMn2As2 was pre- 
pared by solid state synthesis for inelastic neutron scat- 
tering (INS) measurements. Stoichiometric amounts of 
Ba dendritic pieces (Aldrich, 99.9%), Mn powder (Alfa 
Aesar, 99.9%), and As chunks (Alfa Aesar, 99.9%) were 
ground and mixed together in a He-filled glovebox, pel- 
letized, placed in a 50-mL AI2O3 crucible with a lid and 
sealed in a quartz tube under a 0.5 atm pressure of Ar 
gas (99.999%). The tube was placed in a box furnace and 
heated at a rate of 50 °C/h to 575 °C and kept there for 



4 



24 h. The furnace was then heated at 100 °C/h to 800 °C 
and kept there for 48 h before coohng to room temper- 
ature by turning off the furnace. The quartz tube was 
opened inside the glovebox and the product was ground 
and mixed thoroughly and pelletized again. The pellet 
was placed in the same crucible and sealed again in a 
quartz tube. The quartz tube was heated in the box fur- 
nace at 100 °C/h to 850 °C and kept there for 24 h and 
then heated at 100 °C/h to 900 °C and kept there for 
24 h, followed by furnace-cooling at ~ 300 °C/h to room 
temperature. The resulting product was ground and 
pelletized and the above heat treatment was repeated. 
The resulting product was characterized by x-ray pow- 
der diffraction and the majority phase (ss 83%) was found 
to be the desired ThCr2Si2 structure. The major impu- 
rity phase was identified to be tetragonal Ba2Mn3As2 0. 
BaMn3As20 is an insulator, shows low-dimensional mag- 
netic behavior with a broad maximum in x{T) at 100 K 
and antiferromagnetic ordering at « 75 K.^^'^^ Prom x- 
ray diffraction measurements the weight fraction of this 
impurity phase in the INS sample was estimated to be 
« 17%. The INS spectra at 8 K and at 100 K (not 
shown) exhibited no noticeable differences. Since 100 K 
is well above the purported ordering temperature of the 
impurity phase, this eliminates any concern for serious 
contamination of the magnetic INS data by this phase. 

About 20 g of the above material was used for INS mea- 
surements. About 5 g of the polycrystalline material was 
used to grow single crystals. About 3 g of polycrystalline 
BaMn2As2 and 20 g of Sn flux (Alfa Aesar, 99.999%) 
were placed in an AI2O3 crucible and scaled in a quartz 
tube. The crucible was heated at 250 °C/hr to 1000 °C 
and kept there for 24 h and then cooled at 5 °C/h to 
575 °C and kept there for 5 h at which point the Sn fltix 
was centrifuged off to give isolated single crystals of typ- 
ical dimension 5 x 5 x 0.2 mm^. Energy-dispersive x-ray 
(EDX) measurements using a Jcol scanning electron mi- 
croscope confirmed the composition of the crystals to be 
BaMn2As2. 

For the INS measurements, the powder sample of mass 
« 20 g was characterized for phase purity by x-ray pow- 
der diffraction as discussed above. The INS measure- 
ments were performed on the Pharos spectrometer at 
the Lujan Center of Los Alamos National Laboratory. 
Pharos is a direct geometry time-of-flight spectrometer 
and measures the scattered intensity over a wide range 
of energy transfers hco and angles between 1 and 140° al- 
lowing determination of the scattered intensity S{Q,uj) 
over large ranges of momentum transfer UQ and The 
powder sample was packed in a flat aluminum can ori- 
ented at 135° to the incident neutron beam, and INS 
spectra were measured with incident energies Ei of 150 
and 200 meV. The data were measured at a temperature 
T = 8 K, well below the antiferromagnetic ordering tem- 
perature of 625 K.^"^ The time-of-fiight data were reduced 
into Uu) and scattering angle (26) histograms and correc- 
tions for detector efficiencies, empty can scattering, and 
instrumental background were performed. 



The high-temperature anisotropic xC^) measurements 
of a BaMn2As2 single crystal took place in a physical 
properties measurement system (PPMS, Quantum De- 
sign, Inc.) at the Laboratory for Magnetic Measurements 
at the Helmholtz Zentrum Berlin fiir Materialien und En- 
ergie. For these measurements the vibrating sample mag- 
netometer option was used. Data were collected with the 
magnetic field applied both parallel and perpendicular 
to the Mn layers. For field pointing within the afe-plane 
a sample of mass 15.31 mg was used. The sample had 
to be cut for field parallel to c and the sample weight 
was 12.058 mg. For all measurements a constant mag- 
netic field = 3 T was used while the temperature was 
varied between 300 and 1000 K. To achieve these temper- 
atures an oven set-up provided as an option by Quantum 
Design was utilized. The crystal was fixed on a zirco- 
nia sample stick containing a wire system that acts as a 
heating element. The sample was glued on the stick with 
heat-resistant cement and wrapped in low emissivity cop- 
per foil to minimize the heat leak from the hot region to 
the surrounding coil set. The measurements took place 
with heating rates of between 1 and 2 K per minute. The 
magnetic moment of the empty sample holder, 7.63 mg 
of cement and of the copper foil was measured separately 
and subtracted from the data. 

The NMR measurements were carried out on a poly- 
crystalline sample using the conventional pulsed NMR 
technique on ^^As nuclei (nuclear spin 7 = 3/2 and gyro- 
magnetic ratio '^^7/27r = 7.2919 MHz/T) in the tempera- 
ture range 4 < T < 300 K. The measurements were done 
at a radio frequency of about 52 MHz. Spectra were ob- 
tained by sweeping the field at fixed frequency. The ''^As 
nuclear spin- lattice relaxation rate 1 /Ti was measured by 
the conventional single saturation pulse method. 



III. THE J1-J2-JC HEISENBERG MODEL 

A bipartite spin lattice is defined as consisting of two 

distinct spin sublattices in which a given spin on one sub- 
lattice only interacts with nearest-neighbor (NN) spins 
on the other sublattice. In the FeAs-based supercon- 
ductors and parent compounds, when the magnetism is 
analyzed in a local moment model, the magnetic lattice 
is found not to be bipartite."* In addition to the in-plane 
(Ji) and out-of- plane (Jc) NN inter-sublattice interac- 
tions, in-plane diagonal next-nearest-neighbor (NNN) 
intra-sublatticc interactions J2 are also present along 
both diagonals of each square, as shown in Fig. 2. The 
spin Hamiltonian in the Ji- J2- Jc Heisenberg model is 

n = Ji^Si-S,+ J2^Si-S, (2) 

NN NNN 

+ JcJ2Si-Sj+gfiBHj2st, 

c i 

where Si is the spin operator for the ith site, g is the 
spectroscopic splitting factor ((^-factor) of the magnetic 
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FIG. 2: (Color online) CoUinear commensurate in-plane mag- 
netic structures in the J\-J2-Jc model for the square lattice 
antiferromagnet. The top panel shows the G-type (Neel or 
checkerboard) AF structure where nearest-neighbor spins are 
aligned antiparallel. The bottom panel shows stripe-type or- 
dering, along with the definitions of the in-plane exchange 
constants Ji and J2. A J2 interaction is present for both di- 
agonals of each square. According to Eqs. (4), the G-type 
in-plane ordering is favored if J2 < Ji/2, whereas the stripe- 
type ordering is favored if J2 > Ji/2. By examining the 
bottom panel, one sees that the stripe magnetic structure 
consists of two interpenetrating G-type magnetic structures, 
each respectively consisting of next-nearest-neighbor spins. 

moments, /xb is the Bohr magneton and H is the magni- 
tude of the magnetic field which is applied in the -\-z di- 
rection. Throughout this paper, a positive J corresponds 
to an antiferromagnetic interaction and a negative J to 
a ferromagnetic interaction. The indices NN and NNN 
indicate sums restricted to distinct spin pairs in a Mn 
layer, while the index c indicates a sum over distinct NN 
Mn spin pairs along the c axis. This is the minimal model 
needed to explain our INS results below for BaMn2As2. 

The classical energies of collinear commensurate or- 
dered spin configurations in this model with H = are 
analyzed as discussed in Ref. 4. We consider four com- 
peting magnetic structures in the J1-J2-JC model. One 
is the simple FM structure. The other three are two AF 
stripe structures and the G-type (Neel) structure shown 
in Fig. 2. By definition, the NN spins in alternate layers 
are aligned antifcrromagnetically in the G-type AF or- 
dered state, whereas the stripe state can have either AF 
or FM spin alignments along the c-axis which depend on 
the sign of Jc- The classical energies of these states for 
H = are'' 

Efm = NS^{2Ji + J, + 2J2) 
Estripc = NS\-2.h±Jc) (3) 
Eg = NS\-2Ji- JC + 2J2), 

where N is the number of spins S and a factor of 1/2 has 
been inserted on the right-hand sides to avoid double- 



counting distinct pairs of spins. The ± signs in the ex- 
pression for the stripe phase arise due to the possibili- 
ties of either antiferromagnetic (— sign) or ferromagnetic 
(-1- sign) alignment of adjacent spins along the c-axis. 
From these expressions, the in-plane G-type AF mag- 
netic structure observed in BaMn2As2 is lower in energy 
than the stripe structure if 

Ji > (G type AF) 

Ji > 2J2. (4) 

In order that G-type AF ordering occurs along the c-axis, 
one also requires that 

Jc > 0. (G type AF) (5) 

These results place restrictions on the exchange cou- 
pling parameter space that is relevant to the G-type AF 
ordering observed in BaMn2As2. Equations (4) and (5) 
require both Ji and Jc to be positive (antiferromagnetic), 
but J2 can have either sign as long as it satisifies the sec- 
ond of Eqs. (4). The compound BaFe2As2, on the other 
hand, has an in-plane stripe AF state at low temperatures 
(and with the ordered moment in the a&-plane instead of 
along the c-axis as in Fig. 1 for BaMn2As2),'' which in 
a local moment model requires Ji < 2J2 according to 
Eqs. (3). The in-plane stripe phase can be considered 
to consist of two interpenetrating G-type AF sublattices, 
where in this case a sublattice consists of all NNN spins 
of a given spin, and which are connected by an antiferro- 
magnetic interaction J2 (see the bottom panel of Fig. 2). 

IV. INELASTIC NEUTRON SCATTERING 
(INS) MEASUREMENTS AND ANALYSIS 

Figures 3(a) and 3(b) show images of the INS data 
taken at the base temperature of 8 K which share simi- 
lar features at each incident energy. Unpolarized inelas- 
tic neutron scattering contains contributions from both 
magnetic and phonon scattering. The magnetic scatter- 
ing intensity falls off with Q (or scattering angle 26) due 
to the magnetic form factor, while phonon scattering in- 
tensity increases like Q^. One can then observe a large 
contribution from magnetic scattering between 60 and 
80 meV, presumably due to spin wave excitations in the 
magnetically ordered phase, whose intensity only appears 
at small Q. On this intensity scale, strong phonon scat- 
tering is apparent below approximately 40 meV. 

This separation of magnetic and phonon scattering 
is more clearly shown by plots of the Q-dependence of 
the scattering averaged over different energy ranges, as 
shown in Fig. 4. For an energy range from 30 to 40 meV, 
the scattering is dominated by a large phonon contri- 
bution, whose intensity is proportional to Q^, and a 
large constant background due to multiple scattering and 
other background contributions. Q-dependent oscilla- 
tions arise from the powder averaging of the coherent 
phonon scattering and weak magnetic scattering. At the 
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FIG. 3: (Color online) Inelastic neutron scattering data from 
a powder sample of BaMn2As2 as measured on the Pharos 
spectrometer with incident energies a) Ei — 150 meV and b) 
Ei = 200 meV. The white lines delineate scattering angles of 
7 and 30° where the magnetic scattering was estimated. Pan- 
els c) and d) show calculations of the polycrystalline averaged 
spin wave scattering using a Heisenberg model. The calcula- 
tions in c) and d) are identical, however panel c) shows the 
trajectories of the angle summation limits for Ei — 150 meV 
and d) for Ei = 200 meV. 
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FIG. 4: (Color online) Q dependence of the inelastic neutron 
scattering intensity averaged over several energy ranges with 
incident energies a) Ei = 150 meV and b) Ei = 200 meV. The 
red lines show calculations of the polycrystalline averaged spin 
wave scattering using the Ji-J2-Jc Heisenberg model. 



higher energy ranges between 60 and 90 meV, the 
phonon contributions are gone and magnetic scattering 
appears superimposed on a constant background. The 
magnetic scattering intensity falls off with Q as the mag- 
netic form factor for the Mn^+ ion and is no longer visible 
above '-^ 6 A"-'^. Similar to the phonon cross-section, Q- 
dependent oscillations in the magnetic scattering occur 
due to coherent scattering of spin waves. 

The spin wave spectrum can be obtained by averag- 
ing the low-Q (low 29) data to improve statistics. How- 
ever, the magnetic scattering, especially below ~ 50 meV, 
must be separated from the phonon scattering and other 
background contributions. The pure phonon signal can 
be estimated from the high-angle spectra, where mag- 
netic scattering is absent. The magnetic scattering com- 
ponent can then be estimated by subtracting the high 



angle data (averaged from 26 = 50-90°) from low angle 
data (averaged from 29 — 7-30° and indicated by the 
white lines in Fig. 3) after scaling by a constant factor. 
These spectra are shown in Fig. 5 and show a strong 
and broad magnetic peak at ~ 70 meV. At energies 
below 50 meV, the subtraction of the phonon intensity 
is subject to error since the phonon intensity may not 
scale uniformly to low-Q due to coherent scattering ef- 
fects and also due to the different Debye- Waller factors 
for each atomic species. It is difficult to quantify this 
error without detailed phonon models; however, most of 
the magnetic scattering occurs above the phonon cut- 
off. Thus the errors introduced are only a problem below 
50 meV and the isolated magnetic data in this energy 
range can contain large errors. 
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FIG. 5: (Color online) Energy dependence of the inelastic 
neutron scattering intensity averaged over a scattering angle 
range from 7 to 30° with incident energies a) Ei — 150 meV 
and b) Ei — 200 meV. The magnetic intensity was extracted 
from the total scattering as described in the text. The red 
lines show calculations of the polycrystalline averaged spin 
wave scattering using the Ji-J2-Jc Heisenberg model. 



A. Spin Waves in the J\-J2-Jc Heisenbeerg Model 
for a G-type Antiferromagnet 

1. Spin Wave Theory 

In order to analyze the Q and w dependence of the 
magnetic spectra, we utilize a model of the spin wave 
scattering in BaMn2As2. Spin waves in insulators such 
as BaMn2As2 with the ThCr2As2 structure can be de- 
scribed by the Heisenberg Hamiltonian (2) except that 
here we set the magnetic field H in the last term to zero. 

The spin wave dispersions for the G-type AF struc- 
ture are obtained from a Holstein-Primakoff spin-wave 
expansion of the Heisenberg model. When the single-ion 
anisotropy is zero, the dispersions with respect to the 
body-centered-tetragonal (bet) lA/mmm unit cell con- 
taining two formula units of BaMn2As2 are 



hu) (q) 
2SJi 



Ji 



Ji 



[2 - cos(g^a) 



s{qya)] I 



I cos (qx+qy)^ +COS {qx ~ qy) 



Jl 



my 



(6) 



where a — 4.15 and c = 13.41 A are the lattice parame- 
ters of the bet unit cell at our measurement temperature 
of 8 K." 

In the absence of an anisotropy-induced energy gap in 
the spin-wave spectrum, the long-wavelength spin wave 
energies are described for an orthogonal (cubic, tetrago- 



nal, or orthorhombic) antiferromagnetically ordered spin 
lattice by the generic dispersion relation 



=. huj{q) = n^vlql + vlql + vlql, (7) 

where Wa, v\, and are the spin wave velocities (speeds) 
along the respective axes. In our case of tetragonal sym- 
metry we have 



(8) 



where Vah = Va ~ Vh- For G-type AF ordering of a spin 
lattice with our bet unit cell, these velocities are derived 
from the dispersion relation in Eq. (6) as 



(9) 




From the first of Eqs. (9) the in-plane spin wave velocity 
Vab decreases with increasing J2, consistent with expec- 
tation since according to Fig. 2, a positive (AF) J2 is a 
frustrating interaction for G-type AF ordering. Indeed, 
Vab vanishes when J2 = Ji/2, which is the classical crite- 
rion in Eq. (4) for the transition between the G-type and 
stripe-type in-plane AF ordering arrangements. 

In order to make contact with previous spin wave cal- 
culations for isotropic and anisotropic primitive orthog- 
onal spin lattices, one can change unit cell variables to 
those of the primitive tetragonal (pt) spin lattice contain- 
ing one spin at each lattice point with lattice parameters 
a' and c'. Referring to the bet structure with lattice pa- 
rameters a and c in Fig. 1, the pt spin lattice parameters 
are related to these according to 



a = V2a' 
c = 2c. 



(10) 



Furthermore the pt unit cell is rotated about the c-axis 
by 45° with respect to the bet unit cell, so the pt wave 
vectors q'^ , q'y and q'^ are related to those with respect to 
the bet cell qx,qy,qz by 



qx 
-qx 



qy 
qy 
qz 



^q'x 

V2q' 



(11) 



With these conversion expressions, the dispersion rela- 
tion in Eq. (6) becomes 



huj (q') 



2SJi 



{2 + ^-^[l-cos(q;a')cos(<z;a')]}' 

Jc 1 ^ 

cos(g' a') -I- cos(g' a') + cos(g^c') (12) 
Jl J 



2J2 
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Our dispersion relation (12) is identical to that in Refs. 21 
and 22 derived from linear spin wave theory for the J1-J2- 
Jc model. Also, Eq. (12) with J2 set to zero is identical to 
that in Eq. (5) of Ref. 26 and in Eq. (3) of Ref. 27 for the 
anistropic simple cubic G-type bipartite antiferromagnet. 

Using Eqs. (10), for the primitive tetragonal spin lat- 
tice the spin wave velocities in Eqs. (9) become 



S., = 2v^J.Se' + (13) 

In a simple cubic bipartite spin lattice with one spin 
per lattice point and isotropic interactions with c' = a', 
Jcl J\ = 1 and J2 = 0, the spin wave velocity is isotropic 
with magnitude hv' — 2\/iJ\Sa\ which is the same as 
given previously in Table I of Ref. 26 where a! was set to 
unity and is the standard well-known result when mag- 
nctocrystalline anisotropy is negligible. 

2. Application of Spin Wave Theory to BaMn2As2 

Spin wave dispersions using the bet notation in Eq. (6) 
are plotted in Fig. 6 (in units of SJi) for three different 
combinations of the exchange ratios J2/J1 and Jc/Ji- 
The notations in Fig. 6 and Table I for labeling the 
zone boundary reciprocal space positions are given by 
Kovalev.^^ The magnetic excitation wavevector q values 
are in reciprocal lattice units given by (i/, K, L) r.l.u. 
This is a shorthand for q expressed in inverse length units 
of the bet chemical unit cell according to 

2ttH 2nK - 27rL 
q = a H b H c. 

a a c 

In the lA/mmm bet unit cell notation, the magnetic 
propagation vector for G-type AF ordering is tq = 
(1,0,1), which gives tq — (0,0,0) when translated by 
a reciprocal lattice vector to the F-point in the Brillouin 
zone. This corresponds to the more familiar G-type wave 
vector Tq = (5,5,5) in the pt cell according to the 
transformation q' { ^t^ ^ , 7), where H,K,L are 
referred to the lA/mmm crystallographic unit cell, as 
shown in Eqs. (11). 

For J2 = and Jc/Ji — 1 (the top black curve in 
Fig. 6), the dispersion is that of an isotropic G-type anti- 
ferromagnet (similar to LaFeOa)^^ with a maximum spin 
wave energy of QSJi. For the layered BaMn2As2 struc- 
ture, Jc is expected to be much weaker than Ji. When 
J2 = and Jc/Ji =0.1 (the middle red curve in Fig. 6), 
the maximum spin wave energy is reduced to ^ 45Ji 
and the zone boundary M-point spin wave at q = (001) 
is strongly reduced. If we now turn on antiferromag- 
netic NNN interactions with J2/J1 = 0.25 (the bottom 
blue curve in Fig. 6), we observe a further softening of 




r SMRX rMS (101) N r 



Spin Wave Momentum Density-of-States 



FIG. 6: (Color online) The spin wave dispersion of the Ji- 
J2-Jc model along various symmetry directions for three dif- 
ferent combinations of the exchange parameters: J2/J1 = 
and Jc/Ji = 1 (top black curve); J2/J1 = and Jc/Ji =0.1 
(middle red curve); and J2/J1 = 0.25 and Jc/Ji =0.1 (bot- 
tom blue curve). Energies are reported in units of SJi. The F 
point in the Brillouin zone corresponds to wave vector q = 0. 
The q values corresponding to the high-symmetry M, X and 
N points are given in Table I. These wave vector directions 
are written with respect to the lattice translation vectors of 
the lA/mmm chemical unit cell. The labels E, F, U, Y, A 
and A correspond to high symmetry lines in reciprocal space. 
In the far right-hand panel, the energy versus the spin wave 
density-of-states is shown for J2/J1 =0 and Jc/Ji = 1 (top 
gray region) and J2/J1 = 0.25 and Jc/Ji = 0.1 (bottom blue 
region) . 

the spin wave spectrum, most notably at the X-point. 
When J2 > Ji/2, the G-type ordering becomes unstable 
and the stripe AF order is the new ground state with 
ordering wavevector at the X-point as discussed above 
in Sec. Ill, and our spin wave expressions are no longer 
applicable. The spin wave theory for the stripe phase in 
the Fe-based superconductor parent compounds with the 
ThCr2As2 structure is reviewed in Ref. 4. 

The powder-averaged spin wave scattering is closely as- 
sociated with the spin wave density-of-states [SWDOS, 
g{i^)]- The SWDOS is the distribution of spin wave 
energies and is determined by the summation over all 
wavevectors in the Brillouin zone (q), 

g(c.) = ^^[?ic.-?ic.(q)]. (14) 
q 

The SWDOSs versus energy Ulu are shown on the right- 
hand side of Fig. 6. It is observed that the SWDOS 
remains sharply peaked when J2 — 0, and that J2 acts 
to broaden the SWDOS. Table I indicates the energies of 
the various extremal features in the SWDOS (van Hove 
singularities) for ratios J2I J\ and Jd J\ that are associ- 
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TABLE I: Energies of various van Hove singularities in the 
spin wave density-of-states of G-type antiferromagnets with 
NN (Ji), NNN (J2) and interlaycr (Jc) exchange interactions. 



label 



van Hove singularity energy 



(001) 
(^,5.0) 

(io,i) 
f 1 1 1) 

\ 2 ' 2 ' 2 / 



Ex 
En 
Ep 

E' 



45 Ji. 

2SJi(2+i~2^) 
2SJi(2+i-4f) 



ated with zone boundary spin wave energies. 



B. Calculations of the Scattered Intensity 

When performing an INS experiment on a powder, the 
resulting INS intensities arise from the averaging of the 
inelastic scattering structure factor S{Q, oj) over all ori- 
entations of the crystallites. Despite the orientational 
averaging, the spectra can show evidence of the spin 
wave dispersions, especially at low angles (within the 
first Brillouin zone) and in the vicinity of the first few 
magnetic Bragg peaks. Due to the weighting of the spin 
wave modes by coherent scattering intensities, the Q- 
averaged intensity, S{uj), as shown in Fig. 5 does not 
necessarily give the SWDOS. This is only true in the 
incoherent scattering approximation, which does not ap- 
ply to the case of scattering from a magnetically ordered 
system. Therefore, model calculations of the powder- 
averaged spin wave intensities are necessary for accurate 
comparison to the data. 

Numerical calculations of the spin waves in the Heisen- 
berg model give not only the dispersion relation u;„ (q) for 
the n*^ (degenerate) branch [as shown in Eq. (6)], but 
also the spin wave eigenvectors, T„i(q), for the z*"'' spin 
in the magnetic unit cell. The dispersion and associated 
eigenvectors can be used to calculate the spin wave struc- 
ture factor for unpolarized neutron energy loss scattering 



from a single-crystal sample, 5mag(Qj'^)j given by 



'S'mag(Q>'^) 



1 + 



(A-Q)' 



(15) 



n i 

X [n{ijj) + 1] 5[u) — w„(q)], 

where the i^^ spin with magnitude Si pointed in direc- 
tion /i is located at position dj and ai = ±1 is the 
direction of the spin relative to the quantization axis 
fi for a collinear spin structure, as shown in the top 
two rows of Table II. The vector q = Q — r is the 
spin wave wavevector in the first Brillouin zone. Fi- 
nally, the function n(w) is the temperature-dependent 



Bose factor and Fi(Q) = 55i./i(Q)e~^''^'^^ is a product 
of the spectroscopic splitting factor (g- factor), magnetic 
form factor, and Debye-Waller factor for the i^"^ spin, 
respectively. The constant (7ro)^ = 290.6 millibarns 
allows calculations of the cross-section to be reported 
in absolute units of [millibarns steradian"^ meV~^ (for- 
mula unit)~^]. ForBaMn2As2, all Mn ions in the mag- 
netic cell are equivalent. The structure factor can then 
be written 



'S'mag(Q) w) 



1 



1 



(A-Q) 



x5^|$^aiT„,(q)e-^Q-l* 

n i 

x[n(w)-M](5[w-w„(q)]. 



2 



(16) 



In the calculations, we use the isotropic magnetic form 
factor for Mn found in the International Crystallogra- 
phy Tables'^" and the Debye-Waller factor is set equal 
to unity. The differential magnetic cross section that is 
measured in the inelastic neutron scattering experiments 
is proportional to the structure factor. 

To compare Heisenberg model spin wave results to the 
powder INS data, powder-averaging of SmagCQ, '^) is per- 
formed by Monte Carlo integration over 25 000 Q vectors 
lying on a constant-Q sphere, giving the orientationally 
averaged S'mag(Oi'^) which depends only on the magni- 
tude of Q. By a comparison of the total S{Q,uj) in Fig. 3, 
the Q-cuts in Fig. 4, and the energy spectra in Fig. 5, 
we arrive at the following parameters; SJi = 33 meV, 
SJ2 = 9.5 meV (J2/J1 = 0.29), and SJc = 3 meV 
[Jc/Ji = 0.09), as summarized in Table II. Figures 3(c) 
and 3(d) show that calculations of SmagiQ,!^) at 8 K us- 
ing these parameters compare well to the corresponding 
data in Figs. 3(a) and 3(b) and show clearly the coherent 
scattering of the powder-averaged spin waves. The most 
obvious coherent scattering feature is the necking down 
of acoustic spin waves in the vicinity of allowed magnetic 
Bragg reflections. The flrst two observed magnetic Bragg 
peaks are at Q = (101) and (103). Additional coherent 
scattering features can also be seen for zone boundary 
spin waves, where intensities tend to peak in between 
the allowed magnetic Bragg peaks. Figure 3 enforces the 
general agreement of the Heisenberg model calculations 
of the spin wave intensity with neutron scattering mea- 
surements. 

More quantitative estimates of the agreement of the 
calculated spin waves and the data are shown in Figs. 4 
and 5. The calculations can be summed over scattering 
angles in order to compare the equivalent angle-summed 
data, as shown in Fig. 5. The success of the Heisenberg 
model in estimating the measured spin wave intensities 
is better observed by plotting constant energy Q-cuts, 
as shown in Fig. 4. The plots show Q oscillations of 
the experimental magnetic spin wave scattering above a 
background due mainly to phonon scattering and back- 
ground/multiple scattering. A constant background and 
incoherent phonon scattering intensity (proportional to 



10 



TABLE II: Parameters of the Heisenberg model in Eq. (2) 
for BaMn2As2 determined from magnetic inelastic neutron 
scattering measurements at a temperature of 8 K. The Mn 
positions i in the top two rows refer to the crystallographic 
li/mmm unit cell with lattice parameters a = 4.15 and c = 
13.41 A at 8 K. The moment direction is along the c-axis. The 
spin S of the Mn atoms is not independently determined from 
the measurements. Only the products of 5* with the exchange 
constants Ji can be modeled. All exchange parameters are 
positive (antiferromagnetic). Also shown are the low-energy 
spin wave velocities in the afe-plane Vat and along the c-axis Vc, 
each multiplied by h, calculated from the exchange constants 
using Eqs. (9). 



i 






1 


+1 


(0,0,0) 


2 


-1 


(iio) 


exchange constant 


value 


value (K) 


SJi 


(33 ± 3) meV 


380 K 


Ji (S = 2, 5/2) 


16.5, 13.2 meV 


190, 150 K 


SJ2 


(9.5 ± 1.3) meV 


110 K 


J2 {S = 2, 5/2) 


4.8, 3.8 meV 


55, 44 K 


SJc 


(3.0 ± 0.6) meV 


35 K 


Jc {S = 2, 5/2) 


1.5, 1.2 meV 


18, 14 K 


5'(2Ji + Jc) 


69 meV 


800 K 


2Ji + Jc{S = 2, 5/2) 


18.0, 14.4 meV 


400, 320 K 


J2/J1 


0.29 ± 0.05 




Jc/Jl 


0.09 ± 0.02 




spin wave velocity 


value 
(meVA) 




hVab 


180 




hvc 


190 





Q^) are added to the calculated spin wave scattering in 
order to compare to the measured data. The agreement 
confirms the adequacy of the parameters. 

The low-energy spin wave velocities in the a6-plane 
Vab and along the c-axis Vc calculated from the exchange 
constants in Table II using Eqs. (9) are shown in Table II 
for our measurement temperature of 8 K. Remarkably, 
in spite of the layered nature of the spin lattice, the ab- 
plane and c-axis spin wave velocities are seen to have 
nearly the same value hvat ~ f^Vc ~ 180-190 meVA. 
For comparison, the spin wave velocities in the AFe2As2 
compounds are in the ranges hvat w 280-570 meVA and 
hvc « 57-280 meVA.4 



V. MAGNETIC SUSCEPTIBILITY 
MEASUREMENTS 

The anisotropic magnetic susceptibilites x(T) of a sin- 
gle crystal of BaMn2As2 in an applied magnetic field 
H = 3 T are shown in Fig. 7(a) for temperatures of 
300 to 1000 K, together with our previous data^^ below 
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FIG. 7: (Color online) (a) Magnetic susceptibility x versus 
temperature T of single crystals of BaMn2As2 with the ap- 
plied magnetic field parallel to the c-ajcis (xc) and to the ab- 
plane (Xab)- The individual symbols are the data previously 
reported in Ref. 11. The solid curves are the present data 
obtained in an applied magnetic field H = 3 T. The Neel 
temperature Tn is indicated, (b) Expanded plot of Xc(r) for 
temperatures around Tn. The temperature of the maximum 
slope of Xc{T) gives Tn = 618(3) K. 



350 K. Our Xab{T) data are consistent with the pre- 
vious Xab{T) data over the temperature range of over- 
lap (300-400 K),^^ but there is a difference between 
the c-axis data sets over that overlap temperature range 
for reasons that are not clear to us. The temperature 
of the maximum slope of Xc{T) from Fig. 7(b) gives 
the Neel temperature as Tn — 618(3) K, nearly the 
same as the value of 625(1) K determined from the pre- 
vious magnetic neutron diffraction measurements on a 
powder sample. ^"^ Above Tn, the susceptibility is nearly 
isotropic and exhibits negative curvature. The suscep- 
tibility appears to reach a maximum at a temperature 
jnmax ^ 1000 K, where the value of the average suscepti- 
bility is x]^T = X 10~^ cm^/mol and a "mol" refers 
to a mole of formula units (f.u.) unless otherwise stated. 
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TABLE III: Parameters describing the magnetic behaviors of 
BaMn2As2. Here, Tn is the Neel temperature, Xab is the 
magnetic susceptibihty with the magnetic field aligned in the 
a6-planc, Xc is the magnetic susceptibility with the magnetic 
field aligned along the c-axis, and Xave = (Sxab + Xc)/3 is 
the powder-averaged value of the susceptibility. Our x{T) 
results and those of Ref. 11 are on single crystals. The Tn 
value listed for Ref. 13 was obtained from magnetic neutron 
diffraction measurements on a powder sample. 



Property 


Value 


Reference 


Xab{W K) 


1.50(2) X 10"^ cm^/mol 


11 


Xc(10 K) 


0.20(2) X 10"* cm^/mol 


11 


Tn 


625(1) K 


13 


Tn 


618(3) K 


This work 


Xorb 


0.20(2) X 10"* cm*/mol 


This work 


Xab{TN) 


1.66 X 10"* cm*/mol 


This work 




1.60 X 10"* cm*/mol 


This work 


Xavc(rN) 


1.64 X 10"* cm*/mol 


This work 


Xspin(TN) 


1.44 X 10"* cm*/mol 


This work 


XabilOOO K) 


1.81 X 10"* cm*/mol 


This work 


Xc(1000 K) 


1.76 X 10"* cm*/mol 


This work 


Xave (1000 K) 


1.79 X 10"* cm*/mol 


This work 


XBpin(1000 K) 


1.59 X 10"* cm*/mol 


This work 



The values of the anisotropic susceptibihties at several 

distinctive temperatures are summarized in Table III. 

One can partition the measured susceptibility xC^) of ^ 
material into spin Xspin and orbital Xorb parts. Generally 
the orbital part is independent of T but Xspin does depend 
on T, so one obtains 

X{T) = Xorb + Xspin (T). (17) 

The Xorb generally consists of paramagnetic Van Vlcck 
and diamagnetic core contributions, plus the Landau 
diamagnetism of conduction electrons which is not sig- 
nificant in semiconducting BaMn2As2. From Fig. 7(a), 
the measured xC^ > "^n) is (nearly) isotropic. There- 
fore we infer that Xorb is isotropic at all T. For 
a collinear antiferromagnetic insulator (semiconductor) 
such as BaMn2As2, one expects the spin susceptibility 
parallel to the ordered moment direction, Xc spin in our 
case, to be zero at T ^ 0. From Fig. 7(a) we then obtain 

Xorb « 0.20(2) X 10"^ cm^/mol, (18) 

which we have included in Table III. Thus the spin sus- 
ceptibility is given by 

Xspin(T) = X(T) - Xorb. (19) 

We have listed the values of Xspin at T = Tn and T — 
1000 K in Table III. It appears from Fig. 7(a) that x{T) 
reaches a maximum at a temperature T™*" w 1000 K. 
Then one obtains from Table III the product 



Note that this value is for a mole of spins, not a mole of 
formula units. We will use this value later when compar- 
ing theory and experiment. 

The temperature dependence of x{T) above Tn in 
Fig. 7(a) is opposite to that expected for a fully three- 
dimensional antifcrromagnet, where x decreases rather 
than increases above Tn-^^ However, the behavior we 
observe above Tn is common in low-dimensional anti- 
fcrromagncts such as the tetragonal cuprate compound 
Sr2Cu02Cl2 where the intralayer magnetic coupling 
within the Cu+^ spin S = 1/2 square lattice is much 
stronger than the intcrlayer coupling.® Such antiferro- 
magnets exhibit a susceptibility with a broad maximum 
and the corresponding onset of strong short-range AF or- 
dering at a temperature T'^'^^ of order the mean-field AF 
long-range transition temperature [see Eq. (41) below]. 
However, for the compound Sr2Cu02Cl2 one estimates 
ymax ^ l^QQ K but it exhibits long-range AF ordering 
only at a much lower temperature Tn = 250 K <C T™*^. 
The interlayer coupling Jc is much smaller than the in- 
plane coupling Jab in quasi-two-dimensional antiferro- 
magnets. The suppression of Tn with respect to T™**^ 
is due to fluctuation effects associated with the low di- 
mensionality of the system. 

In the following we consider what can be learned about 
the signs and strengths of the exchange interactions in 
BaMn2As2 from analysis of our experimental data on 
this compound in terms of molecular field theory. Later 
in Sec. IX we develop the theory for fitting the experi- 
mental data taking into account the intralayer magnetic 
correlations that are present above Tn , which we will then 
apply to fit the x(T > Tn) data in Fig. 7(a) in Sec. X. 

VI. MOLECULAR FIELD THEORY (MFT) 

We will be analyzing various experimental data 
for BaMn2As2 using the Weiss molecular field theory 
(MFT). To introduce the MFT, we first consider the 
known results for a local magnetic moment model on 
a bipartite spin lattice with equal numbers of spins S 
in the two spin sublatticcs i and j interacting with the 
same nearest-neighbor (NN) exchange constant J with 
the Heisenberg Hamiltonian 

n = jJ^Si ■ Sj + giiBHY^St, (21) 

NN i 

where g is the spectroscopic splitting factor ((;- factor), 
/Ub is the Bohr magneton and H is the magnitude of the 
applied magnetic field which is in the ^;-direction. For 
such a quantum local moment system of identical spins 
interacting by NN interactions, if the susceptibility in the 
absence of J follows the Curie law xo = C/T, then in 
MFT the x(^) above the magnetic ordering temperature 
follows the Curie- Weiss law^^ 



(20) 



(22) 
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where the Curie constant C is 



1) 



3fcB 



(23) 



N is the number of spins and /cb is Boltzmann's constant. 
Taking N to be Avogadro's number A^a and g = 2 gives a 
useful expression for the Curie constant per mole of spins 

as 



C = 0.50020 ,S(,S+1) 
The Weiss temperature 9 is 



cm^K 
mol spins 



zJS{S+l) 



(24) 



(25) 



where z is the coordination number of each spin. Here, 
positive 9 corresponds to the case when J is positive (AF 
interactions), whereas a negative 9 corresponds to the 
case when J is negative (FM interactions). If 9 is positive, 
then the magnetic ordering temperature is Tn — 9 for AF 
ordering. On the other hand, if 9 is negative, then FM 
ordering occurs at the Curie temperature Tq = \9\. 

As discussed in Appendix A, the Curie- Weiss law is 
not simply a mean-field expression. ^'"'^"'^^ It arises from 
the first (l/T) term in the exact quantum mechani- 
cal high-temperature series expansion of the nearest- 
neighbor two-spin correlation function and is accurate 
in the limit that higher order 1/T" terms in the two- 
spin correlation functions are negligible. Thus the Curie- 
Weiss law, and hence our scaling expressions in Eqs. (82) 
and (87) below, begin to fail when l/T^ and higher or- 
der terms in the two-spin correlation functions become 
significant compared to the 1 /T term with decreasing T. 

Another important conclusion from Appendix A is that 
the Weiss temperature in the Curie- Weiss law results 
from all the spins that a given spin interacts with, ir- 
respective of the dimensionality of the spin lattice, of 
whether or not the spin lattice is bipartite (see Sec. VII) 
or whether all those interactions are the same, but where 
all spins are equivalent. Thus if there are different in- 
teractions present of a given spin i with other spins j, 
in Eq. (25) for the Weiss temperature one can make the 
replacement zJ — > '^y j where z is the total num- 

ber of spins that spin i has interactions with, giving the 
Weiss temperature as 



3fcB 



(26) 



VII. THE Ji-J'z-Jc HEISENBERG MODEL 
TREATED IN MOLECULAR FIELD THEORY 



respectively. Consider a specific spin i in sublattice 1. 
This spin i has four in-plane NN in sublattice 2 coupled 
by Ji and two out-of-planc NN in sublattice 2 coupled 
by Jo. Within the same sublattice 1, spin i is coupled 
to four in- plane NNN by J2. Since there are multiple 
exchange constants present from a given spin to its NN 
and NNN spins, we have 
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Jij = 2(2 Ji + Jc + 2 J2) 



and the Weiss temperature (26) becomes 

^ _ 2(2 Ji + Je + 2J2) S{S + 1) 

3fcB 



(27) 



We cannot measure 9 for BaMn2 As2 because according to 
Fig. 7 the temperature range required for the susceptibil- 
ity measurments to be in the Curie- Weiss regime would 
be far above 1000 K. 

In MFT, the magnetic induction B = Bk seen by each 
sublattice 1 and 2 is the sum of the applied field H = iJk 
and the respective exchange field H = ilexchk, i.e., 



Bl = H + Hi exch 
B2 = H + H2exch- 



(28) 



The MFT exchange field Hexch seen by each sublattice is 
respectively 



B^l exch 
H2 exch 



AsMi + AdM2 
AdMi + AsM2, 



(29) 



where As is the net molecular field coupling parameter 
for coupling within the same sublattice and Ad is the net 
molecular field coupling parameter for coupling between 
the two different sublatticos. We will obtain in Eq. (35) 
below expressions for these A values in terms of the J 
parameters in Hamiltonian (2). 

Wo only consider here the limit of low applied fields H. 
In MFT, the magnetization of each sublattice 1 and 2 is 
given by the response to the applied field plus the ex- 
change field as 



MiiT,H) 



M2{T,H) 



Xo{T)Bi 



2 



{H + AsMi + AdM2) (30) 



Xo{T)B2 
2 

Xo{T) 



{H + X,Mi + XdM2), 



The Hamiltonian (2) represents a situation where there 
is coupling both between the two spin sublattices and 
within each sublattice, where the two sublattices 1 and 
2 correspond to the red (up-pointing) and blue (down- 
pointing) magnetic moments in the top panel of Fig. 2, 



where Xq{T) = liniH^oM/H is the temperature- 
dependent spin susceptibility of the whole system in the 
absence of the explicit exchange fields, the factors of 1/2 
are there because each sublattice only has half of the total 
number of spins, and Mj is the z-axis magnetization of 
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the system induced by a magnetic field in the z-direction 
with magnitude H. In the paramagnetic state, M2 = Mi 
and Eqs. (30) yield 



Mi{H,T) 



'T)II/2 



l-(Xo/2)(Ad + As)' 

where i = 1,2. Since M ~ 2Mj, one obtains the spin 
susceptibility x{T) = 2Mi/H as 



X{T) = 



Xo{T) 



(xo/2)(Ad + As)- 



(31) 



The inverse susceptibility is 



1 



Ad + As 



X{T) Xo{T) 



This is typical of molecular field theory, where the molec- 
ular exchange field just shifts the inverse susceptibility up 
or down by a temperature-independent amount that de- 
pends on the sign and magnitude of the net moleciilar 
field coupling constant. It is important to note, with re- 
spect to fitting experimental data by molecular field the- 
ory, that the presence of molecular fields cannot change 
the temperature of peaks in the susceptibility Xo{T) that 
is assumed in the absence of explicit exchange couplings. 
For example, one could take Xo{T) to be the suscepti- 
bility of the isotropic square lattice Heisenberg antiferro- 
magnet such as in Fig. 17 below, which has a broad peak 
at T ~ J/fce- If one uses a molecular exchange field to 
magnetically couple the square lattice layers, this molec- 
ular field cannot change the temperature of the broad AF 
short-range ordering peak. 

To determine the magnetic ordering tempcrature(s) 
Tm, we set the applied field H to zero in Eqs. (30) and 
solve for nonzero Mi and M2. For the general case one 
obtains 



(As ± Ad) = 1, 



(33) 



so Tm depends on the assumed Xo{T)- From Eqs. (29), 
we see that for G-type AF ordering, we need to have Ad 
to be negative, so we take the minus sign in Eq. (33) to 

get 



Xo(Tn) 



(As - Ad) = 1, 



(34) 



where now is the antiferromagetic ordering (Neel) 
temperature Tn- Now we can use the solution for a A 
in terms of the related J vahie(s) from Ref. 31 to get 



As 
Ad 

which yield 

As — Ad = 



(4J2) 



^)(4J. + 2J,), (35) 



(4Ji + 2Je-4J2). 



(36) 



Inserting this expression into Eq. (34) for G-type antifer- 

romagnets gives 



Xo{Tn 



4Ji + 2.7, - 4J2 



= 1. 



(37) 



This is a constraint on the exchange parameters in 
BaMn2As2 in addition to those in Eqs. (4). If Xo{T) 
is the spin susceptibility per mole of spins, then A'' is 
Avogadro's number A'a- Taking g = 2 we have 



fee 



1.500 



mol 



(32) and Eq. (37) becomes 



4Ji -I- 2Jc - 4J2 1.500 cmVmol 



Xo(Tn) 



(38) 



(39) 



In the following sections we will assume that the spin 
susceptibility in the absence of any explicit exchange 
fields follows a Curie law, xo{T) = C/T. Then Eqs. (24) 
and (39) yield 



2Ji Jc — 2J2 



3Tn 



or 



2S{S+1)' 
2(2 Ji + Je - 2 J2) 



3fcF 



■S{S+1). 



Substituting Eq. (41) into (36) gives 



As — Ad 



Ng'^SiS + l)/i| 



2 ■ 



(40) 



(41) 



(42) 



It is useful to express differently how the NNN intra- 
sublattice interaction J2 affects Tn. From Eq. (41), one 
obtains 



Tn(J2) 
Tn(J2 = 0) 



2J2 



2Ji + Jc' 



(43) 



which is independent of the spin S and only depends 
on the ratio of the intrasublattice exchange constant J2 
to the net intersublattice exchange constant 2Ji -|- Jc. 
From Fig. 2 and Eq. (43), an antiferromagnetic J2 > 
is frustrating for G-type AF ordering and hence lowers 
Tn, whereas a ferromagnetic J2 < is nonfrustrating for 
G-type AF ordering and instead enhances Tn. 



A. Neel Temperature Reduction Factor / 

One can define a Neel temperature reduction factor / 
for antiferromagnets by 



/ = 



Tn' 



(44) 
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where 9 is the positive AF Weiss temperature in the 
Curie- Weiss law in Eq. (22). For molecular-field bipar- 
tite antiferromagnets with only nearest-neighbor inter- 
actions, / = 1.^^ However, there are four classes of AF 
materials in which / can be much different from unity: 
(1) materials in which fluctuation effects associated with 
a low-dimensionality (0, 1 or 2) of the spin lattice are 
strong, (2) three-dimensional materials in which geomet- 
ric frustration for AF ordering occurs, (3) spin lattices 
in which the signs of the exchange interactions of a spin 
with its neighbors frustrate the ordering, and/or (4) spin 
lattices that are not bipartite; i.e., interactions between 
spins on the same sublattice occur. In each of these 
classes of materials, Tn can be strongly suppressed, some- 
times to T = 0, which gives / » 1. Alternatively, it can 
occur that second neighbor interactions can enhance Tn 
but suppress 1^1 as we will see below in Eq. (45) if J2 
is negative (ferromagnetic). It can occur that a given 
compound belongs to more than one class. 

One of us has discussed class (1) in the context of 
low- dimensional copper oxide compounds such as quasi- 
two-dimensional La2Cu04 containing a Cu"*"^ (fi spin- 
1/2 square lattice and quasi-onc-dimcnsional Sr2Cu03 
containing Cu+^ spin-1/2 chains.^ In these mate- 
rials the AF correlation length ^ grows with decreas- 
ing T. In La2Cu04, long-range AF ordering occius at 
Tn ~ 7r($/a)^ Jc/A:b, where the number of spins within 
an AF correlated area in the plane is ~ 7r(^/a)^, 
Jc is the interplanc nearest-neighbor exchange coupling 
constant and a is the square lattice parameter. A large 
number of spins within a correlated area amplifies the 
effect of a small Jc- In Sr2Cu03, grows much more 
slowly with decreasing T than in La2Cu04 because what 
is relevant here is the number of spins within a correla- 
tion length rather than within a correlation area, and the 
former is much smaller than the latter at the same tem- 
perature. Hence, one expects / for Sr2Cu03 to be much 
larger than for La2Cu04, as observed. The ^ and the 
Weiss temperature 6 are determined by the in-chain or 
in-plane exchange coupling Jab ^ Jc, respectively, and 
hence 't^ for both compounds. 

Ramirez has extensively discussed class (2).^^ In frus- 
trated three-dimensional antiferromagnets, the suscepti- 
bility follows a Curie- Weiss-like temperature dependence 
down to temperatures much less than 9. One can describe 
the physics in two equivalent ways. In one view, the AF 
correlation length ^ does not grow as fast as one would 
predict from the Curie-Weiss law where one expects ^ to 
diverge at the mean- field Tn = 9- An alternate equivalent 
explanation is that because the Curie- Weiss law holds to 
low temperatures T 9, which results in / ^ 1, the co- 
efficients of the higher-order 1/T" terms (n > 1) in the 
high temperature series expansions of the two-spin cor- 
relation functions in Eqs. (A4) and (A5) in Appendix A 
are strongly suppressed in frustrated antiferromagnets. 
BaMn2As2 likely belongs to classes (1), (3) and (4). 

Using Eqs. (27) and (41) which assume xq = C/T and 
Ji, Jc > 0, the ratio of the Weiss temperature to the Neel 



temperature for G-type antiferromagnets in the Ji- J2- Jc 
model within MET is 



^ Tn 
which gives 



2 Ji + Jc + 2J2 
2</i -\- Jc — 2^/2 



2J2 



1 



2.72 

2.h+.Ja 



2.J2 

2Ji + Jc 



^ /-I 

2Ji + Jc f + l 



(45) 



(46) 



Thus / depends on the sign and magnitude of the NNN 
in-plane interaction J2. For an antiferromagnetic J2 > 0, 
one gets / > 1, whereas for a ferromagnetic J2 < 0, 
one gets / < 1. The constraint on J2 in Eqs. (4) that 
J2 < Ji/2 still applies, giving an upper limit (for which 
Jc = 0) of 

/ < 3. (for G type AF ordering) (47) 
Using Eq. (46), one can rewrite Eq. (43) as 
Tn(J2) 2 



Tn(J2 = o) f + r 



(48) 



In MET in the paramagnetic state, the spin suscepti- 
bility (22) follows the Curie- Weiss law x{T) = C/{T+9), 
and x(r) reaches a maximum at T = Tn- Therefore we 
obtain the product 



maxrriinax 
Aspin 



Xspin 

(Tn)Tn 



c 



Tn+i 



■Tn 



Tn(1 + /) 1 + / 

0.5002 5(5 -h 1) cm^K 



(49) 



I + f mol spins 

where we used Eq. (24) in the last equality. This gives 
(0.5002 cm3 K/molspins)5(5 -hi) ^ 

Aspin 

From Eqs. (40) and (49) one can solve for 2Ji + Jc and 
J2 to obtain 



2Ji + Jc (3 cm3 K /mol)TN 



o.^max'T'n 
Aspin 



(51) 



-spin 

J2 2Ji -|- Jc 3Tn 



kn 2kB 45(5 + 1) 



Additional useful expressions include the following. 
From Eqs. (27), (41) and (44) one obtains 



2 Ji + Jc 



3fcBTN(/ + l) 



45(5 -hi) 
Then from Eq. (46) one gets 

3fcBTN(/ - 1) 



2J2 = 



45(5 -hi) 



(52) 



(53) 
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Now using Eqs. (52) and (53) we can rewrite the molec- 
ular field coupling constants in Eqs. (35) in the simple 
symmetric forms 



As = - 
Ad = - 



W-i) 

c 

TnU + 1) 
C ' 



(54) 



where C is the Curie constant in Eq. (23). 



B. Anisotropic xiT) below Tn 

We would like to compare our experimental anisotropic 
x{T) data below Tn with the MFT predictions using 
the Ji-J2-Jc model. We discuss first the perpendicu- 
lar susceptibility x_l and then the parallel susceptibility 
X||, where x± refers to the susceptibility with the ap- 
plied magnetic field perpendicular to the easy axis of the 
coUinear antiferromagnetic structure and X|| to the sus- 
ceptibility when the applied magnetic field is parallel to 
it. For BaMn2As2, X|| = Xc and x± = Xab- In the Heisen- 
berg model, above Tn the susceptibility is isotropic and 
hence X|| = X-L- Below Tn, x± and X|| are no longer the 
same. 

Below Tn of a coUinear antiferromagnet, one always 
has x\\ < X± (sec also Fig. 13 below). Since the magnetic 
energy of the system at low fields is — (l/2)x-f^^, if the 
field is aligned along the ordered moment axis the spin 
system can lower its energy via a "spin-flop" transition 
where the ordered moment axis rotates to be perpendic- 
ular to the applied field. To prevent this from happening, 
one needs to have an anisotropy energy present that is not 
included in the Heisenberg Hamiltonian. Otherwise one 
could never measure X|| • An important example of such 
an anisotropy energy is the axial single ion anisotropy 
energy with the form DS^ (for 5 > 1/2) and T> < 0, 
and/or higher order forms, that arise from the spin-orbit 
interaction of the magnetic moments with the crystalline 
electric field of the lattice. Here we assume that an in- 
finitesimal axial anisotropy is present with sufficient mag- 
nitude to prevent the ordered moment axis from flopping 
from the parallel to the perpendicular orientation when 
we are measuring the parallel magnetization in the limit 
of an infinitesmal field. We will not further consider the 
spin-flop transition in this paper. 

The X-l(T) and X||(^) are derived for the Ji-J2-Jc 
model at T < Tn in Appendix B. For the perpendic- 
ular susceptibility, one obtains the constant value 



X± 



1 



c 



c 



Tn(1 + /) Tn + ^ 



x(7^n), (T<rN), 

(55) 

using / This result is similar to that for a 

bipartite lattice, "^^ except in that case one has = Tn 
whereas in our case we have 6 = /Tn with, in general, 
/ 7^ 1 from Eq. (45). The estimated values of / from 
Eqs. (60) below are ~ 3-5 in BaMn2As2, i.e., Tn is much 



smaller than 9, but within MFT the susceptibility still 
follows the Curie- Weiss law x = C/{T + 6) down to Tn. 
This interesting behavior is the result of bond frustration 
for AF ordering (the antiferromagnetic NNN interaction 
J2 frustrates the occurrence of G-type AF ordering) and 
has been noted as a property of geometrically frustrated 
antiferromagnets . 

The dependence of X||(T)/x||(7n) on t = T/T^ deter- 
mined by solving Eqs. (B13), (B14) and (B18) is shown 
in Figs. 8(a) and 8(b) for spins S = 1/2 and S = 5/2, 
respectively, for various values of / = ^/Tn- The value 
/ = 1 corresponds to the conventional nonfrustrated bi- 
partite stacked square spin lattice as in the top panel of 
Fig. 2 with J2 = 0. Figure 8 shows that the presence of 
a nonzero diagonal coupling J2 has a strong influence on 
X\\{T)/x\\{Tn)- Complementary plots of X|| (T)/x|| (Tn) 
versus T/Tn at flxed / = 0, 1 and 3 for S = 1/2, 5/2 
and 10 are shown in Fig. 9. 



C. Ordered Moment versus Temperature below Tn 

The ordered moment in the antiferromagnetic state 
of BaMn2As2, which is the staggered moment n\ in 
Eq. (B16), has been previously measured, but not 
modeled. In Appendix C we determine the MFT pre- 
dictions on the basis of the Ji-J^-Jc model. In Fig. 10 are 
plotted the solutions of Eq. (CI) for the nonzero ordered 
moment /i^ = /ij versus reduced temperature T/Tn for 
classical spins and for four values of quantum spins. In 
contrast to quantum spins for which n\ approaches the 
respective saturation moment /igat = gSfiB exponentially 
fast for T due to an energy gap between the ground 
state and the lowest excited states, the low-temperature 
classical behavior is linear. This results in a magnetic 
heat capacity Cmag constant ^ as T ^ for classi- 
cal spins, which violates the third law of thermodynam- 
ics, whereas for quantum spins Cmag — >■ as T — >^ (see 
Fig. 11 below). 

Interestingly, the parameter / = that character- 

izes the influence of J2 on the magnetism has disappeared 
from the expression for A*|(T) in Eq. (B16) when the tem- 
perature scale is normalized by Tn. Thus Eq. (CI) and 
the plots in Fig. 10 are identical to the corresponding 
MFT predictions for an AF bipartite spin lattice with 
J2 = 0. However, in our case with J2 7^ 0, we must keep 
in mind that J2 has already manifested its influence on 
the magnetism by changing Tn- 



D. Zero-Field Magnetic Heat Capacity Cmag and 
Entropy 5mag below Tn 

The zero-field magnetic heat capacity Cmag(T') is de- 
rived in MFT in Appendix D as 



R ~ S + dt ' 



(56) 
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FIG. 8: (Color online) Parallel susceptibility x\\ versus tem- 
perature T through the Neel temperature Tn for G-type AF 
ordering using the Ji-J2-Jc model in molecular field theory 
(MFT) for various values of f — as listed, for spins (a) 

5 = 1/2 and (b) S = 5/2. The order of the curves from top 
to bottom is the same as in the figure legends. At tempera- 
tures T > Tn, X is isotropic. For T < Tn, the perpendicular 
susceptibility is constant, x± = X^C^^n) (not shown). The 
G-type AF state is unstable against the stripe AF state for 
/>3. 



where t — T/T^ is the reduced temperature and fil = 
fil/fisat is the reduced ordered (staggered) moment. 
Since ftl does not explicitly depend on J2 as discussed 
above, neither does Cg,ve{T/T^), but rather implicitly 
via the dependence of Tn on J2. The fll{t) is deter- 
mined by numerically solving Eq. (CI). Inserting this 
result into (56), Cinag(T) was calculated for several spin 
values as plotted in Fig. 11. One observes a triangular 
shape for Cmag(i) near Tn for each S, which is charac- 
teristic of the mean field solution, with a discontinuous 
increase ("jump") in Cmag(r) upon decreasing T through 
Tn given by^^ 



ACniagCTN) 

R 



(25-^1)^-1 
(25+1)2 + 1 



(57) 




— s = 


10 
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FIG. 9: (Color online) Parallel susceptibility x\\ versus tem- 
perature T through the Neel temperature Tn for the Ji- J2- Jc 
model in molecular field theory (MFT) for spins S = 1/2, 5/2 
and 10 and f = O/Tt^ values of (a) 0, (b) 1 and (c) 3. The 
order of the curves in each panel from top to bottom is the 
same as in the figure legends. The value / = 1 corresponds 
to the conventional bipartite lattice with J2 — 0. 



upon varying the spin S. From Eq. (57) one obtains 

V 



ACmag(7N) _ 3 

2 

5 

" 2' 



R 

ACmagCTk) 



R 



S = 



{S = 00) 



where 2S + 1 is the Zeeman degeneracy in zero field for 
a spin S. There is not a large range of ACmag possible 



consistent with Fig. 11. 

The evolution in Fig. 11 of the low temperature 
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FIG. 10: (Color online) Ordered moment fi^ = U'l versus tem- 
perature T from molecular field theory of the Ji- J2- Jc model 
for a collinear antiferromagnet for classical spins and for sev- 
eral quantum spins as listed, where the saturation moment is 
Msat = gSfiB- The order of the curves from top to bottom is 
the same as in the figure legend. Remarkably, the results are 
independent of J2. 
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FIG. 11: (Color online) Magnetic component of the heat ca- 
pacity Cmag, divided by the molar gas constant K, versus the 
ratio of the temperature T to the Neel temperature Tn ac- 
cording to the molecular field theory Eq. (56) of the J\-Ji-Jc 
model for a collinear antiferromagnet for classical spins and 
for several quantum spins as listed. The order of the curves 
from top to bottom is the same as in the figure legend. As in 
Fig. 10, the results are independent of J2. 



CmagCT") with increasing spin 5* is interesting. It develops 
a hump at a temperature that decreases vifith increasing 
5, until in the classical limit 5 — > 00 the hump merges 
into the classical finite- value behavior for T — >■ 0. The 
hump is required in order that the entropy of the disor- 
dered spin system increase vifith increasing S (see below) , 
since C,nag (T) is bounded from above by the classical pre- 
diction. For quantum spins, the heat capacity approaches 
zero exponentially at sufiiciently low temperatures irre- 
spective of the (finite) spin value, whereas for classical 
spins the heat capacity approaches a nonzero finite value 
for T -> 0. 



The magnetic entropy S'mag 
magnetic heat capacity via 



is determined from the 



*^mag (/) 

'r 



(58) 



The magnetic entropy obtained from Eq. (58) and from 
the data in Fig. 11 is plotted versus temperature for quan- 
tum spins in Fig. 12. The constant values for T > Tn 
as indicated by the notations on the right-hand ordinate 
agree with the values expected for disordered spins given 
by the molar magnetic entropy S'mag = R\n{2S + 1). 
For classical spins the calculated entropy for T is 
S'mag(r ^ 0+)/i? = limT^oln[r/(0 K)] = 00, which 
violates the third law of thermodynamics. 
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FIG. 12: (Color online) Magnetic entropy Smag/i? for quan- 
tum spins versus reduced temperature T/Tk according to 
molecular field theory for the quantum spins S indicated. The 
order of the curves from top to bottom is the same as in the 
figure legend. In the disordered state at T > Tn, the mag- 
netic entropy is the constant value S'mag ~ R\n{2S + 1) for 
each S, as indicated along the right-hand ordinate. 
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TABLE IV: Parameters determined from a fit of magnetic sus- 
ceptibility data by molecular field theory of G-type antiferro- 
magnetic BaMn2As2 with NN (Ji), NNN (J2) and interlayer 



2.0 



(t/c) 6xcli£Liig6 intGrcLctions. 







c 

CO 


1 

1 


Quantity 


5 = 2 


5 = 5/2 


/ = e/T^ 


2.75 


4.47 


E 



1 


(2Ji + - 2J2)/feB (K) 


156 


107 


CO 

1 




(2Ji + Jc)/fcB (K) 


293 


293 







2Ji + Jc (meV) 


25.2 


25.2 







S{2.h + Jc)/fcB (K) 


586 


733 






S{2Ji + Jc) (meV) 


50.5 


63.2 






J2/fcB (K) 


68 


93 







J2 (meV) 


5.9 


8.0 






SJz/fcB (K) 


136 


233 






572 (meV) 


11.7 


20.1 







VIII. COMPARISON OF THEORETICAL 
PREDICTIONS WITH EXPERIMENTAL DATA 
FOR BaMn2As2 

A. Comparisons with Molecular Field Theory 

We expect the Mn+^ ion in BaMn2As2 to have 
the high-spin configuration with spin S = 5/2. On 
the other hand, the observed ordered moment is /i = 
3.9(1) /xe/Mn/'^ suggesting from the relation ^ — gSfiB 
with g = 2 that S — 2. Therefore in the following we will 
consider both of these possibilities. 



1. Neel Temperature 

Using Eq. (40) and Tn = 625 K for BaMuaAsa, one 
obtains 



2J\ -\- Jc — 2J2 



= 156 K (5 = 2) 

= 107 K. (5* = 5/2) (59) 
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FIG. 13: (Color online) Comparison of the MPT predic- 
tion below Tn of the anisotropic susceptibilities in Eqs. (55) 
and (B18) versus temperature with the experimental data 
from Fig. 7(a). 



z = 4, J/fcB = 1600 K,S^ 1/2, and Tn = 325 K,^ which 
yields a Weiss temperature 6 — 1600 K and / = 4.9, and 
the magnetism of this compound is known to be described 
very well by two-dimensional physics in the temperature 
range above Tn.^ As a further comparison, the quasi- 
one-dimensional spin- 1/2 chain compound Sr2Cu03 has 
z = 2, J/kB = 2200 K, 5* = 1/2, and Tn = 5.4 K,^ which 
yields = 1100 K and / = 200. This large / value is the 
reason that Sr2Cu03 is often considered to be a model 
quasi-one-dimensional Heisenberg antiferromagnet.^ 

Again using the experimental X^f^T™^'' value from 
Eq. (20), Eqs. (51) yield 



2 Ji + Jc 



fee 



= 293 K (5 = 2) 

= 293 K (5 = 5/2) 

= 68 K (5 = 2) 

= 93 K. (5 = 5/2) 



(61) 



2. Magnetic Susceptibility 

Inserting the experimental Xsp^n"^"^^^ value from 
Eq. (20) into (50) gives 

/ = 2.75 (5 = 2) (60) 
= 4.47. (5 = 5/2) 

According to Eq. (47), the value of / for 5 = 5/2 is not 
possible for G-type AF ordering and hence 5 = 5/2 is 
ruled out by this criterion. The / value for 5 = 2 suggests 
that interlayer coupling might have a significant effect on 
the observed magnetic susceptibility above Tn. On the 
other hand, for the layered cuprate La2Cu04 one has 



The above results, summarized in Table IV, are only ap- 
proximate qualitative constraints on the exchange pa- 
rameters in BaMn2As2, because they assume that the 
susceptibility follows the Curie- Weiss law above Tn, 
which Fig. 7 shows is not accurate. In particular, if 
Jc/Ji ~ 0.1 as determined from the neutron scattering 
results and the theoretical predictions in Sec. XII, one 
obtains unrealistically large J2/J1 = 0.50 and 0.67 for 
5 = 2 and 5 = 5/2, respectively. The problem stems 
from the fact that T"^'^^ and Tn do not coincide, which 
is an inconsistency in the analysis. 

A comparison of the MFT predictions below Tn of 
the anisotropic susceptibilities in Eqs. (55) and (B18) 
with the experimental data from Fig. 7(a) is shown in 
Fig. 13. For the MFT dashed-line predictions we used 
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FIG. 14: (Color online) Ordered moment versus temper- 
ature T measured for BaMn2As2 from Ref. 13 (filled red cir- 
cles). The Neel temperature Tn and saturation moment /isat 
are given in the figure. Also shown are the molecular field 
theory (MFT) predictions for quantum spins S — 1/2, 2 and 
5/2 from Eq. (CI) (solid and dashed curves), where the order 
of the curves from top to bottom is the same as in the figure 
legend. 



Eq. (17) withxorb = 0.2x10"^ cm^/mol and Xspin(7N) = 
1.35 X 10-3 cmVmol. We used the value Tn = 625 K 
and the MFT parameter / = 2.75 hsted in Table IV 
for S — 2. The temperature dependences of the MFT 
predictions for the anisotropic susceptibilities are seen to 
be in semiquantitative agreement with the experimental 
data. We do not consider the case S = 5/2 because the 
large / = 4.47 > 3 in Eqs. (59) and Table IV for 5 = 5/2 
makes the G-type AF structure unstable with respect to 
the stripe AF structure in Fig. 2. 



3. Ordered Moment 

The theoretical MFT results for the ordered moment 
versus temperature in Fig. 10 for S = 3/2 to S = 5/2 
are nearly the same, so we do not expect to be able to 
differentiate between the two possibilities of S = 2 and 
S = 5/2 for the Mn spins in BaMn2As2 on the basis 
of the observed temperature dependence of the ordered 
moment. This expectation is confirmed in Fig. 14 where 
we compare the MFT predictions for S = 1/2, 2 and 5/2 
from Eq. (CI) with the experimental data from magnetic 
neutron diffraction measurements in Ref. 13. Although 
the overall temperature dependence of the data agrees 
with MFT, the data are not quantitatively fitted by the 
prediction for any particular fixed S value. 
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FIG. 15: (Color online) Heat capacity C versus tempera- 
ture T. The experimental Cp{T) data obtained for a single 
crystal of BaMn2As2 (Ref. 11) are shown as the filled black 
circles. The lattice heat capacity for a Debye temperature 
Od ~ 300 K is plotted versus T as the solid red curve. The 
sum of the lattice heat capacity and the magnetic heat capac- 
ity for spin S = 2 is shown as the dashed blue curve. 



4- High- Temperature Magnetic Heat Capacity 

Here we will compare our experimental heat capacity 
Cp data for BaMn2As2 single crystals at temperatures up 
to 350 K with the prediction of MFT for the magnetic 
heat capacity Cmag at high temperatures, i.e., near room 
temperature. To do this we will need to estimate the 
lattice heat capacity contribution using the Debye model. 

The heat capacity at constant pressure Cp{T) for a 
single crystal of BaMn2As2, previously reported by Singh 
et al.,^^ is shown in Fig. 15 for the measured temperature 
range 2-350 K. We fitted the data by the Debye function 
for the molar lattice heat capacity of acoustic phonons at 
constant volume, given by'^^ 

Cocbyc - 9ni? j .^——-^dx, (62) 

where n is the number of atoms per formula unit (n — 5 
here) for various values of the Debye temperature Qd. 
In order that CDcbyo(T) does not lie above the experi- 
mental data over any temperature range, the minimum 
value of 0D is about 300 K, for which the Debye function 
is plotted as the solid red curve in Fig. 15. The &b is 
evidently temperature-dependent because the deviation 
of the curve from the experimental data varies nonmono- 
tonically with temperature. From the same set of exper- 
imental Cp(T) data,^^ at low temperatures T < 5 K a 
value 8d = 246(4) K was deduced using the Debye 
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law [the low-temperature limit of Eq. (62)] given by 



CDebye = 0T>T 



The experimental data at the highest temperatures lie 
above the lattice heat capacity curve for Gd = 300 K 
in Fig. 15, suggesting the presence of one or more heat 
capacity contributions in addition to that due to acoustic 
phonons. 

We calculated the difference Cp — Cy for the lat- 
tice heat capacity for the compound BaFe2As2, where 
Cy is the lattice heat capacity at constant volume, 
according to the thermodynamic relation Cp — Cy = 
VmPy{T)B{T)T, where Vm is the molar volume, /?v 
is the volume thermal expansion coefficient, and B 
is the bulk modulus. For the 200-300 K tempera- 
ture range, using the values jiy « 4.8 x 10~^ K^^,'^^ 
B = 6.6 X 10^2 dyne/cm2,38 and Vm = 61.5 cm^/mol,'* 
we obtained Cp - Cy = (9.3 mJ/molK2)r. This gives 
(Cp - Cy)(300K) = 2.8 J/molK, which is about a factor 
of two too small to account for the difference between 
the data and the Dcbye curve. It was not possible to 
calculate a value of Cp — Cy specific to BaMn2As2 be- 
cause Pv{T) and B{T) have not been measured for this 
compound. 

The magnetic contribution Cmag(T') to the heat capac- 
ity at high temperatures was calculated using the MFT 
prediction in Eq. (56). We chose to calculate it for spin 
5 = 2 because the 5 = 5/2 possibility was ruled out by 
the large value of / > 3 for spin 5 = 5/2 in Eq. (60). 
Using Tn = 625 K, the C„iag(T') from two moles of 
spins 5 = 2 per mole of BaMn2As2 was added to the 
Debye heat capacity and is plotted as the dashed blue 
curve in Fig. 15. Now the calculated curve lies above the 
experimental data around room temperature, indicating 
that the magnetic heat capacity is smaller than predicted 
by MFT. The x(T) data in Fig. 7(a) appear to be ap- 
proaching a maximum at a temperature T"^^^ w 1000 K 
that is far above Tn = 625 K, indicating the occurrence 
of strong short-range AF ordering above Tn (see also 
Sec. IX below) . This removes spin entropy and decreases 
Cmag below the value expected from MFT at tempera- 
tures below Tn- This may be the reason for the sup- 
pression of Cniag(^) in our measurements around room 
temperature. 



B. Comparison of Experiment with Spin Wave 
Heat Capacity Theory at Low Temperatures in the 
Ji-Ji-Jc Model 

1. Theory 

The lack of significant susceptibility anisotropy above 
Tn in Fig. 7 indicates that single-ion anisotropy is small. 
This anisotropy, if present, gives rise to an energy gap in 
the spin wave excitation spectrum. Here we assume that 
the anisotropy gap is infinitesmally small and calculate 
the low-temperature magnetic heat capacity of AF spin 
waves in the Ji-J^-Jc model. This is an extension of the 
standard treatment for simple cubic spin lattices with 
isotropic NN exchange interactions. 

The original 1952 papers by Anderson"^^ and by 
Kubo^° give a clear prescription of how to do this us- 
ing a spin wave model with two AF sublattices 1 and 2 
containing a total of N spins 5. Their starting Heisen- 
berg Hamiltonian is 

U = jY^Si-Sj (64) 

m 

where there is only a single J and the sum is over dis- 
tinct nearest-neighbor spin pairs. In zero field and in the 
absence of significant anisotropy the diagonalized spin 
wave Hamiltonian contains the following term involving 
the excitation energies of spin waves 

S = ^ (niqTiWiq -I- n2q?iW2q), (65) 

q 

where q is the wave vector of a spin wave excitation, 
riiq is the occupation number of the mode for sublattice 
i, and the two terms correspond to excitations on the 
two degenerate spin wave branches wiq and Ci;2q associ- 
ated with the two spin sublattices, respectively. Since 
riiqTiwiq = n2qhLU2q = nqfiwq are degenerate, the excita- 
tion energy of the system can be written 

E = 2^nqliojq. (66) 
q 

The thermal-average energy of the spin waves is then 

^ave = 2E ,.Jr._l - (67) 
q 

where (riq) = l/(e^'^o/'=BT _ i) ig the Planck distribution 
function for the thermal-average number of quanta in an 

oscillator at energy hujq. One converts the sum into an 
integral over q for a three-dimensional spin lattice via 

E^f^/ / / '^q. (68) 

where Vs-pin is the volume per spin. The factor of N/2 
arises because each spin sublattice has A^/2 spins. Then 
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FIG. 16: (Color online) Spin wave dispersion relation 
huj/iJiS of the isotropic two-dimensional square lattice over 
the Brillouin zone of the primitive tetragonal space lattice. 
The dispersion relation is doubly degenerate everywhere. At 
low temperatures, there are two distinct doubly degenerate 
spin wave branches that are relevant, one at the F point at 
(0,0) and the other at (f , (and equivalent points). 



Eq. (67) becomes 



(27r)3 



dqx / dqy 

-Tv/a J—Tv/b 



IT / C 



dq2 



(69) 



Note that the integration in Eq. (69) is over the entire 
Brillouin zone of the primitive direct lattice (containing 
a single spin), not over the Brillouin zone of the mag- 
netic lattice. The reason for this important fact is that 
integrating over the Brillouin zone of a primitive space 
lattice with one spin in the basis sums up the response 
of a single spin, whereas if one were to integrate over 
an antiferromagnetic Brillouin zone, this zone would in- 
clude the response of more than one spin. Indeed, the 
average energy per spin calculated this way does not de- 
pend on the type of magnetic ordering at all, even if 
the magnetic ordering is ferromagnetic or incommensu- 
rate. The only relevant difference between the thermal 
average energy per spin of different magnetic ordering 
configurations is the difference between the specific coq 
functions and their degeneracies over the Brillouin zone 
of the primitive space lattice. 

The dispersion relation for a general spin lattice is 



hujq = zJSJl - 7^ 



(70) 



make a point that will be illustrated using the spin wave 
spectrum of an isotropic two-dimensional square spin-S* 
lattice {z — A). In this case Eq. (71) yields 



1 



7q 



cos(fcxa) + cos(fcj,a)] 



and Eq. (70) gives the doubly degenerate dispersion re- 
lation as 



hujq = AJSJ 1 — [cos(A;j;a) + cos{kya)]^ / 4. 



(72) 



This dispersion relation is plotted in Fig. 16. One sees 
that uiq has doubly degenerate branches arising from zero 
energy at the F point (0,0), as expected, but also at the 
corners of the Brillouin zone at (-^j f ) and equivalent 
points. In a three-dimensional spin lattice with Jc ^ 0, 
using the dispersion relation in Eq. (12), one sees that the 
low-energy points of the dispersion relation move from 
the (f 7 f , 0) points in the corners of the two-dimensional 
Brillouin zone to the (f : -f j f-) and equivalent points at 
the other four corners of the three-dimensional Brillouin 
zone. Thus in either case there is another multiplicative 
factor of two to include in Eq. (69) if we only integrate 
over the two degenerate F point branches for T — > 0. 

Equation (69) is evaluated in Appendix E to yield the 
magnetic heat capacity per mole of spins at low temper- 
atures due to the spin waves as 



4^^i?fc|Kpin 



rj-io 



(T«rN) (73) 



where R is the molar gas constant, Kipin is the volume 
per spin, and Vx, Vy, Vz are the spin wave velocities along 
the a-, h- and c-axes, respectively. This expression in- 
cludes the contribution of the low energy spin waves at 
the Brillouin zone corners, and can be written in a form 
analogous to Eq. (63) for phonons as 



Cmag — kwT^, 

Psw — ^ "7771 



(74) 



By writing the Debye temperature in Eqs. (63) in terms 
of its constituent quantities, "^^ one obtains the lattice heat 
capacity coefficient /3d per mole of atoms as 



27r2flfc|Ktom 



(75) 



where 



7q 



1 ^ 

z ^ 



(71) 



z is the coordination number of a spin on one sublat- 
tice by spins on the other sublattice, and rj is a vector 
from a spin to one of its z neighbors. We now need to 



where v is the sound wave speed, assumed isotropic, and 
Kitom is the volume per atom. This expression is similar 
to Eq. (74) except that the prefactor is three instead of 
two, due to the three sound wave polarization directions 
for each sound wave mode (two mutually perpendicular 
transverse polarizations and one longitudinal polariza- 
tion) which are assumed to have the same wave speed v 
in the Debye model. 
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2. Application of the Spin Wave Theory for the Magnetic 
Heat Capacity to the Ji-J2-Jc Heisenherg Model and 
BaMn2As2 

From the expressions for the spin wave velocities in the 
Ji-J2-Jc model in Eq. (9), one has 



VxVyV;, = V^f,Vc = ^3 i ^ + 2Ji 



3/2 




For BaMn2As2, there are two formula units, or four Mn 
atoms, per unit cell with volume a^c. The volume per 
spin is thus 



^Din — 



(77) 



Dividing Eq. (76) by Eq. (77) gives 



VxVyVz ^ 16\/2(Ji5)3 / J, y^/^ / 2J2\ fZ 




(79) 




Inserting Eq. (78) into (74) gives 

= ""'^ 

60\/2(Ji5/fcB)3 



Ftom Table II, the exchange constants from the neu- 
tron data arc Ji^/fce = 380 K, J2/J1 = 0.29 and 
Jc/Ji = 0.09. Inserting these values into Eq. (79) gives 
the calculated value 

^gw = 0.13 mJ/mol spins for BaMn2As2. (80) 

From Eq. (63), the observed /3 value per mole of Mn 
spins is 0.325 mJ/mol spins K"'. The calculated /3sw 
value is thus 40% of the measured value, so the observed 
/3 value contains a significant magnetic contribution if the 
anisotropy gap in the spin wave spectrum is negligible. 
However, an anisotropy gap would reduce the spin wave 
contribution to the heat capacity to exponentially small 
values at low temperatures. 



IX. MONTE CARLO SIMULATIONS OF THE 
MAGNETIC SUSCEPTIBILITY AND MAGNETIC 
HEAT CAPACITY IN THE Ji-Ji-Jc MODEL 

Both our classical and quantum Monte Carlo simula- 
tions were carried out within the framework of the Ji- 
J^-Jc Heisenberg model introduced above in Sec. III. We 
have calculated the magnetic heat capacity and magnetic 
spin susceptibility versus temperature for various size lat- 
tices of quantum spins S = 1/2, 1, 3/2, 2, and 5/2, and 



for the classical model. We first motivate the scaling of 
the axes of our theoretical plots of x(T), remark on the 
temperature regime over which this scaling is expected 
to hold, and then present our Monte Carlo simulation 
results. Then we will compare our predictions for the 
magnetic susceptibility with the experimental suscepti- 
bility data for BaMn2As2 above Tn in Fig. 7 to obtain 
additional estimates of the exchange constants in this 
compound. 



A. Scaling of the Theoretical x(T) Axes 

Using Eqs. (23) and (25) in the Heisenberg "J model" 

for a bipartite spin lattice with equal NN exchange, the 
Curie- Weiss law (22) can be rewritten as 



xJ 



1 



JS(S+1) 



(81) 



The quantity on the left-hand side of Eq. (81) is the the- 
orist's definition of "x", which is the susceptibility per 
spin, in units of 1/J, with ^/zb set equal to 1. On the 
right-hand side, we sec that if we use a temperature scale 
defined by k^T /[JS{S+1)\, then all spin lattices with the 
same coordination number z but with different J and/or 
S will all follow the same universal curve at high temper- 
atures. Therefore in this paper we scale the calculated 
susceptibilities when J2, Jc = as 



versus 



JiS{S + l) 



(82) 



This is the same scaling of the temperature axis as for 
the magnetic heat capacity in Eq. (AlO). 

In the Ji-J2-Jc model, according to Fig. 2 there are 
four in-plane next-nearest-neighbor interactions 



J2 = aJi, 
two NN interactions along the c-axis 

Jc = iJi, 



(83) 



(84) 



in addition to the z\ = A = z nearest-neighbor interac- 
tions Ji. When these additional interactions are present, 
according to Eq. (26) the Weiss temperature becomes 



^Ji(l + a + 7/2)g(g+l) 

3fcB 



(85) 



and the form of the new Curie Weiss law corresponding 
to Eq. (81) is 



xJi(l + a + 7/2) 



3fcBr , ■ 

Ji(l+OL+-i/2)S{S+l) ^ 



(86) 



A more accurate high-temperature scaling is obtained in 
this case by replacing Ji in Eq. (82) by Ji -|- J2 + </c/2 = 



23 



Ji(l + a + 7/2) and scaling the data according to 



Semiclassical Square Lattice 



xJi(l + a + 7/2) 



versus 



Ji(l + a + 7/2)S'(5+l)' 
(87) 

The scahngs in Eqs. (82) and (87) are expected to be 
universal with respect to the spin and the exchange con- 
stants only at "high" temperatures. Appendix A shows 
that the calculations begin to deviate from the Curie- 
Weiss behavior when 1 /T^ and higher order terms in the 
two-spin correlation functions become significant com- 
pared to the 1/T terms with decreasing T. 



B. Classical Monte Carlo Simulations 

The classical Monte Carlo (CMC) simulations were 
performed on periodic L x L clusters for Jc = and 
ioi L X L X Lc clusters for Jc 7^ using a hybrid al- 
gorithm that combines Metropolis and over-relaxation 
sweeps.'*^ In order to obtain statistically reliable data we 
have generated ^ 10'^ configurations at each tempera- 
ture and then averaged the results over 50 independent 
annealing runs. 

The spin Hamiltonian for our classical Monte Carlo 
simulations is the classical analogue of the quantum spin 
Hamiltonian (2), given by 

classical — 

+J,S^ ^^ ■ Sfc + gf^BH SI 

(il) i 



where S is the magnitude of the spin, S is a classical spin 
unit vector, and Si - Sk — cosOij. According to Eq. (88), 
the exchange parameters Jq are always combined with 
the classical spin magnitude S in the combination J^S^ . 

In the following, we first consider our simulations for 
Jc — Q and then for Jc 7^ 0. 



1. Jc = Q 
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FIG. 17: (Color online) (a) Normalized magnetic spin 
susceptibility xJ^I 9^ A) versus normalized temperature 
k^T /[J\S(S ~\- 1)] for the classical spin S Heisenberg square 
lattice with Jc = and J2/J1 = -0.4 to 0.4, with S{S + 1) 
replacing 5*^. The lattice size is 80 x 80 in each case, (b) 
Expanded plot at low temperatures (open circles) of the data 
for J2 = in (a). The error bars are shown and are inside 
the open circles. A linear fit to the data up to a reduced tem- 
perature of 0.14 is shown as the solid straight line, and the 
dashed line is an extrapolation of the fit. The coefficients of 
the fit are listed in Eq. (89). 



The semiclassical magnetic spin susceptibilities x ver- 
sus T for the square lattice calculated using CMC sim- 
ulations on 80 X 80 spin lattices are shown in Fig. 17(a) 
for Jc = and J2/J1 — —0.4 to 0.4. Here, the term 
"semiclassical" means that S*^ in the final result of the 
classical simulations is replaced by the quantum mechan- 
ical expectation value (5^) ~ 5(5 + 1). This replacement 
allows the classical simulations to merge smoothly with 
the quantum Monte Carlo simulations (see Fig. 23 be- 
low) . We carried out simulations of various other L x L 
lattice sizes with L — 10-100 for J2/J1 — and 0.2 and 
found that finite-size corrections to both the calculated 
magnetic susceptibility and magnetic heat capacity are 
negligible for L > 50. 



The x{T) data in Fig. 17(a) show two interesting 
trends. First, at high temperatures the Curie- Weiss law 
C/{T + 6) is obtained, in which the (positive) Weiss tem- 
perature 9 is proportional to the sum of all interactions 
of a given spin with its neighbors according to Eq. (26). 
Thus for a negative (ferromagnetic) J2 that partially can- 
cels the positive Ji , the susceptibility increases at a fixed 
T, and for a positive J2 it decreases. Second, at low 
temperatures this trend is reversed. A negative ferro- 
magnetic J2 is nonfrustrating with respect to Ji, and 
reinforces the short-range ordering that causes the peak 
in xiT). This moves the peak up in temperature and 
suppresses the susceptibility in the short-range ordered 



24 



0.00 



Curie-Weiss Law 




J /J = -0.4, -0.2, -0.1, 

2 1 , , , 

0.0, 0.1, 0.2, 0.3, 0.4 



l<J/[(J^+J^)S(S + 1)] 



FIG. 18: (Color online) Normalized magnetic spin suscep- 
tibility x{Ji + >^i)/(-^5^Mb) versus normalized temperature 
ksT/KJi + J2)S{S + 1)] determined from classical Monte 
Carlo simulations for the spin S Heisenberg square lattice 
with Jc — and J2/J1 = —0.4 to 0.4, as indicated, with 
S{S + 1) replacing the classical S^. The lattice size is 80 x 80 
in each case. The axis scaling is superior at high tempera- 
tures to that in Fig. 17(a). A plot of the Curie- Weiss law in 
Eq. (86) is shown by the black dashed curve. 



state at low temperatures below the peak temperature. 
The opposite behavior is found for a positive AF J2 which 
is frustrating with respect to Ji. This J2 suppresses the 
short-range AF ordering, which decreases the peak tem- 
perature and increases the susceptibility below the peak 
temperature compared to the case when J2 = 0. 

These trends are illustrated in a different way if the 
best high-temperature scaling for these plots, given in 
Eq. (87), is used, as shown in Fig. 18. In addition, the 
Curie- Weiss law from Eq. (86) is plotted in Fig. 18 as the 
blue dashed line. From a comparison of the simulation 
data with the Curie- Weiss prediction, one sees that the 
two-spin correlations higher order than present in the 
Curie- Weiss regime ^/T) begin to become observable 
on the scale of the figure for T < 5{Ji/kB)S{S + 1). 
According to Eq. (85), this latter value is about four 
times the Weiss temperature 9, which has the value 4/3 
on the horizontal scale in Fig. 18. 

The data in Fig. 17(a) for J2 = were obtained down 
to a reduced temperature of 0.01 as shown in the ex- 
panded plot in Fig. 17(b). The lowest temperature data 
are linear in T. A linear fit yielded 



xJi 



0.08333(2) + 0.0218(3) 



Ji5(5 + 1) 



(89) 



as shown by the solid line in Fig. 17(b). According to 
Takahashi's modified spin wave theory for the AF square 
lattice, the classical limit (A9) in Ref. 42 reads 
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FIG. 19: (Color online) Classical Monte Carlo (CMC) simu- 
lations on 80 X 80 spin lattices of the magnetic heat capacity 
Cmag divided by the molar gas constant R versus the scaled 



temperature ksT / J-i_S{S + 1) for Jc = and J2/J1 
to 0.4 as shown. 



-0.4 



The zero-temperature reduced susceptibility 1/12 « 
0.08333 in Eq. (90) is the same as our value in Eq. (89) 
to within the errors of our Monte Carlo data, but the 
theoretical initial slope l/(247r) w 0.01326 is too small 
compared to our Monte Carlo value in Eq. (89). On the 
other hand, in a 1/D expansion where D is the dimen- 
sionality of the spins [D = 3 here), for the classical square 
spin lattice at low T Hinzke et al.''^ obtained 



1 

12 



1 kBT 
32 JiS'2' 



(91) 



The zero temperature susceptibility is the same as our 
and Takahashi's value but Hinzke et al.'s initial slope is 
1/32 = 0.03125, which this time is larger than our Monte 
Carlo value in Eq. (89). Thus our value of the initial slope 
is bracketed by the predictions of the modified spin wave 
theory and the l/D expansion. 

The magnetic heat capacity is plotted in Fig. 19 ac- 
cording to Eq. (AlO) versus the scaled temperature for 
exchange constant ratios J2/J1 — —0.4 to 0.4 on 80 x 80 
spin lattices. The broad peaks in the curves decrease in 
temperature with increasing J2. This is understandable 
in terms of the enhancement of short-range AF order for 
ferromagnetic (negative) J2 , which increases the temper- 
ature of the peak, and the frustration effect for antiferro- 
magnetic (positive) J2, which decreases the temperature 
of the peak. The peak values and the temperatures at 
which they occur are listed in Table V. It is interest- 
ing that the variation in Cniag(r) with J2 depends on 
the sign of J2 , in contrast with expectation from the first 
term in the HTSE in Eq. (A9) in which the uniform J ap- 
pears as the square and is hence independent of the sign. 
Thus one cannot replace zJ"^ in Eq. (A9) by Jf^. This 
constraint is not present when calculating the Weiss tem- 
perature in the Curie- Weiss law from Eq. (26), in which 
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TABLE V: Parameters obtained from semiclassical (SC) 
Monte Carlo and quantum Monte Carlo simulations of the 
magnetic heat capacity Cmag(r) for the square lattice with 
no interlayer coupling (Jc = 0). Here, C™ag is the value of 
Cmag at a magnetic ordering peak, T"^^^ is the temperature at 
which the maximum occurs, and SC means we have replaced 
by S{S + 1) in the temperature scaling of the classical 
Monte Carlo data. Note that the value of T"""'' is different 
than the temperature of the maximum in the magnetic sus- 
ceptibility. 



5" 


lattice size 


J2/J1 


R 


JiS{S+l) 


1/2 


32 X 32 





0.4606(7) 


0.801(2) 


1 


32 X 32 





0.885(2) 


0.690(4) 




64 X 64 





0.879(2) 


0.700(3) 


3/2 


32 X 32 





1.159(2) 


0.674(1) 


2 


32 X 32 





1.325(2) 


0.673(1) 




64 X 64 





1.295(2) 


0.684(2) 


5/2 


32 X 32 





1.428(2) 


0.673(2) 


SC 


80 X 80 


-0.4 


1.801(5) 


1.055(3) 


SC 


80 X 80 


-0.2 


1.752(3) 


0.861(3) 


SC 


80 X 80 


-0.1 


1.699(1) 


0.777(2) 


SC 


80 X 80 





1.666(1) 


0.678(1) 


SC 


80 X 80 


0.1 


1.621(2) 


0.575(1) 


SC 


80 X 80 


0.2 


1.567(2) 


0.471(2) 


SC 


80 X 80 


0.3 


1.498(3) 


0.352(3) 


SC 


80 X 80 


0.4 


1.378(3) 


0.232(2) 



one includes the interactions of a given spin with all of 
its neighbors algebraically and on the same footing. 



2. Jc / 

The classical Monte Carlo simulations do not produce 
the same results as the molecular field theory does be- 
cause the interaction between a spin and its neighbors 
is not approximated by the interaction of the spin with 
the average spin of its neighbors as in the molecular field 
theory. In particular, according to the Mermin- Wagner 
theorem,'** a Heisenberg spin system in one or two di- 
mensions, as in the J1-J2 model with only intraplanar 
exchanges, should not show long-range magnetic ordering 
at finite temperature. This theorem is respected in our 
classical simulations, but not in molecular field theory. 
On the other hand, when the simulations are carried out 
with Jc 7^ 0, we find that long-range AF ordering does 
occur, as expected. Because a uniform magnetic field 
does not directly couple to the AF order parameter [the 
staggered moment, see Eq. (B16)], these AF phase transi- 
tions have rather subtle effects on the calculated uniform 
susceptibility. They are much more clearly manifested in 
the magnetic heat capacity which we will also present, 
and would also be clearly delineated in calculations of 
the staggered susceptibility in which the applied mag- 



2.50 




kJ/[J^S(S + 1)] 




kJ/[J^S(S + 1)] 

FIG. 20: (Color online) (a) Magnetic heat capacity and (b) 
magnetic susceptibility x versus temperature T from classical 
Monte Carlo simulations for the spin S Heisenberg square 
lattice with J2/J1 = 0.1 and with Jc/Ji = 0, 0.02, 0.05 and 
0.1, as indicated. The order of the curves from top to bottom 
is the same as in the respective figure legends. The lattice 
size in each case is 20 x 20 x 10. 



netic field has opposite directions for the two sublattices. 

Throughout this section, we replace the classical vari- 
able by its quantum-mechanical counterpart S{S + 
1). We show our results in two formats. First, in 
Fig. 20 are shown the magnetic heat capacity Cmag and 
the spherically-averaged magnetic susceptibility x versus 
temperature T for fixed J2/ Ji — 0.1 and variable Jc = 0, 
0.02, 0.05 and 0.1, where Ji, J2 and Jc are all antiferro- 
magnetic. From 20(a), one sees that Cmag{T) for Jc — 
just shows a broad peak characteristic of short-range AF 
order. However, the Cmag quickly and clearly shows a 
cusp-like behavior with increasing Jc at temperatures Tn 
corresponding to long-range AF order. Second, in Fig. 21 
are shown Cmag(T) and xi^) for fixed Jc/Ji — 0.1 and 
variable J2/J1 — 0-0.4, where Ji, J2 and Jc are again 
all antiferromagnetic. For all combinations of J2/J1 and 
Jc/ Ji we have studied, to within the errors the peak in 
Cmag at Tn coincides in temperature with the peak in 
d{xT)/dT on the low-T side of the broad peak in x{T), 
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FIG. 21: (Color online) (a) Magnetic heat capacity Cmag and 
(b) spherically averaged magnetic susceptibility x versus tem- 
perature T from classical Monte Carlo simulations for the 
spin S Heisenberg square lattice with J2/J1 = to 0.4, as in- 
dicated, and with fixed Jc/Ji = 0.1. The order of the curves 
from right to left is given in the figure (a) legend. The lattice 
size in each case is 20 x 20 x 10. 



in agreement with the Fisher relation. '^^ We note that the 
Tn and shape/magnitude of Cmag at Tn determined in 
our simulations may be affected by finite size effects. 

The temperatures Tn of the peaks in Cmag {T) and the 
values of Cmag(?N) at the peak, versus the exchange con- 
stant ratios J2/J1 and Jd Ji, are listed in Table VI and 
plotted in Figs. 22(a) and (b). For J2 = 0, our Tn values 
are lower by ^ 1% than the values obtained by Yasuda et 
al. for Jc/Ji = 0.02, 0.05 and O.l.^^ For J2 = 0, a good 
fit to Tn versus for various spin values was obtained 
in Ref. 46 using the expression'*^ 



/crTn 



A 



JiS[S+l) B~\n{Jc/Jiy 



(92) 



where different values of the constants A and B were 
required for different spin values. We fitted the classical 
Monte Carlo data for J2 = in Fig. 22(a) using Eq. (92) 



TABLE VI: Temperatures Tn of the magnetic ordering peaks 
in the molar magnetic heat capacity Cmag and the values of 
Cmag(rN) at the peak, versus the exchange constant ratios 
J2/J\ and Jc/Ji, obtained from classical Monte Carlo simu- 
lations. The Tn and Cmag values may be affected by finite 
size effects. Here R is the molar gas constant. 



J2/J1 



Jc/Jl 



k^Tn/[JlS{S + 1)] Cmag(rN)/i? 






0.02 


0.717(2) 


1.890(3) 




0.05 


0.773(1) 


2.139(5) 




0.10 


0.842(1) 


2.29(2) 


0.1 


0.02 


0.611(1) 


1.879(4) 




0.05 


0.667(1) 


2.151(6) 




0.10 


0.735(1) 


2.318(5) 


0.12 


0.06 


0.661(1) 


2.190(6) 


0.2 


0.02 


0.502(1) 


1.880(8) 




0.05 


0.557(1) 


2.14(1) 




0.10 


0.619(1) 


2.32(1) 


0.3 


0.02 


0.390(1) 


1.839(2) 




0.05 


0.439(1) 


2.125(9) 




0.10 


0.497(1) 


2.343(7) 


0.4 


0.02 


0.263(1) 


1.784(8) 




0.05 


0.313(1) 


2.06(1) 




0.10 


0.360(1) 


2.26(1) 



and obtained a good fit with the values 

A = 7.70 
B = 6.87, 



(93) 



which are both about a factor of two larger than the 
respective values A — 3.96 and B = 3.01 obtained in 
Ref. 46 for the classical limit S = co. The fit is shown as 
the solid red curve for J2/J1 = in Fig. 22(a). 

From Fig. 22(a), a positive antiferromagnetic J2 frus- 
trates the G-type AF ordering and depresses Tn approx- 
imately linearly with J2, and the Tn for each Jc value 
extrapolates to zero at J2/J1 ~ 0.6, which is close to the 
value of 0.5 from Eq. (4) at which one classically expects 
the G-type AF order to become unstable with respect to 
the stripe AF order. Therefore, we fitted the dependence 
of Tn on J2/J1 of ah the data for JjJx = 0.02, 0.05 
and 0.1 together using the expression 



Tn(Jc, J2) 



1 



J2/J1 
Ri 



R2 



(94) 



Tn(Jc, J2 = 0) 
(not shown) and obtained the values 

Ri = 0.644 (95) 
i?2 = 1.082. 

The fits are shown as the solid curves in Fig. 22(b). Thus 
the global function to fit all of our Tn( Jc/ Ji, J2/J1) data 
is 

R2' 



/cbTn 



A 



JiS'(S'+l) B -\n{Jc/Ji) 



1 



J2/J1 
Ri 



(96) 
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FIG. 22: (Color online) Antiferromagnetic ordering tempera- 
ture Tn versus (a) the interlayer coupling Jc for fixed values 
of J2/J1 and (b) J2/J1 at fixed Jc/ Ji- In (b), the order of the 
curves from top to bottom is the same as in the figure legend. 
The solid curves are a global fit to all the data by Eq. (96), 
using the parameters in Eqs. (93) and (95). 



The fits for Tn versus Jc at fixed J2/J1 = 0.1-0.4 
are shown as the solid curves in Fig. 22(a). We see 
that Eq. (96), together with the four parameters in 
Eqs. (93) and (95), provide a good global fit to all fif- 
teen Tn(Jc/Ji, J2/J1) data points in Fig. 22(a) from 
our CMC simulations. 



C. Quantum Monte Carlo Simulations 

1. Magnetic Susceptibility 

Our quantum Monte Carlo (QMC) simulations 
were carried out with the ALPS'*^''** directed loop 
application"*^ in the stochastic series expansion 
framework^° using version ALPS 1.3. Up to about 
1 X 10^ sweeps were carried out for the 32 x 32 lat- 
tice and the sign- free situation J2 — 0. In order to 
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FIG. 23: (Color online) Normalized magnetic spin sus- 
ceptibility xJ'i-K^g^ A) versus normalized temperature 
ksT /[J-iS{S + 1)] determined from quantum Monte Carlo 
(QMC) simulations for quantum Heisenberg square lattices 
with (a) J2/J1 = and spins S = 5/2, 2, 3/2, 1, and 1/2, 
and (b) J2/J1 = 0.1 and spins S — 5/2 and 2. The data for 
the semiclassical model (green) with J2/J1 =0 and 0.1 are 
included in the two panels, respectively. The lattice size for 
each simulation is indicated in the respective figure. 



compensate for the sign problem introduced by J2 > 
this was increased to about 2 x 10^* sweeps on the 
6x6 lattice for J2 — 0.1. QMC simulations have been 
previously reported for 5* = 1 over the temperature 
range 0.5 < kBT/JiS{S + 1) < 5 by Harada et al.^^ 
We have extended these simulations to much lower 
temperatures (see Fig. 24). 

Our QMC simulations of the magnetic spin suscep- 
tibilities versus temperature for the square lattice with 
quantum spins 1/2 to 5/2 with J2 — are shown in 
Fig. 23(a), and for J2/J1 = 0.1 and 5 = 2 and 5/2 in 
Fig. 23(b). Various parameters obtained from these and 
the above semiclassical data are listed in Table VII as 
described in the caption. We checked by comparison of 
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FIG. 24: (Color online) Normalized QMC magnetic spin 
susceptibility X>^i/(-'^5^Mb) versus normalized temperature 
ksT /[JiS{S + 1)] data (symbols) at low temperatures from 
Fig. 23(a) for spins S = 1/2 (bottom), 1, 3/2, 2, and 5/2 (top). 
The error bars for the QMC data are also plotted. The corre- 
sponding Hasenfratz-Niedermayer+spin wave theory predic- 
tions (HN + SWT) for the low-temperature behaviors^^'^^ 
are also shown for these spin values as solid curves. 



et al. as^^ 
X^(0)Ji 



1 0.034446 959 42 

8 



S 

0.002 040 06(7) 

^2 



+ 0(3-^) 



(99) 



and 



Ji 



-0.117 628 25445- 0.010207 9873 
0.00316(2) 



S 



+ o{s-'). 



(100) 



The low-temperature QMC data from Fig. 23(a) are 
shown in Fig. 24 together with the above predictions of 
Hasenfratz and Niedermayer (HN) combined with the 
SWT results of Hamer et al. The lowest-temperature 
QMC data for spins 1/2 to 5/2 in Fig. 24 are all seen to 
be in good agreement with the HN + SWT predictions. 
High resolution calculations of x{T) from the literature 
for the S' = 1/2 square lattice Heisenberg antiferromag- 
net also confirm the form of Eq. (97).^'*'^^ For S ~ I 
our value of x(0) from Fig. 24 disagrees with the value 
0.07197 given in Ref. 46. It was claimed in Ref. 42 that on 
the basis of spin wave theory, the next-order term above 
the term in Eq. (97) is 0{T^), as in Eq. (90) above, 
in disagreement with Eq. (97). However, the next higher 
order term is indeed the term.^^ 



the data for the 32 x 32 lattice with 64 x 64 lattice data 
for 5 = 1 and 5 = 2 (see Table VII) that the 32 x 32 
lattice data in Figs. 23 and 24 (below) are close to the 
thermodynamic limit; i.e., finite size effects are smaller 
than the size of the symbols (except probably the lowest 
temperature datum for 5 = 1/2 in Fig. 24). 

Hasenfratz and Niedermayer obtained the low temper- 
ature limit of the spin susceptibility of the Heisenberg 
antiferromagnet on a square lattice from chiral perturba- 
tion theory, given by^'^ 



X(T) 



2X±(0) 



2Trps 



(ML 

\27rps 



(97) 

where ps is the spin wave stiffness, X-l(O) is the zero- 
temperature perpendicular susceptibility given by 



X_l(0)Ji psJia'^ 



Ng'pl {hey 



(98) 



c is the spin wave velocity, and a is the square lattice 
parameter. The X-l(O) ^-nd ps depend on the spin S and 
were calculated using spin wave theory (SWT) by Hamer 



2. Magnetic Heat Capacity 

The magnetic heat capacity Cmag versus temperature 
data from our QMC simulations for the square lattice 
with only NN interactions (J2 = Jc = 0) are shown 
in Fig. 25 for spins S = 1/2, 1, 3/2, 2 and 5/2. Also 
shown for comparison are our CMC heat capacity data 
for J2 = Jc = from Fig. 20(a) and the first term in 
the HTSE Cmag oc for the magnetic heat capacity 

from Eq. (A9) using the nearest-neighbor coordination 
number z = 4. The CMC and QMC data exhibit this 
HTSE behavior for temperatures T > 2Ji5(5 + l)/fcB. 
The values of the heat capacities of the peaks in the sim- 
ulation data for the spin values 5* = 1 to 5 = 5/2 and the 
temperatures at which they occur are listed in Table V 
above. 

We checked finite-size effects associated with the QMC 
data by simulating Cmag for 64 x 64 5 = 1 and S* = 2 
lattices for comparison with the 32 x 32 lattices in Fig. 25. 
On the scale of the figure, the 64 x 64 data (not shown) 
were close to the 32 x 32 lattice size data. For example, 
the peak heights differ by less than 2 percent between the 
simulations for the different size lattices (see Table V). 

According to Eq. (41), the Neel temperature in MET 
occurs in Fig. 25 at a value of 4/3 on the horizontal scale 
and with a heat capacity jump on cooling below Tn given 
by Eq. (57) as ACmag/i? = 3/2 for S = 1/2 and = 5/2 
for S — 00. The data in Fig. 25 are very different from 
these MET predictions due to the presence of short-range 
magnetic ordering and the lack of long-range magnetic 
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TABLE VII: Calculated values of (i) the maximum spin susceptibilities x™"^ ^md (ii) the temperatures T™*'' at which they 
occur for different lattice sizes, J2/J1 values and quantum and semiclassical S values, (iii) the product T^ma'=j"nax^g^ where C 
is the Curie constant per mole of spins in Eq. (23), (iv) the value of the product •^ma=<:j'ma^ for the listed S assuming g — 2, and 
(v) the value of the product for an alternate semiclassical value of 5 = 5/2 assuming g — 2. By "semiclassical" (SC) is meant 
that in the final result of the classical calculation is replaced by the quantum mechanical expectation value (S^) = S{S + 1). 



In the last column is 


listed the molecular field Neel temperature Tn, normalized by JiS{S + l)/kB, according to Eq. (41). 


S 


lattice size 


T / T 
J2/J1 


T / T 

Jc/Jl 


^max 


j,^ ymax 


^max /jnmax 


m a,5f /Tim ax 

(cm^K/mol) 


m ax mm pt.'x 
^maxymax 

(cm^K/mol) 






JlS{S+l) 


3C 


JlS(S+l) 


1/2 


32 X 32 








0.09370(3) 


1.248(3) 


0.1169(3) 


0.1316(4) 




1.33 


1 


32 X 32 








0.10438(5) 


1.090(3) 


0.1138(3) 


0.3415(11) 




1.33 




64 X 64 








0.10424(3) 


1.085(3) 


0.1131(4) 


0.3394(11) 






3/2 


32 X 32 








0.10790(3) 


1.050(4) 


0.1133(5) 


0.638(3) 




1.33 


2 


32 X 32 








0.10952(4) 


1.030(4) 


0.1128(5) 


1.016(4) 




1.33 




64 X 64 








0.10957(6) 


1.038(5) 


0.1137(6) 


1.024(6) 






5/2 


32 X 32 








0.11040(3) 


1.018(3) 


0.1144(8) 


1.50(1) 




1.33 


2 


6x6 


0.1 





0.10882(6) 


0.966(5) 


0.1051(6) 


0.946(6) 




1.20 


5/2 


6x6 


0.1 





0.10966(6) 


0.959(4) 


0.1052(5) 


1.381(6) 




1.20 
















S = 2 


5 = 5/2 




SC 


80 X 80 


-0.4 





0.11115(5) 


1.453(4) 


0.1615(5) 


1.454(5) 


2.121(8) 


1.87 


SC 


80 X 80 


-0.2 





0.11173(5) 


1.225(4) 


0.1369(5) 


1.232(5) 


1.797(8) 


1.60 


SC 


80 X 80 


-0.1 





0.11210(4) 


1.120(3) 


0.1256(4) 


1.130(4) 


1.649(5) 


1.47 


SC 


80 X 80 








0.11235(3) 


0.999(2) 


0.1122(3) 


1.011(3) 


1.470(4) 


1.33 


SC 


80 X 80 


0.1 





0.11274(3) 


0.876(3) 


0.0988(4) 


0.889(3) 


1.297(5) 


1.20 


SC 


20 X 20 X 10 


0.1 


0.02 


0.11216(6) 


0.874(5) 


0.0980(6) 


0.883(6) 


1.287(8) 


1.21 


SC 


20 X 20 X 10 


0.1 


0.05 


0.1113(1) 


0.878(7) 


0.0977(9) 


0.880(8) 


1.28(1) 


1.23 


SC 


20 X 20 X 10 


0.1 


0.1 


0.1098(1) 


0.896(5) 


0.0984(6) 


0.886(6) 


1.29(1) 


1.27 


SC 


20 X 20 X 10 


0.12 


0.06 


0.11126(8) 


0.853(7) 


0.0949(9) 


0.854(8) 


1.25(1) 


1.21 


SC 


80 X 80 


0.2 





0.11325(2) 


0.750(2) 


0.0849(3) 


0.765(2) 


1.115(3) 


1.07 


SC 


20 X 20 X 10 


0.2 


0.02 


0.11275(5) 


0.749(10) 


0.0844(12) 


0.76(1) 


1.11(1) 


1.08 


SC 


20 X 20 X 10 


0.2 


0.05 


0.11187(10) 


0.755(10) 


0.0845(12) 


0.76(1) 


1.11(1) 


1.10 


SC 


20 X 20 X 10 


0.2 


0.1 


0.11029(5) 


0.767(7) 


0.0846(8) 


0.762(7) 


1.11(1) 


1.13 


SC 


80 X 80 


0.3 





0.11391(3) 


0.616(2) 


0.0702(3) 


0.632(3) 


0.921(3) 


0.93 


SC 


20 X 20 X 10 


0.3 


0.02 


0.1134(10) 


0.614(10) 


0.070(1) 


0.63(2) 


0.91(3) 


0.95 


SC 


20 X 20 X 10 


0.3 


0.05 


0.11257(7) 


0.627(7) 


0.0706(8) 


0.635(8) 


0.93(1) 


0.97 


SC 


20 X 20 X 10 


0.3 


0.1 


0.1108(10) 


0.652(7) 


0.072(2) 


0.65(1) 


0.95(2) 


1.00 


SC 


80 X 80 


0.4 





0.11509(4) 


0.468(2) 


0.0539(3) 


0.485(2) 


0.707(3) 


0.80 


SC 


20 X 20 X 10 


0.4 


0.02 


0.1147(10) 


0.468(7) 


0.054(2) 


0.48(1) 


0.70(2) 


0.81 


SC 


20 X 20 X 10 


0.4 


0.05 


0.1136(10) 


0.467(10) 


0.053(2) 


0.48(1) 


0.70(2) 


0.83 


SC 


20 X 20 X 10 


0.4 


0.1 


0.11175(5) 


0.497(3) 


0.0555(4) 


0.500(3) 


0.729(5) 


0.87 



ordering'*^ in these two-dimensional spin lattices at finite unit cell. The spin wave velocity v for the AF Heisenberg 
temperatures. square lattice is^^'^'^ 

The expression of Hasenfratz and Niedermayer for the 
low temperature magnetic heat capacity of the Heisen- 
berg antifcrromagnct on a square lattice from chiral per- 
turbation theory, per mole of spins, is^^ 



6C(3)i? „ , ^,^4^ 

— = 2V2SJ1 



hv 

where C(.t) is the Riemann zeta function with ^(3) « a 



2S {2S) 

1.20206 and a is the length of an edge of the square lattice (102) 



0.157947421 0.02152 
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From Eqs. (101) and (102) one obtains 
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FIG. 25: (Color online) Magnetic heat capacity Cmag divided 
by the molar gas constant R versus normalized temperature 
k-BT/[JiS{S + 1)] for the square lattice with only nearest- 
neighbor couplings Ji (exchange constants J2 = Jc = 0). 
Quantum Monte Carlo (QMC) data (filled circles) for spins 
S = 1/2 (bottom), 1, 3/2, 2, and 5/2 (top) are shown, to- 
gether with classical Monte Carlo (CMC) data from Fig. 20(a) 
(solid black curve at the top) and the first term in the high- 
temperature series expansion HTSE for the magnetic heat ca- 
pacity from Eq. (A9) using z = 4 (dashed black curve). The 
error bars for the QMC data are also plotted, and the lines 
connecting the data points are guides to the eye. The order 
of the Monte Carlo data from top to bottom is the same as 
in the figure legend. 
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The factor A{S) is 0.2063, 0.2441, 0.2578, 0.2649 and 
0.2692 for S* = 1/2, 1, . . . , 5/2, respectively. Figure 26 
shows our low-temperature Cnia.g{T) data for S = 1/2, 
1 and 3/2 plotted according to Eq. (103). The approxi- 
mate extrapolated zero-temperature values are in accord 
with the above A(S) values to within the data error 
bars. After setting A{S) to the above respective fixed 
values, the initial slopes were estimated by fitting the 
data by Eq. (103), yielding D{S) = 1.8, 0.31 and 0.14 
for S — 1/2, 1 and 3/2, respectively, as shown by the 
respective dotted lines in Fig. 26. The slope D decreases 
significantly with increasing S but remains positive from 
S = 1/2 up to 5 3/2. The sign of D{S) was indeed 
predicted by Hofmann to be positive using the effective 
Lagrangian method.^* 



X. COMPARISON OF MONTE CARLO 
SIMULATIONS OF THE MAGNETIC 
PROPERTIES WITH EXPERIMENT 
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FIG. 26: (Color online) Low-temperature magnetic heat ca- 
pacity Cmag times 5*^ divided by the molar gas constant R 
and versus normalized temperature squared, {ksT / Jif , 
for the square lattice with only nearest-neighbor couplings Ji 
(exchange constants J2 = Jc = Q)- Quantum Monte Carlo 
(QMC) data (filled symbols) for spins S = 1/2, 1, and 3/2 
are shown, together with error bars. The dotted lines are 
fits of the data by Eq. (103) assuming the specific A{S) y- 
intercepts given in the text below Eq. (103). Only the slopes 
D{S) were fitted. The fitting ranges of were 0-0.06 for 
S = 1/2, 0-0.25 for 5 = 1 and 0-0.5 for 5 = 3/2. 



A. Neel Temperature 

Using Eq. (96), we can predict the Neel temperature 
from the values JiS/k^ = 380 K, J2/J1 = 0.29, and 
Jcl J\ = 0.09 in Table II obtained from the fit of the in- 
elastic neutron scattering data by spin wave theory. Note 
that the neutron fit only provides products of the J val- 
ues with S. Using the above-given parameters, Eq. (96) 
predicts 



1. Tn 



550 K for S 



2. Tn = 640 K for S* = 5/2. 

A comparison of these values with the experimental value 
Tn — 620-625 K clearly favors spin 5/2 over spin 2 for 
the Mn ions. Indeed, the Tn calculated for S ~ 5/2 is 
remarkably close to the observed value. 



B. Magnetic Susceptibility 

Our tables of calculated susceptibities are in the form 
of Eq. (82). A very useful quantity for comparison 
with experimental susceptibility data is the product of 
the scaled maximum susceptibility X™'*'' Ji/A^g^/il and 
the scaled temperature at which the maximum occurs 
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FIG. 27: (Color online) Product of the maximum suscepti- 
bility Tj^™''" per mole of spins and the temperature y™^^ at 
which it occurs for the J1-J2 model on a Heisenberg square 
spin lattice for spins (a) 5 = 2 and (b) S — 5/2, versus the 
ratio J2 / Ji of the diagonal next-nearest-neighbor coupling to 
the nearest neighbor coupling. Here a g-factor g = 2 was as- 
sumed. The data were obtained using quantum Monte Carlo 
(blue, QMC) and classical Monte Carlo (red, CMC) simula- 
tions on the lattice sizes indicated. The i^maxymax yg^j^gg 
tained from CMC for the three-dimensional lattices are very 
insensitive to the coupling Jc between layers up to at least 
Jc/Ji ~ 0-1 as shown in Table VII. The solid black curves in 
(a) and (b) are least square fits to the respective data by the 
second order polynomials in Eqs. (105). 



A:bT™^''/ JiS'(S'+l). Setting N equal to Avogadro's num- 
ber Na so that X is the susceptibility per mole of spins, 
the product of these two variables is 



JiSiS + 1) 



X 



maxT^max 



3C 



(104) 



where C is the Curie constant in Eq. (23). This prod- 
uct does not contain any exchange constants and hence 
is a potential diagnostic for the value of the spin S from 



TABLE VIII: Parameters determined from a fit of magnetic 
susceptibility data for BaMn2As2 by classical Monte Carlo 
simulations. The value Jc/Ji = 0.1 was chosen to agree 
with the inelastic neutron scattering and band theoretical val- 
ues which are both close to 0.1. The values of J2/J1 were 
fixed at the listed values by comparing the predicted values 
of XspTn^™'"' with the experimental values for S = 2 and 
S = 5/2. Then the spin value S was determined more pre- 
cisely by fitting the respective CMC simulation to the exper- 
imental value of X^Tn^™'"'. Finally Ji was found by fitting 
the experimental value of Xsptn to the theoretical value for the 
respective simulation. 



Quantity 


Fit 1 


Fit 2 


5 


2.06 


2.64 


Ji 


207 K = 17.8 meV 


210 K = 18.1 meV 


J2 


= 41.4 K = 3.6 meV 


= 85 K = 7.3 meV 


J2/J1 


= 0.2 


= 0.4 


Jc 


= 21 K = 1.8 meV 


= 21 K =1.8 meV 


Jc/Jl 


= 0.1 


= 0.1 



experimental data. One cannot hope to obtain a good 
fit to an experimental x(T) data set by the theoretical 
predictions unless one can at least fit the experimen- 
tal ;^'naxj.max jatum. The quantities x"'''''Ji/Ng^^il, 
kBT"'^^/S{S+l)Ji, and ^'n^^T"''^^/3C are listed in Ta- 
ble VII for both the classical and quantum Monte Carlo 
simulations. Using the values of C obtained from Eq. (23) 
using g = 2, the predicted values of ^maxymax direct 
comparison to our experimental datum are listed in Ta- 
ble VII for the quantum value of S in the QMC simula- 
tions and for classical values S = 2 and S = 5/2 in the 
CMC simulations, respectively. One sees from the table 
that the value of -^^maxjimax very sensitive to the ratio 
J2/J1 but that it hardly changes for a given J2/J1 as the 
interlayer coupling ratio Jc/ Ji is changed over the range 
from to 0.1. 

Plots of •jj/maxjimax ygj-g^g J2/J1 are shown in 
Figs. 27(a) and 27(b) for = 2 and S = 5/2, respectively. 
The data were fitted by the second-order polynomials 
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(105) 

2 



for S — 2 and S = 5/2, respectively, as shown by the 
respective solid curves in Fig. 27, where the units of the 
fits are cm'^K/mol. The rms deviations of the fits from 
the data are 0.005 and 0.010 for 5 = 2 and ^ = 5/2, 
respectively. 

A comparison of the calculated values of T^maxjnmax 
in Fig. 27 and Eqs. (105) with the observed value 
^maxymax ^ Q gg cm^K/molMu in Eq. (20) indicates 
that the local moment model can reproduce the observed 
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FIG. 28: (Color online) Two fits (Fit 1 and Fit 2, dashed 
curves) of the high-temperature magnetic susceptibility x of 
BaMn2As2 from Fig. 7(a) by CMC simulations. The fits are 
only valid above Tn but are extrapolated to lower tempera- 
tures. The fit parameters are given in Table VIII. The tem- 
peratures of the breaks in slope of the fits are discernable and 
denote the predicted Neel temperatures in Eq. (96) for the re- 
spective parameters, which are somewhat above the observed 
value. 

-^maxjnmax yg^j^g -ftrith the following Combinations of pa- 
rameters 

• g = 2, 5* = 2 and J2/J1 w 0.17 

• .g = 2, 5* = 5/2 and J2/J1 ~ 0.36. 

In the J1-J2 model, the G-type AF magnetic structure 
that is observed in BaMn2As2 is stable against the stripe 
state as long as J2/J1 < 1/2 [Eq. (4)], which is satis- 
fied by both of these estimates. We do not have simula- 
tion data for precisely these two values of J2/J1. Also, 
the parameter set {S, Ji, J2, Jc} is underdetermined by 
the experimental susceptibility data, so we have to make 
choices about some of the parameters when we fit the 
experimental data by the available classical Monte Carlo 
data. We choose Jc/Ji = 0.1 because this value is in- 
dicated both from the neutron scattering fit in Table II 
above and from the theoretical results in Table IX below. 
For each of the two potential values S = 2 and 5/2, we 
use the respective CMC J2/J1 simulation in Table VII 
that shows the closest agreement with the experimental 
Xsptn^™'*'^ for that spin value, namely 

1. Fit 1: J2/J1 = 0.2, Jc/Ji = 0.1, 5 = 2 

2. Fit 2: J2/J1 = 0.4, JjJi = 0.1, S = 5/2. 

Next, we have a choice of how to obtain a precise fit to 
the magnitude of the experimental data by the simula- 
tion data for J2/J1 — 0.2 and 0.4, both with Jc/ Ji = 0.1. 
We could adjust the (/-factor, the orbital contribution to 
the susceptiblity, and/or the spin value. At this stage 



such changes are just fitting parameters, so we arbitrar- 
ily choose to adjust the spin value slightly to obtain a 
good numerical fit of the particular simulation to the ex- 
perimental value of X^in^™^'' given in Eq. (20). Then we 
fix the value of Ji by substituting the experimental value 
of x™t? = 0-80 X 10"^ cm^/mol Mn and 5 = 2 into the 
expression x^^^Ji/Np^g^^^ and equating that with the 
X'^'^^Ji/Ng'^fj.Q value Hsted in Table VII. The parame- 
ters obtained from the two fits are listed in Table VIII. 
Remarkably, the value of Ji is not sensitive to the val- 
ues of S, J2 or Jc, and a consistent value Ji « 210 K 
= 18 meV is obtained for both fits. The two fits are 
compared with the experimental data from Fig. 7(a) in 
Fig. 28. When plotting the fits, the calculated spin sus- 
ceptibility per mole of spins has to be multipled by two 
(two atoms of Mn per formula unit) and then added to 
the orbital contribution given in Eq. (18). These fits are 
only valid in the paramagnetic regime above Tn = 625 K, 
but they are extrapolated to lower temperatures. The 
quality of the fits to the experimental data is reasonable 
for both fits. Thus we cannot distinguish between the 
two possibilities S — 2 and S — 5/2 for the Mn spins on 
the basis of magnetic susceptibility measurements alone. 



XI. ^5 As NMR MEASUREMENTS AND 
ANALYSIS 

A. ^^As NMR Spectrum 

As shown in Fig. 1, each As atom is coupled to four 
Mn atoms. Thus through ^^As NMR one can probe the 
magnetism of the Mn sublattice in BaMn2As2. Figure 29 
shows typical ''^As NMR spectra in the magnetically or- 
dered state at different temperatures T < Tn for a poly- 
crystalline sample of BaMn2As2. At low temperatures, 
along with the most intense central line the spectrum 
contains extra shoulder-like features on either side. The 
broad linewidth is attributed to the random orientation 
of the internal field with respect to the external field in 
the powder sample. ^^As has an electric quadrupolar mo- 
ment that interacts with the local electric field gradient 
(EFG) in the crystal giving rise to the splitting of the 
NMR line. Thus in principle one should see in the '''^As 
spectra three allowed transitions: an = — i -s-> +i 
central transition, and two Iz — ±i f-s- ±| satellite 
transitions. Therefore in an attempt to fit the experi- 
mental spectra taking into account both the EFG and 
the isotropic spin shift effects, we find that the spec- 
trum at 4.2 K can be fitted reasonably well with iso- 
shift A'iso — 0.38%, quadrupolar frequency vq ~ 2.1 
MHz, width of central peak 0.43 kOe, width of satel- 
lite 1.13 kOe, and EFG asymmetry parameter 77 ~ 0.0. 
The fit is shown as the solid red curve through the data 
at 4.2 K in Fig. 29. The value of vq is comparable to 
that reported for BaFe2As2 in the ordered state. 

The linewidth and position were found to be almost 
temperature independent. As shown in Ref. 13 from 
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FIG. 29; (Color online) ''''As NMR spectra for polycrystalline 
BaMn2As2 at different temperatures. The solid red line is a 
fit to the spectrum at 4.2 K. 



magnetic neutron diffraction data, the sublattice mag- 
netization is nearly saturated at 300 K. Since the NMR 
linewidth in the ordered state is a measure of the sub- 
lattice magnetization, the independence of the linewidth 
over our temperature range is consistent with the neutron 
diffraction results. 

The internal field at the ^^As site can be analyzed by 
taking the crystal symmetry into consideration, which 
has been adopted in an analysis of the hyperfine field 
at the ^^As site in BaFe2As2 by Kitagawa et al.'^^ Ac- 
cording to their analysis, for a G-type antiferromagnetic 
spin structure the internal field at the ^^As site is zero 
due to a perfect cancellation of the off-diagonal hyperfine 
field produced by four in-plane NN Mn spins when the 
spin moments are parallel to the c-axis. Thus the spin 
components along this axis do not produce any magnetic 
broadening in the ^^As NMR spectra. Only the a6-plane 
components of the ordered Mn spin can produce an in- 
ternal field perpendicular to the c-axis at the ^'^As site. 
On the other hand, for a stripe-type AF spin structure, 
a c-axis component of the spin moments produces an in- 
ternal field Hint = z'BS along the a-axis, where z' is the 
number of nearest neighbor Mn spins of the ^^As site, B 
is the off-diagonal hyperfine coupling constant and S is 
the Mn spin. 

Assuming that the broadening of the NMR spectra 
originates from Hi^t at the ^^As site, Hmt is estimated 
to be ~ 215 Oe from the spectral width. Using the 
ordered moment ^ = 3.9 /xb/Mu and z' — 4, the off- 
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FIG. 30: (Color online) Nuclear spin-lattice relaxation rate 
(l/Ti) measured at the '^As site versus temperature T. The 
solid and dashed lines represent T'^ and behaviors, respec- 
tively. The slope of the former line fitted to the 50-300 K 
data is 2.51 x 10~''s-^K-^ Inset: (rir)-^ versus T. 



diagonal hyperfine coupling constant B is estimated to be 
^ 14 Oe/ for the case of stripe-type AF order. Such 
a small B is of the order of the nuclear-nuclear dipo- 
lar field and is not likely due to transferred hyperfine 
couplings. For the G-type AF structure, the a6-plane 
components can be produced by a canted component of 
the Mn spins when the magnetic field is applied perpen- 
dicular to the ordered moment axis, i.e., perpendicular 
to the c-axis. Using the perpendicular component of 
the spin susceptibility x± ~ 1-3 x 10^"^ cm'^/mol from 
Fig. 7(a) and H = 7.1 T that we used for measurements 
of the spectra, the a6-component of the ordered Mn mo- 
ment fiab is evaluated to be 0.0083 /ub/Mu in this field. 
Now using Hint = 215 Oe and fiab = 0.0083 /ib, the 
off-diagonal hyperfine coupling constant is calculated to 
he B ^ 6.5 kOe//iB. This value of B is comparable to 
B = 4.3 kOe/^B reported in BaFe2As2.^^ Thus our '''^As 
NMR spectra observed in the AF ordered state are con- 
sistent with the G-type AF structure reported from the 



neutron experiment 



13 



B. Nuclear Spin-Lattice Relaxation Rate 

The longitudinal nuclear magnetization recovery curve 
following saturation was fitted by the double exponential 
function^'' 



1 



M{t) 
M(oo) 



= 0.1e"*/^i +0.9e-^*/^\ 



as expected for the center line of the spectrum of the 
^^As nuclear spin / = |, where 1/Ti is the ^^As nuclear 
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spin-lattice relaxation rate and M{t) and M(oo) are the 
nuclear magnetization at time t after saturation and the 
equilibrium nuclear magnetization at time t = oo, respec- 
tively. The extracted l/Ti as a function of temperature 
is shown in Fig. 30, where 1/Ti is seen to increase rapidly 
with increasing temperature. In the AF state, this rapid 
increase in 1/Ti with T is a clear signature of relaxation 
due to scattering of magnons by the nuclear spins. Ac- 
cording to Beeman and Pincus,*'^ in the AF state for mag- 
netic insulators, 1/Ti is mainly driven by such magnon 
processes, leading to a power law T-dependence.^^"^^ For 
T ^ A/fce, where A is the anisotropy gap in the spin- 
wave spectrum, it either follows a behavior due to a 
two-magnon Raman process or a behaviour due to a 
three-magnon process, while for T <C A/fce, it follows a 
thermally activated behaviour 1/Ti oc r2g-A/fcBT^ 
seen from Fig. 30, our ''^As l/Ti data in the T-range 
50 < T < 300 K follow a behavior rather than a 
behavior indicating that the relaxation is mainly gov- 
erned by the two-magnon Raman process. A fit over 
this T range yields 



1 



(2.51 X lO"'' s-^K-^)T^ 



(106) 



as shown in Fig. 30. The lack of activated behavior down 
to 50 K indicates that A/fce is smaller than 50 K. 

For the two-magnon process, 1 /Ti is determined by the 
slopes of the spin wave dispersion relations at u; and 
thus by the spin wave velocities. The spin wave veloci- 
ties within the afe-plane and along the c-axis in terms of 
the exchange constants in the Ji-J2-Jc Heisenberg model 
arc given above in Eqs. (9). Since the spin wave veloc- 
ity depends on the direction of propagation, 1/Ti should 
also depend on the spin wave direction. Based on the 
1 /Ti expression for the two-magnon process reported by 
Beeman and Pincus,^^ we have calculated the 1/Ti for 
BaMn2As2 with the body-centered-tetragonal structure 
{lA/mmm) arising from the two spin wave velocities as 



(27r)3 



- 1 



dx, (107) 



where 2:' = 4 is the number of Mn nearest neighbors to 
a given '^^As site (see Fig. 1), 2; = 4 is the number of 
nearest-neighbor Mn spins interacting with a given Mn 
spin, i = ab with a = a or i = c with a = c, and 9 
denotes the angle between the local hyperfine field at the 
^^As site and the anisotropy axis (c-axis) which is the 
Mn ordered moment axis. We also have 

n 

where g is the electronic 5-factor, j4hf is the hyperfine 

coupling constant and 7n is the ^^As nuclear gyromag- 
netic ratio given by '^^7n/27r = 7.2919 MHz/T. Using 



g = 2 and A\^[ = 6.5 kOc//iB that was estimated from 
the spectrum analysis, one obtains 



- = 6.0 X 10^ s-\ 
a 



(108) 



For T ^ A/Zcb as in our temperature range 50-300 K 
where 1/Ti oc T^, the integral in Eq. (107) approaches 
its maximum value 7r^/6, so Eq. (107) reduces to 



ab 

-1 



A 

n 

= Cab 



z'zha^k^ sin^ 6 



rp6 



(109) 
(110) 



The ratio Rab/c of the relaxation rates for a6-plane and c- 
axis spin waves should be independent of T. If we assume 
that the hyperfine coupling A of the electronic spins to 
the nuclear spins is isotropic, then Rab/c obtained using 
Eqs. (9), (109) and (110) is given by 



Rab/c — 



{l/Ti)ab 
(l/ri)c 

6 



/cy / vc 

^aJ \Vab 



(111) 



4 \aJ 



{JJ..h)(l + Jj2.h} 
(1-2J2/Ji)(l + Jc/2Ji) 



Taking o = 4.15 and c = 13.41 A for the lattice parame- 
ters at 8 K," and J2/J1 = 0.29 and J^/Ji = 0.09 from 
the neutron scattering fit in Table II, Eqs. (Ill) yield 



R. 



ab/c 



13. 



(112) 



Thus the nuclear spin-lattice relaxation rate due to spin 
waves traveling in the a6-plane is much larger than that 
due to c-axis spin waves and we will therefore assume 
that Eq. (109) gives the observed 1/Ti to a good approx- 
imation. 

Since A^f was estimated in Eq. (108), one can obtain 

information on the exchange constants from the coeffi- 
cient of the T'^ fit in Eq. (106). Inserting hvab from 
Eqs. (9) into (109) gives 



AV = 1 (AV ^cY z'zHsm'e) 

kBj Cab\nJ \aJ 1927rA;B5'4 ^ ' 



2J2 
Ji 



A 

2Ji 



Our single fit parameter Cab in Eq. (106) can only be 
used to determine a single exchange constant or a sin- 
gle combination of them. We therefore estimate Ji using 

the above values J2/J1 = 0.29 and Jc/Ji = 0.09 de- 
rived from our inelastic neutron scattering experiments. 
In Eq. (113), we also use z' = 4, z = we take 5 to be 
the ordered spin (S) = 2 (from magnetic neutron diffrac- 
tion experiments),^"^ and (sin^ 0) = (1/2) JJ^ sin"^ 6* ^6* = 
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2/3 (i.e. considering an average over all angles). Using 
Eq. (113), we then obtain 

^ = 160K (S=2) (114) 
Kb 

Ji = 14 meV. 

This value is close to the value Ji = 16 nicV estimated 
in Table II for S = 2 from our neutron scattering data. 
If we take the spin to be S — 5/2, the value of Ji from 
Eq. (113) would be a factor of (5/4)"* = 2.4 times smaller. 

The overall temperature dependence of l/Ti in 
BaMn2As2 in Fig. 30 is similar to that reported for 
KMuFg.^* In KMnFa, a deviation from power law behav- 
ior was observed at low temperatures and l/Ti shows a 
broad maximum. This broad feature at low temperature 
was attributed to the effects of defects or extrinsic im- 
purities. Thus in BaMn2As2, the deviation of the data 
below 50 K from the higher-temperature fit in Fig. 30 
is likely due to relaxation associated with defects and/or 
extrinsic impurities. 

For a metallic system, one would expect a Korringa- 
like behaviour [{TiT)~^ = constant] as has been observed 
in (Ba,Ca)Fe2As2 (Refs. 59, 65, 66) and i?FeAsOi_:,F^ 
{R = La, Pr) (Refs. 67, 68) in the paramagnetic state. 
In these compounds, (TiT)~^ is also constant at low 
temperature below Tn due to their metallic character 
and increases sharply near Tn. In contrast, (TiT)~^ 
in BaMn2As2 (inset of Fig. 30) shows a gradual in- 
crease with increasing temperature signifying the insu- 
lating ground state of the compound. 

XII. BAND-THEORETICAL ESTIMATES OF 
THE EXCHANGE COUPLINGS 

The quantitative analysis of the magnetic interactions 
in real magnets is based mostly on density functional the- 
ory. To a large extent this theory is very similar to the 
Fermi-liquid theory of Landau, however, strictly speak- 
ing, it allows one to obtain only the total energy of the 
ground state, the distribution of charge and spin densi- 
ties, and other quantities that can be directly determined 
by these. Several notable exceptions (Mott insulators, 
rare earth systems, systems near quantum critical points) 
have been revealed but currently it is believed that for 
magnets of the Fe group the accuracy of the commonly 
used local density approximation (LDA) is acceptable for 
the description of the ground state properties including 
the equilibrium magnetic moments at T = K. 

While the numerical agreement between experimental 
and theoretical magnetic moments is often very good, 
there are certain cases when the local approximation 
numerically violates quantum mechanical laws. For in- 
stance, even in insulating systems or magnetic molecules 
where the total magnetic moment is close to an integer 
number of Bohr magnetons, the value obtained from den- 
sity functional theory is usually not an integer. This dis- 
crepancy is related to the fact that the wave function of 
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FIG. 31: (Color online) Total energy iJ of a system when the 
polar angle 9ij between two local magnetic moments i and j is 
varied in band theory. The energies of the ferromagnetic Efm 
and antiferromagnetic Eaf magnetic structures and their dif- 
ference are indicated. Here the antiferromagnetic state is the 
ground state. For a ferromagnetic ground state, the minimum 
in energy would be at 9ij = 0. 

the density functional method is often not an eigenfunc- 
tion of the square of total spin (even without relativistic 
effects). This effect of 'spin contamination' usually can- 
not be eliminated or easily resolved. While non-integer 
values of the moment in metallic systems are traditionally 
explained by itineracy of the system and partial occupa- 
tion numbers, the problematic issue of whether or not 
is an integral of the motion is usually ignored with the 
hope that such errors are small. The relationship between 
the single-particle spectrum obtained in the density func- 
tional theory and the physical properties of the magnetic 
excitations is not clearly defined. Nevertheless, the re- 
search of the last 20-25 years revealed that LDA often 
provides good agreement between theory and experiment 
for the magnetic excitation spectra if the ground state is 
properly described. 

The description of the intersite magnetic interactions 
represents a typical problem within the topic of magnetic 
excitations. By itself, the determination of a pairwise ex- 
change parameter Jij between atoms i and j in an arbi- 
trary magnetic material is not well-posed. For instance, 
in very itinerant systems the effective spin Hamiltonian 
can have very non-Heisenberg behavior. However, from 
the phenomenological theory of ferromagnetism^^'^" the 
energy of any weak and smooth variation of spin den- 
sity can be described by the effective classical Heisenberg 
Hamiltonian for equivalent classical spins Si and Sj with 
magnitudes 5*^ = Sj = S given by 

E = J2 Jv^^ ■ ^3 =^^J2 ^"-i ■ (11^) 

Traditionally, a set {Jij} of exchange coupling con- 
stants in the density functional theory can be calculated 
using two approaches. In the first approach, using the ef- 
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TABLE IX: Parameters of BaMn2As2 with NN (Ji), NNN (J2) and interlayer (Jc) exchange interactions obtained from density 
functional theory. Here S is the calculated spin, E is the total energy per Mn atom, FM means ferromagnetic structure and 
G-type AF structure is the Neel (checkerboard) antiferromagnetic structure with an in-plane structure shown in the top panel of 
Fig. 2. The two stripe structures have the in-plane antiferromangetic structure shown in the bottom panel of Fig. 2, where the 
Stripe- AF structure has AF stacking and the Stripe-FM structure has FM stacking along the c-axis. The estimated exchange 
constants J in rows 5 and 6 are calculated from Eqs. (119)-(121) using the total energy values in column 4. Our exchange 
constants using the LDA and GGA in the last two rows were calculated from the excitation energies from the magnetically 
ordered ground state. In the reference (Ref.) column, "PW" means "present work". 



Magnetic 


M 


S 


E/Mn 


2Ji + Jc 


Jc 


Ji 


J2 


Jc/Jl 


J2/J1 


Ref. 


Structure 


(Mb/Mu) 




(meV) 


(meV) 


(meV) 


(meV) 


(meV) 








FM 


2.74 


1.37 


-330 














12 


G-type AF 


3.20 


1.60 


-660 














12 


Stripe- AF 


= 3.20 


= 1.60 


-515 














12 


Stripe-FM 


= 3.20 


= 1.60 


-505 














12 










41.0" 


2.0" 


19.5" 


-5.4" 


0.10" 


-0.28" 


12, PW 










28.5'' 


1.2'' 


13.7'' 


-2.5'' 


0.09'' 


-O-IS*" 


12, PW 


FM (LDA) 


2.8 


1.4 






wO 


-9.1 


-2.2 






PW 


FM (GGA) 


3.0 


1.5 
















PW 


G-type AF (LDA) 


3.3 


1.65 




27.2 


1.03 


13.1 


2.8 


0.08 


0.21 


PW 


G-type AF (GGA) 


3.6 


1.8 




26.2 


1.0 


12.6 


2.7 


0.08 


0.21 


PW 



"Calculated using Eqs. (119)-(121) as written. 

''Calculated by replacing 5^ by 5(5 + 1) in Eqs. (119)-(121). 



fective Heisenberg model (115) one can solve for the 
from the set {Fa/3{Jij)} of equations for the differences 
of the energies between different magnetic structures a 

and (3 

^a0{Jij) ^^^2 V^^"^ 

obtained from band structure calculations as shown in 
Fig. 31. This is the usual way to obtain {Jij} for highly 
localized magnetic insulators and is usually the most suit- 
able method for the calculation of magnetic phase tran- 
sition temperatures. 

Another approach is based on the definition of Jij as 
the second derivative of the total energy in Eq. (115) 
with respect to rotation of moments from their magnetic 
alignment in a given magnetically ordered ground state 

1 d^E 

^ ^^^^^ 

which is proportional to the curvature of the total en- 
ergy E versus angle Oij near the minimum for an anti- 
ferromagnet at 9ij = 180° in Fig. 31. This definition of 
Jij corresponds to a linear response scheme and is usu- 
ally the most suitable technique for the analysis of the 
excitations above the ground state (spin waves) and is 
directly related to the dynamical magnetic susceptibil- 
ity measured in inelastic neutron scattering experiments. 
The procedure for evaluating Eq. (117) depends on the 
band structure methods and the specifics of the linear re- 
sponse method employed. This technique has been used 
for many magnetic materials in the past.^^ This approach 



can be understood as a static limit of the dynamic linear 
response technique which has been used for calculations 
of spin waves and Stoncr excitation spectra in magnets. 
One can analytically obtain an expression for the on- 
site stability parameter Jq which should be the same as 
^ Jij . A comparison of Jq and ^ Jij is a check on the 
consistency of the calculations and the reliability of the 
numerical scheme. 

We first consider the exchange constants obtained from 
total energy differences and then in the subsequent sec- 
tion from the energy of excitations from the magnetically 
ordered ground state. 



A. Exchange Interactions from Total Energy 
Calculations 

Using density functional theory in the LDA, An et al. 
correctly deduced from total energy calculations, prior 
to the availability of the experimental results, that the 
G-type antiferromagnetic structure of BaMn2As2 has a 
lower energy than cither the FM striicturc or of two types 
of stripe structure. Their predicted ordered moment for 
the G-type AF structure was = 3.20 /xe/Mn, somewhat 
smaller than the value of = 3.9(1) /ie/Mn observed 
later. Their LDA total energies and ordered moments 
for the FM and G-type AF structures are listed in Ta- 
ble IX, together with their total energies of two com- 
mensurate coUinear stripe states with the in-plane stripe 
structure shown in the bottom panel of Fig. 2.^^ The 
Stripe- AF structure has AF alignment of the ordered mo- 
ments along the c-axis, whereas the Stripe-FM structure 
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has FM alignment along the c-axis. As seen in Table IX, 
the ordered moment ji of the Mn in the FM structure is 
not the same as the value of /U in the G-type AF struc- 
ture. 

From the LDA total energies and ordered moments cal- 
culated by An et al.-^^ for the magnetic structures listed 
in Table IX, one can obtain estimates of the exchange 
couplings in BaMn2As2 using the value for the spin S 
obtained from the ordered moment fj, as 

S = — = — 

with g = 2. The classical energies per spin of the mag- 
netic structures in Table IX obtained using Eqs. (3) and 
Fig. 2 are 

Efm = S|m(2Ji + Jc + 2J2) (118) 

= 5|(-2Ji - Je + 2J2) 
-E'StripeAF = Sq{Jc — 2.J2) 
-E'StripeFM = Sq{ — Jc — 2J2). 

We have taken the Mn spin in the two stripe phases to 
be the same as in the G-type AF structure, since they 
were not given by An et al. Because the total energy 
contains a constant term proportional to the square of the 
magnetization, we solve for the exchange constants using 
only differences between these total energies according to 
Eq. (116). From the last two expressions in Eqs. (118) 
we obtain 

Jc = ^^(-BstripeFM — -BstripeAp)- (HQ) 

From the first two expressions in Eqs. (118) we obtain 
2J, + Je = ife-|^). (120) 

^ \ "^FM '-'G / 

Thus Eqs. (119) and (120) determine the two exchange 
constants Ji and Jc- Then from the second and fourth 
of Eqs. (118) we solve for J2 according to 

J Jl -E^StripeFM — -Eg /■, oi ^ 

•^^ = T 45| • 

It is not clear whether to retain S"^ or to insert the 
quantum mechanical expectation value S{S + 1) of (5*^) 
in place of 5^ in Eqs. (119) (121), so we calculate two sets 
of exchange parameters based on these two assumptions, 
which are given in Table IX. Using the second assump- 
tion, the values of Jc and Ji are respectively about the 
same as the values in Table II deduced from our inelastic 
neutron scattering experiments, but J2 has the opposite 
sign in the theory and experiment. 

We studied the properties of BaMn2As2 using den- 
sity functional calculations of the electronic structure and 
magnetic interactions in the FM and G-type AF struc- 
tures. For consistency, we used the experimental values 



of the lattice parameters a = 4.15 A and c = 13.47 A 
and the theoretically optimized value of the internal As 
parameter zas = 0.3524 utilized by An et al.^^ Our elec- 
tronic structure calculations were performed using the 
recently developed full-potential linear muffin tin orbital 
program. The accuracy of the exchange couplings ob- 
tained is about 2-3%. The studies of the exchange cou- 
plings were done using the static linear response tech- 
nique described in Refs. 73 and 74. 

Our results using LDA and the generalized gradient 
approximation (GGA) are very similar to those reported 
by An et al.'^^ We find that that BaMn2As2 has a rela- 
tively large ordered moment in the G-type AF structure 
with fj, = 3.3/XB/Mn in LDA and 3.6/XB/Mn in GGA, as 
listed in Table IX, with a small band gap (0.15 eV) in the 
electronic spectrum as observed, and with the G-type AF 
ordering having the lowest energy among all considered 
magnetic structures. The ferromagnetic structure has no 
charge gap, i.e., the compound would be metallic. Our 
total energy E-pM — Eq differences were 350 meV/Mn 
(LDA) and 375 meV/Mn (GGA), which are similar to 
the values of 330 meV/Mn and 380 meV/Mn obtained in 
Rcf. 12, respectively. 

While the magnetic moments are relatively large, they 
show a significant dependence on the magnetic structure, 
in agreement with the results of Ref. 12. For instance, 
the ferromagnetically ordered BaMn2As2 has an ordered 
moment of 2.8 ^b/Mu in LDA and 3.0 ^b/Mu in GGA, 
which deviate from the corresponding values for G-type 
AF ordering (Table IX) by about 20%. Due to this rel- 
atively strong dependence of the ordered magnetic mo- 
ment on the magnetic structure, estimates of exchange 
couplings from total energy calculations should be used 
with caution. 



B. Exchange Interactions from Excitations from 
the Magnetically Ordered Ground State 

For comparison with the exchange constants dediiced 
from inelastic magnetic neutron scattering experiments, 
calculations using the linear response technique^^ are 
preferable to the total energy technique, as noted above. 
Our calculations of the parameters of the Heisenberg 
model are: Ji — 13.1 mcV, J2 = 2.8 meV, Jc = 1.03 mcV 
(LDA) and Ji = 12.6 meV, J2 = 2.7 meV, Jc = 1.0 meV 
(GGA), as summarized in Table IX. These values are 
quite comparable with the; values deduced from our in- 
elastic neutron scattering measurements in Table II, and 
roughly similar to those in Table IV obtained from molec- 
ular field analysis of our magnetic susceptibility data in 
Fig. 7. 

The G-type AF ordering temperatures obtained from 

our spin value and exchange parameters in Table IX using 
the molecular field expression (41) are 

Tn = 730 K (G type AF, LDA) (122) 
Tn = 810 K. (G type AF, GGA) 



38 



These mean-field Nccl temperatures are somewhat larger 
than the observed value of 625 K, as expected, and indeed 
are approaching the temperature of the maximum of the 
measured susceptibility in Fig. 7(a) which from Table VII 
is of order the mean-field transition temperature. 

The other longer-range pair exchange parameters ap- 
pear to be much smaller, suggesting very short-ranged 
exchange interactions in this material. In particular, the 
difference between the above parameter Jg and ^ 
over six NN and four NNN is only about 5%, suggesting 
very short-ranged exchange interactions in BaMn2As2. 
This is different from the corresponding results for many 
Fe pnictides, where the exchanges with further neighbors 
are not so small and definitely provide a finite contri- 
bution to the spin wave spectrum. We attribute this 
difference to the metallic character of the Fe pnictides 
and the semiconducting character of BaMn2As2. 

To check the dependence of {Jij} on the type of mag- 
netic order we also performed linear response calculations 
of the {Jij} for the FM phase. The stability parameter 
Jo for this phase appears to be negative confirming the 
instability of such order with respect to the deviation 
of a single spin from the ordered moment direction (i.e., 
from 9ij — 0). This directly supports the qualitative 
behavior of the total energy versus 6ij in Fig. 31. The 
pair exchanges in this phase are Ji = —9.1 meV and 
J2 — —2.2 meV with a very weak coupling along the z- 
direction, which are compared in Table IX with the other 
exchange constant values discussed above. These results 
indicate that not only the ordered moments are different 
in the different magnetic phases, but the exchange cou- 
pling parameters depend on the type of magnetic order 
even in materials with a relatively large (3-4 /xb) mag- 
netic moment. 

Overall the localized Heisenberg model with four NN 
interactions Ji and four NNN interactions J2 in the ab- 
planc and two NN interactions Jc along the c-axis is suf- 
ficient to theoretically describe the magnetic properties 
of BaMn2As2 quite well. 



XIII. ORDERED MOMENT IN THE Ji-J-z-Jc 
HEISENBERG MODEL FROM SPIN WAVE 
THEORY 

As previously noted, an ionic picture suggests that the 
spin of the Mn^+ ion in BaMn2As2 is S = 5/2, yielding 
for g = 2 an ordered moment of 

(m) = gSuB = 5 MB . (123) 
On the other hand, the observed ordered moment 
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is only 3.9(1) /xb/Mu (Ref. 13) implying a substantial 
spin reduction (5). In view of the sizable frustrating 

AF ncxt-nearest-neighbor exchange J2 discussed above, 
it is natural to ascribe the moment reduction to enhanced 



quantum fluctuations. In the following, we shall use 
the conventional spin- wave theory to examine the quan- 
tum spin reduction for the layered Ji- J2- Jc square-lattice 
Heisenberg antiferromagnet. 

The spin- wave theory provides an expansion of the sub- 
lattice magnetization in powers of 1/S: 



S - (S) = no 



2S (25)2 + 



(125) 



where the leading correction no is determined by non- 
interacting spin waves, while higher order corrections 
nk>i come from magnon interactions. For the nearest- 
neighbor square-lattice Heisenberg antiferromagnet the 
two versions of the spin-wave expansion based either on 
the Dyson-Maleyev''® or the Holstein-Primakoff^^ repre- 
sentation of spin operators yield identical results: no = 
0.19660, ni = 0, and 712 = —0.0035, such that the se- 
ries (125) rapidly converges and compares extremely well 
with existing numerical results. 

Chandra and Doucot''* used the harmonic spin-wave 
theory to investigate the quantum renormalization of 
ordered moments for the next-nearest-neighbor J1-J2 
square lattice Heisenberg antiferromagnet. They found 
that the leading order correction no diverges as J2 — > 
Ji/2 due to a softening of the excitation spectrum seen 
in the first of Eqs. (13) above. This fact is considered 
as an indication of a quantum spin-liquid state around 
the strongly frustrated point J2 = Ji/2. Subsequently, 
Chakravarty, Halperin, and Nelson^*^ calciilated the next- 
order correction m, which becomes finite for J2 ^ 0, has 
an opposite sign compared to no, and also diverges at 
J2 = Ji/2. Below, we extend the results of Chakravarty 
et al. to a finite coupling between frustrated antiferro- 
magnetic layers. Comparison of two consecutive terms 
in the series (125) is necessary to judge the accuracy of 
the spin- wave expansion for large J2 . 

In the spin- wave calculations for the Ji- J2- Jc model (2) 
with H = we use a single-rotating-sublattice basis*° 
for the Neel structure with Q — {it, tt, tt) and apply the 
Holstein-PrimakofF transformation for spin operators ex- 
panded to first-order in 1/S. The quantum reduction of 
ordered moments in the harmonic approximation is given 
by an integral over the paramagnetic Brillouin zone 



no 



- — V 



2wk ' 



(126) 



where A'' is the number of spins. 



1 



Ak = l-j2(l-72k) + |jc , (127) 
-Bk = 7ik + ^ jc COS kz , Wk = Bl 

with j2 = J2/J1, jc = Jc/Ji, where Jioj^ is the magnon 
energy in units of 4 Ji5, and 



39 



Here the positive x and t/-dircctions arc defined to be _ J_ -Bk cos fcz (129) 

in the directions of the a and b primitive square-lattice N ^—^ 2wk 



translation vectors, respectively, the positive z-direction 
is perpendicular to the layers in the direction of the c 
lattice translation vector, and we have set a = 6 = c = 1. 

To treat the effect of magnon interaction one needs to 
introduce various Hartree-Fock averages of the bosonic 
operators compatible with the harmonic spectrum. For 
the present model this procedure yields in addition to no 
in Eq. (126) three other integrals: 



1 ^k72k 1 -Bk7ik Then, the leading nonlinear correction to the sublattice 

— ■ ^ — , Ai = — ^ — , magnetization in Eq. (125) is expressed as 



k 



1 



j2(m - Ai)7ik(l - 72k) + ^hicim - A2)(l - 72k) cosfc^ + ijc(Ai - A2)(7ik - cos A:^) 



(130) 



r 



TABLE X: The first-order and the second-order corrections, 
Eq. (125), to the ordered moment in the J\-Ji-Jc stacked 
square-lattice Heisenberg antiferromagnet. 





Jo 


= 


Jc = 


0.05 


Jc 


= 0.1 


J2/J1 


no 


ni 


no 


ni 


no 


ni 


0.00 


0.1966 


0.0000 


0.1427 


0.0087 


0.1260 


0.0088 


0.05 


0.2124 


-0.0047 


0.1529 


0.0071 


0.1347 


0.0075 


0.10 


0.2312 


-0.0117 


0.1648 


0.0045 


0.1447 


0.0055 


0.15 


0.2542 


-0.0227 


0.1790 


0.0003 


0.1566 


0.0024 


0.20 


0.2828 


-0.0405 


0.1963 


-0.0065 


0.1709 


-0.0028 


0.25 


0.3198 


-0.0710 


0.2177 


-0.0182 


0.1885 


-0.0115 


0.30 


0.3698 


-0.1278 


0.2455 


-0.0392 


0.2111 


-0.0270 


0.35 


0.4423 


-0.2485 


0.2832 


-0.0809 


0.2414 


-0.0573 


0.40 


0.5605 


-0.5691 


0.3389 


-0.1784 


0.2853 


-0.1264 


0.45 


0.8074 


-1.9842 


0.4354 


-0.5045 


0.3594 


-0.3483 


0.49 


1.6005 


-23.899 


0.6257 


-2.6777 


0.5008 


-1.7482 



For vanishing interlayer coupling Jc 
comes 



0, Eq. (130) be- 
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7?k(l-72k) 



Ak72k -7fk 



2wk 



7?k 
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as in Eq. (A3) of Chakravarty et al7^ To obtain Eq. (131) 
we replaced the factor (m — Ai) in Eq. (130) by its inte- 
gral representation from Eqs. (129) and i?k by 71k from 
Eqs. (127), and set jc = in ^j^. 

Numerical values of uq and ni for a range of ratios 
J2/J1 = to 0.49 and for Jc/Ji = 0, 0.05 and 0.1 
are listed in Table X. These allow one to compute 
the spin reduction S — {S) for any value of S. The 
spin reduction for S — 5/2 is plotted in Fig. 32 for 



Jc/J2 = and 0.1. Three-dimensional effects generally 
suppress quantum fluctuations, as seen in a comparison 
of Figs. 32(a) and 32(b), and extend the validity of the 
spin-wave expansion (125) to somewhat larger values of 
J2/J1 ^ 0.44, although the series remain divergent at 
J2/J1 - 0.5. 

The above fits to our experimental neutron scatter- 
ing and magnetic susceptibility results for BaMn2As2 
and the band theoretical estimates of the exchange pa- 
rameters in this compound yielded Jc/Ji ~ 0.1 and 
J2/J1 = 0.2-0.4. According to the second-order cal- 
culations in Fig. 32(b), together with Eq. (124) with 
(7 = 2, this parameter regime predicts an ordered mo- 
ment reduction of « 0.34-0.52 ^b/Mu due to quan- 
tum fluctuations.*^ This result appears to rule out the 
possibility that the spin of the Mn is S = 2 because 
the ordered moment would then be a maximum of ~ 
3.66 /xb/Mu for 5 = 2, which is significantly smaller 
than the observed value^^ of 3.9(1) /itB/Mn. On the other 
hand, if 5* = 5/2, then the corresponding predicted or- 
dered moment is ^ 4.66 /xb/Mu, which is too large com- 
pared to the observed value. 

It seems likely that charge and/or magnetic moment 
amplitude fluctuations which arise from both on-site and 
intersite interactions, can account for the additional re- 
duction needed to reach agreement with the observed or- 
dered moment for a Mn spin S = 5/2. For example, the 
ordered moments of the Mn atoms are 3.50(4) /xs/Mn in 
SrsMngAsaOa,"'^^ 4.15(3) mb/Mu in La2Mn2Se203,*2 
and 4.04(8) hb/Mu in Ba2MnMo06,*^ all containing 
Mn+^ ions with nominal spin S = 5/2. Such reductions 
are also often attributed to covalency. On the other hand, 
for the more ionic compound MnF2 , the ordered moment 
of 4.82 /iB/Mn+^ is much closer to the value of 5 /xe/Mn 
expected for S = 5/2 with g = 2,** consistent with ex- 
pectation. 
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FIG. 32; (Color online) Reduction S — {S} in the ordered spin 
(5') from its value S in the absence of quantum fluctuations 
versus J2 / Ji according to linear spin wave theory to first order 
(solid black curves) and second order (filled circles) in 1/S for 
(a) Jc = and (b) Jc/Ji = 0.1. 



XIV. SUMMARY AND CONCLUSIONS 

Our anisotropic magnetic susceptibility x versus tem- 
perature T measurements from 300 to 1000 K of sin- 
gle crystals of BaMn2As2 yielded the Neel temperature 
Tn = 618(3) K, close to the value of 625(1) K previously 
determined from neutron diffraction measurements on a 
polycrystalline sample. ^'^ The xC^) above Tn is nearly 
isotropic, indicating that single-ion anisotropy effects are 
small and that a Heisenberg model for the spin interac- 
tions is appropriate. Below Tn, the x becomes strongly 
anisotropic, with x± nearly independent of T and x\\ 
dropping nearly to zero for T — 0, which corresponds 
qualitatively to the textbook behavior for collinear an- 
tiferromagnets in molecular field theory (MFT). How- 
ever, the temperature dependence of x above Tn con- 
tinues to increase, rather than decrease as expected 
from MFT, indicating the presence of strong short-range 
AF order above Tn. Such short-range AF order above 
Tn is expected for a quasi-two-dimensional spin lattice 
as in BaMn2As2. Magnetic inelastic neutron scatter- 
ing measurements were carried out on a polycrystalline 



BaMn2As2 sample at 8 K with momentum transfers up 
to 6 and energy transfers up to 140 meV. These 
data allow estimates of the magnetic exchange interac- 
tions in this compound to be made using appropriate 
models. We also report ^^As NMR measurements in the 
antiferromagnetically ordered state of a polycrystalline 
BaMn2As2 sample from 4 to 300 K. The nuclear spin- 
lattice relaxation rate is found to obey the power law de- 
pendence 1/Ti oc T"^ from 50 to 300 K which we interpret 
in terms of the exchange interactions in this compound. 

We developed various theories for the J1-J2-JC Heisen- 
berg model in order to model our experimental data 
and extract values of the exchange constants between 
Mn spins and the value of the spin. Our inelastic neu- 
tron scattering measurements indicate that this is the 
minimal model needed to understand these data. For 
G-typc antifcrromagnetic ordering shown for the square 
spin lattice in the top panel of Fig. 2, and where the c- 
axis alignment is also antifcrromagnetic, linear spin wave 
theory at T ^ Tn was used to calculate the magnon 
dispersion relations in Sec. IV A 1. The in-plane spin 
waves soften as J2 increases, and become unstable for 
J2 > Ji/2, signaling a phase transition to the in-plane 
stripe state shown in the bottom panel of Fig. 2. Thus 
the G-type AF ordered state requires J2 < Ji/2. This 
theory is used in Sec. IV B to fit our inelastic magnetic 
neutron scattering data at 8 K, a temperature far below 
the Neel temperature of ~ 625 K, and obtain estimates 
of SJi, SJ2 and SJc for BaMn2As2, where S is the spin 
on the Mn atoms that is not determined separately in 
the spin wave fit to the data. From the ordered mo- 
ment fi = 3.9(1) i-iB = gSfiB, one would estimate = 2 
for g = 2. On the other hand for the ion Mn+^ one 
would estimate a high-spin 5* = 5/2. In Sec. VIII Bl 
we also calculated the spin wave contribution to the low- 
temperature heat capacity for comparison in Sec. VIII B 2 
with our previously published^ ^ experimental heat capac- 
ity data for a single crystal of BaMn2As2. We also used 
spin wave theory to extend the nuclear spin-lattice re- 
laxation rate 1/Ti calculations of Beeman and Pincus 
for the isotropic cubic Heisenberg spin lattice^^ to the 
J1-J2-JC model. We obtained Eqs. (109) and (110) that 
were used to analyze the fit to our ^^As 1/Ti NMR data 
with 1/Ti cx T3 from 50 to 300 K for BaMn2As2. 

A molecular field theory (MFT) treatment of the Ji- 
J2- Jc Heisenberg model was described in Sec. VII. In the 
paramagnetic state the system follows the Curie- Weiss 
law C/{T + 6) for T > 6, which has the same form as 
described in many textbooks for the Ji-only model. The 
ratio / of the Weiss temperature 6* to Tn, f — O/Tj^i, is 
found to be / = 1 for J2 = 0, as expected for a bipartite 
spin lattice, but is / > 1 for frustrating AF J2 > and 
is / < 1 for nonfrustrating reinforcing FM J2 < 0, which 
are intrasublattice interactions so the spin lattice is no 
longer bipartite. Thus for J2 > 0, the Curie- Weiss law 
continues to be followed below T = 6 down to T = Tn, a 
characteristic already noted by Ramirez for geometrically 
frustrated antiferromagnets.'^^ As shown in Fig. 8, we find 
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that X\\{T) for T < Tn strongly depends on J2, whereas 
X_L is independent of T and J2 at T < Tn, apart from 
the implicit influence of J2 on Tn- We further find that 
the staggered moment and the magnetic heat capacity 
versus T/Tn at T < Tn are also independent of J2, again 
apart from the implicit infiuence of J2 on Tn- 

We carried out quantum (QMC) and classical Monte 
Carlo (CMC) simulations of both x{T) and the magnetic 
heat capacity Cniag(2^) in H = versus Jc/Ji and J2I Ji- 
Most of the QMC simulations were for J2 = due to 
severe negative sign problems when J2 was taken to be 
positive, which is antiferromagnetic and frustrating for 
G-type AF order. When we replaced the square of the 
spin, 5^, in the CMC simulations by the quantum me- 
chanical expectation value S{S + 1), the QMC simula- 
tions for J2 = for increasing S merged smoothly with 
the CMC simulation (which corresponds to S ^ 00) as 
shown in Fig. 23, so we used the CMC simulations to 
fit the experimental x{T) data for T > T^. The CMC 
simulations of Cniag(T') as a function of ,/c showed AF 
phase transitions at temperatures Tn that increased with 
Jo > but decreased with J2 > 0, as shown in Fig. 22. 
The Tn( Jc, J2) data are well-fitted by Eq. (96). 

We also carried out band theoretical estimates of the 
exchange couplings in BaMn2As2. There are two generic 
ways to do this. The first is to take the differences be- 
tween the total energies for different spin configurations 
such as in Eqs. (119-121), where the lowest energy spin 
configuration is the ground state. This method is of- 
ten used to determine exchange constants to be used in 
the calculation of magnetic transition temperatures and 
yields the exchange constants in rows 5 and 6 of Ta- 
ble IX which will not be further discussed. The second 
is to measure the change in the total energy due to small 
deviations of the spin directions from the magnetically 
ordered ground state via Eq. (117), which gives the ex- 
change values in rows 9 and 10 of Table IX. These values 
are considered to be more reliable for comparison with 
values extracted from inelastic neutron scattering exper- 
iments. 

Our exchange constants from the MFT fit to our 
anisotropic x(^) data for BaMn2As2 below Tn in Fig. 28 
and Table VIII are probably not reliable because that fit 
assumes that x(^) follows the Curie- Weiss law above Tn 
(i.e., that there are no two-spin AF correlations stronger 
than 1/T above Tn), which is strongly violated by the 
data in Fig. 7(a). Similarly, although our fits by MFT 
to the ordered moment Jlz{T) from neutron diffraction 
measurements^^ on BaMn2As2 in Fig. 14 and to the heat 
capacity Cp(T) (Ref. 11) in Fig. 15 are reasonable, they 
are not sufficient to distinguish between the possible spins 
S = 2 and S = 5/2 discussed above. Furthermore, the 
fit to the Cp (T) data near room temperature by the sum 
of the Debye lattice heat capacity and the MFT pre- 
diction of the magnetic heat capacity indicated that the 
measured magnetic heat capacity is too low. This dis- 
crepancy suggests the presence of strong AF fluctuations 
above Tn that reduce the magnetic entropy and magnetic 



TABLE XI: Summary of our most reliable exchange constants 
in BaMn2As2 obtained from band theory and from fitting our 
experimental data by our predictions of the Ji-J^-Jc Heisen- 
berg model. The notation "=" means that the value that 
follows it was assumed, not fitted. Two spin values for the 
neutron fit are listed because the neutron fit gives the product 
of the spin S and the respective exchange constant, not the 
two separately. For the x fit, a range of spin values can fit 
the data. We chose two spin values close to 2 and 5/2, cor- 
responding to J2/J1 = 0.2 and 0.4, respectively. Two band 
theory estimates were obtained using the LDA and the GGA, 
as noted. 



Data 


S 


Ji 
(meV) 


J2/J1 


Jc/Jl 


neutrons 


2 


16.5 ± 1.5 


0.29 ±0.05 


0.09 ±0.02 




5/2 


13.2 ± 1.2 


0.29 ± 0.05 


0.09 ± 0.02 


X(T > Tn) 


2.06 


17.8 


= 0.2 


= 0.1 




2.64 


18.1 


= 0A 


= 0.1 


■^'^As NMR 


= 2 


14 


= 0.29 


= 0.09 


(T < Tn) 










band theory 










LDA 


1.65 


13.1 


0.21 


0.08 


GGA 


1.8 


12.6 


0.21 


0.08 



heat capacity below Tn, consistent with the behavior of 
the magnetic susceptibility above Tn in Fig. 28. The cal- 
culated T^ contribution to the heat capacity at low tem- 
peratures from spin waves, without an anisotropy gap in 
the spin wave spectrum, is about 40% of the measured 
value. However, an anisotropy gap would suppress the 
spin wave contribution exponentially to zero at low tem- 
peratures. 

We have gathered together in Table XI our most re- 
liable exchange constants in BaMn2As2 from our band 
theory calculations and from the theoretical fits to our ex- 
perimental data by the Ji- J2- Jc Heisenberg model. Sev- 
eral features arc noteworthy. First, all three exchange 
constants are consistently positive (antiferromagnetic). 
Second, the estimates give similar values of Ji « IS- 
IS mcV for S values in the range from 2 to 5/2. For the 
susceptibility fits, nearly the same Ji was obtained from 
the two fits despite the significant differences between 
the respective S and J2/J1 values. Third, the estimates 
of J2/ Ji from band theory and from analysis of the neu- 
tron and magnetic susceptibility measurements are in the 
range 0.2-0.4, which are below the value of 0.5 at which 
the in-plane G-type AF order would classically become 
unstable with respect to the stripe AF order [sec Fig. 2 
and Eq. (4)], and are therefore consistent with the ob- 
served G-type AF order. From our classical Monte Carlo 
simulations of the heat capacity of stacked square lattice 
layers, the exchange parameters from the neutron scat- 
tering fit predict Tn « 640 K if the Mn spin is 5* = 5/2, in 
close agreement with the experimental value of « 625 K. 

Finally, with the above range of exchange parame- 
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tcrs, our second-order spin wave calculations in Sec. XIII 
show that the ordered moment reduction chie to quan- 
tum fluctuations alone is at least ~ 0.4 //b/Mh. Because 
the measured ordered moment is (/u) = 3.9(1) /xb/Mh, 
this argues against assigning a spin S = 2 to the Mn+^ 
ions which gives (n) = gSfiB = 4 /xb/Mu and favors 
S = 5/2 for which one would obtain {fx) = 5 /xb/Mh 
for g( = 2 in the absence of quantum fluctuations. The 
additional reduction needed to reach the experimental 
value is likely due to charge and/or magnetic moment 
amplitude fluctuations which arise from both on-site and 
intcrsite interactions, and/or from hybridization eff'ects, 
consistent with the reduced ordered moment measured 
for other Mn"*"^ compounds. ^^'^^^^^'^"^ This effect cannot 
be described in the Heisenberg model formalism used in 
this paper. For instance, in Sec. XII we discussed that 
despite the large moment of the Mn, the magnitude of 
the magnetic moment can vary by ^ 20% depending on 
the specific magnetic configuration, suggesting that am- 
plitude fluctuations of the magnetic moment may indeed 
be relevant. 

As noted in the introduction, the view that BaFe2As2 
is an itinerant antiferromagnet is not universally held. 
Furthermore, the results of many magnetic inelastic neu- 
tron scattering measurements on BaFe2As2 have been an- 
alyzed in terms of local moment Heisenberg models, even 
when the authors of this modeling believe that the itin- 
erant model is valid. The reason for this latter analysis, 
as has been stressed in the literature, is that the mag- 
netism of itinerant models can often be parametrized by 
local moment Heisenberg models. Further review and 
discussion of this issue is given in Ref. 4. 

We therefore now compare the exchange constants in 
BaMn2As2 with those in the isostructural (at room tem- 
perature) high-Tc AFe2As2 parent compounds {A = Ca, 
Sr, Ba) within the context of the Ji-J2-Jc local mo- 
ment Heisenberg model. The j4Fe2As2 compounds or- 
der into an in-plane stripe-type antiferromagnetic struc- 
ture below ~ 200 K (lower panel of Fig. 2) and the 
lattice distorts to orthorhombic symmetry at or above 
Tn-^ Within the orthorhombic structure and assuming 
S = 1/2 and gS^B w 1 /xb, one defines the average 
(Ji) = (Jia + Jib)/2, yielding (Ji)/J2 = 0.7 1.4 for a 
variety of AFe2As2 compounds,'' which is in the regime 
(Ji)/J2 < 2 expected for in-plane stripe- type order- 
ing [see Eq. (4)]. These values can be compared with 
those for BaMn2As2 in Table XI where J1/J2 ~ 3. In 
BaMn2As2 and also in the AFe2As2 compounds, the in- 
terlayer coupling is weak compared to the in-plane 
couplings and the systems should be considered to have 
strongly spatially anisotropic exchange, but not nec- 
essarily two-dimensional. Thus, the ^T2As2 systems, 
where T is a 3rf-transition metal element, appear to 
be ideal systems to study the physics of the Ji-J2-Jc 
Heisenberg model, including the possibility of tuning 
the system to and through the quantum critical point 
J1/J2 = 2 by doping. Indeed, doping-dependent stud- 
ies of Ba(Fei_a;Cra:)2As2 have recently revealed a tran- 



sition from stripe-type AF order to G-type AF order at 
X « 0.3,*^ and studies of Ba(Fei_2;Mnj;)2As2 have re- 
vealed a transition to a new state, possibly due to com- 
petition between G-type and stripe-type AF ordering, at 



Acknowledgments 

We are grateful to Andreas Kreyssig for insights about 
spin wave theory and to Ferenc Niedermayer for confirm- 
ing the calculation of xiT) (Ref. 52) in Eq. (97). Work 
at the Ames Laboratory was supported by the Depart- 
ment of Energy-Basic Energy Sciences under Contract 
No. DE-AC02-07CH11358. R.J.M. would like to thank 
F. Trouw for assistance with Pharos. The work has ben- 
efitted from the use of the Los Alamos Neutron Science 
Genter (LANSE) at Los Alamos National Laboratory. 
LANSGE is funded by the U.S. Department of Energy 
under Contract No. W-7405-ENG-36. The work of A.H. 
was supported by the Deutsche Forschungsgemeinschaft 
through a Heisenberg Fellowship (Grant No. HO 2325/4- 
2). M.E.Z. acknowledges Grant No. ANR-09-Blanc- 
0211 SupraTetrafer from the French National Research 
Agency. 

Appendix A: High- Temperature Series Expansions 

From quite general considerations, one can show that 
the results of MET at high temperatures for the mag- 
netic susceptibility (the Curie- Weiss Law) is an exact 
result arising from a quantum mechanical treatment of 
local moment magnetism using a high-temperature se- 
ries expansion HTSE. We also discuss a complemen- 
tary "Curie- Weiss Law" for the magnetic heat capacity 
at high temperatures, which is useful when discussing our 
Monte Carlo simulations of the magnetic heat capacity 
in Sec. IX. 



1. Magnetic Susceptibility and the Curie- Weiss 
Law 

Using the fluctuation-dissipation theorem, one can ex- 
press the diagonal aa components Xa (o^ = y, z) of the 
magnetic susceptibility tensor in terms of the two-spin 
correlation functions 

^ {s^s:) (Ai) 

where r is the distance measured in the number n of 
bonds, including zero, of spin Sr from a typical cen- 
tral spin 5*0.^'^^ In the isotropic Heisenberg model, one 

obtains^^ 

r 
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If one only considers the single-spin autocorrelation func- 
tion (r = 0), then one has = {S^) = {S'^)/3 = 
S{S + l)/3 which gives the Curie law 



iVgV|5(5+l) 



(A3) 



which in turn is the Curie- Weiss law (22) with ^ = and 
Curie constant (23). 

Writing r = n = 0, 1, 2, in terms of the dis- 
tance of a spin in number of bonds from the central spin 
at position (i.e., n means the n"^-nearest-neighbor of 
the spin at the origin in terms of the minimum number 
of bonds between them), one can express the temper- 
ature dependences of the two-spin correlation functions 
as high-temperature series expansions in 1/T with the 
general form 



-I- 



n,n-\-l 



(fceT/J)" (/cbT/J)" 



-1-1 



-I- 



(A4) 



where = S{S + l)/3 as noted above and the lowest- 
order term for a given Ff^ is Substituting the first 
three terms of Eq. (A4) into Eq. (A2) gives 



X = 



S{S+1) , Ff^i ^ Ff 



+ 



+ 



C 
T 

+ 



1 + 



knT/J {kBT/jy 



ri,; 



S{S + 1) V^bT/J {kBT/jy 



(A5) 



where the Curie constant C is the same as in Eq. (23). If 

one keeps only the first two terms in the square brackets 
and uses the Taylor series expansion l-t-j;~ 1/(1 — a;) 
for small x to put the quantity in square brackets into 
the denominator, one gets the Curie- Weiss law (22) with 
Weiss temperature 



S{S+l)kB' 

Comparing Eqs. (25) and (A6) gives the coefficient 
_ z[S{S + l)]^ 

1,1 ~ n • 



(A6) 



(A7) 



2. The HTSE for the Magnetic Heat Capacity 



capacity Cmag- We show this to be the case [Eq. (A9)] 
and apply the result in Sec. IX. 

From Hamiltonian (21), the thermal-average magnetic 
configuration energy in zero field only depends on the 
nearest-neighbor two-spin correlation function as^^ 

i^mag(T) = ^NzJ{So ■ Si)T(r), (A8) 

where the factor of 1/2 is introduced to avoid double 
counting the distinct AF NN bonds and (• • •)t denotes a 
thermal average of the quantum mechanical expectation 
value. The magnetic specific heat Cmag(T') is obtained 
by diflferentiating Eq. (A8) with respect to T. Thus the 
first 1/T HTSE term of E^^giT) gives the first term in 
the HTSE of Cmag(r) as a l/T^ term. 

Rushbrooke and Wood showed that the first six terms 
of the HTSE of x{T) and Cn,^g{T) of a Heisenberg spin 
lattice containing equivalent spins S can be expressed 
in terms of the identical NN exchange couplings J and 
the bond connectivity ( "lattice parameters" ) of the spe- 
cific spin lattice. With respect to the present discussion, 
they found that the first (l/T^) term of the HTSE for 
Cniag(T') is independent of the type of spin lattice and of 
the spin lattice dimensionality and only depends on z, S, 
and J according to^^ 



mag 



R 



JS{S + 1) 



(A9) 



This term is the same for FM and AF interactions be- 
cause the exchange constant is squared, which gives a 
positive-definite result for Cniag(T). Thus when compar- 
ing calculated Cmag(^) data for lattices with the same 
coordination number z but different exchange constants 
and/or spins, a universal high-temperature behavior is 
obtained if the data are scaled according to 



a 



mag 



k^T 



R 



versus 



JS{S+1) 



(AlO) 



The Curie law for the magnetic susceptibility arises be- 
cause there is a nonzero susceptibility for isolated spins. 
This is modified at high temperatures by the addition of 
a 1/r^ term in Eq. (A5) arising from spin interactions, 
yielding a Curie- Weiss law with nonzero Weiss tempera- 
ture 0. On the other hand, the magnetic heat capacity 
of isolated spins is zero, and hence there is no equivalent 
Curie law for Ej^agiT) or Cmag(T): the values are just 
zero. Equation (A9) can therefore be considered to be a 
"Curie- Weiss law" for Cmag(T)- 



We discussed above that the Curie- Weiss law for the 
magnetic susceptibility of equivalent spins is rigorously 
derived from the first {1/T) term in the HTSE of the 
nearest-neighbor two-spin correlation function and hence 
does not depend on the particular crystal structure or 
spin lattice dimensionality. This suggests that there is 
an analogous term in the HTSE of the magnetic heat 



Appendix B: Anisotropic Susceptibility below Tn for 
the Ji-Ji-Jc Model in Molecular Field Theory 

1. Perpendicular Susceptibility xx below Tn 

Setting the external field H to zero and using the ex- 
change fields in Eq. (29), the average exchange energy of 
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FIG. 33: (Color online) Influence of a perpendicular mag- 
netic field Hx on the sublattice magnetizations of an ordered 
antiferromagnet. The Hx tilts the ordered sublattice mag- 
netizations Ml and M2 that are initially pointed along the 
2-axis by an angle 7 towards the applied field Hx along the 
X-axis. The angle 7 is greatly exaggerated for clarity. 



the spin system is 
EcKch — —Ml • Bi 



-(A,M2 + AdMi-M2). (Bl) 



We do not add a second term — M2 • B2 to this, because 
that would double-count the exchange interactions be- 
tween the spins which occur pairwise. 

We apply an external magnetic field Hj_ that is perpen- 
dicular to the ordered moment direction in the antiferro- 
magnetically ordered state that induces a perpendicular 
magnetization Mj, in the sample, where M_l = Mi +M2 
is the vector sum of the sublattice magnetizations. From 
Fig. 33, the angle between Mi and M2 is 

cos(f>^ cos(180° - 27) = - cos(27) « - (l - 27^) , 

where we used cos{A — B) = cos A cos -B -I- sin A sin i? and 
on the right-hand side we used the small angle approxi- 
mation cos a; « 1 — 2:^/2. We assume that 7 is very small 
because x± = ^^^n^iH^^o M±/ H± by definition. Thus the 
exchange energy in Eq. (Bl) becomes 



£^cxch = M'[-As-fAd(l-272)], 



(B2) 



where Mi = M2 = Mi. 

A perpendicular magnetic field = H±i causes 

the ordered AF spin sublattices to tilt towards the ap- 
phed field direction, away from the ordered moment z- 
direction, as shown in Fig. 33. The magnetic energy due 
to the perpendicular field is 



-2M„; ■ H 



-2M,H±_ sin 7 ss -2M,H^j, (B3) 



where we have used the small-angle approximation 
sinx w X. Thus the total magnetic energy is 

E = i;cxch + E^= Mf[-K + Ad(l - 27^)] - 2M,Hj_-t. 

(B4) 



The stable configuration minimizes the energy. Taking 
the derivative of E with respect to 7 and setting it to 
zero and setting Ad = — |Ad| because Ad is negative gives 



7 = 



2|AdlM, 



Thus the interactions within the same sublattice (As, i.e., 
J2) have no influence on this equilibrium condition. The 
equilibrium value of the component Mj_ of the total mag- 
netization in the direction of H^, using the small-angle 
approximation smO 9, is 

M_L = 2M,sin7 « 2Ma = 

|Ad| 

which gives the perpendicular susceptibility as 

Mi_ 1 



XI- = 



Hi_ \Xd\' 



Note that Mi and hence also its temperature dependence 
have dropped out, so that x± in this treatment is inde- 
pendent of T below Tn . 
From Eqs. (54) we have 



I Ad 



W + f) 
c 



Thus we obtain 

1 _ C 

M ~ Tj,{l + f) ~ Tn + i 



c 



(B5) 



x(7^n), (T<Tn) 



which is Eq. (55) in the text. Although J2 is not present 
explicitly, its influence is expressed through the implicit 
dependence of Tn on J2. 



2. Parallel Susceptibility X\\ below Tn 

Here again we assume that the susceptibility in the ab- 
sence of explicit exchange couplings is xo = /T, which 
is isotropic above Tn- We apply a small magnetic field H . 
Below Tn a large exchange field develops as seen by each 
sublattice because of the ordered moments. Therefore we 
must use a Brillouin function to describe the magnetiza- 
tion of each sublattice. 

The saturation magnetic moment of a spin S is 

For N spins, the saturation magnetization is therefore 
Msat = NgStJL^. (B6) 

The magnetization of the spins is written 

M, = Ms^tBsiy) (B7) 

where Bs{y) is the Brillouin function given by 



Bs{y) 



1 

2S 



{(25 



Dcoth 



(25+1 



y 



coth 



■)}• 

(B8) 
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where < Bs(y) < 1 for y > 0, 



y 



(B9) 



and the jf-factor is usually set to the value g = 2. 

In MFT, we replace H in the Brillouin function by 
the magnetic inductions Bi in Eqs. (28) and (29), which 
include the exchange fields. Thus we have 



Ml 



Mo = 



(BIO) 



^M,^tBs 



giiB{H + XdMi+X,M2) 



Substituting for Ag and Ad from Eqs. (54) gives 



Ml 11 ff^B 

= —Bq 

M,at 2 



H 



TnU+I) 

c 



Mo 



knT 



(Bll) 



Ma _ 1 f^/^B 

M,at " 2 ^ 



H- 



C 



Ml 



Tn(/-1) 



Mo 



To simplify the notation and the solution to Eqs. (Bll) 
we define reduced magnetic field, temperature and or- 
dered moment variables, respectively, by 



h 
t 



gfJ-BH 



Mi 



kBTn{H = 0) 
T 

Tn{H - 0) 

Msat 



(B12) 



gSjiB 



Then using M, = Nii,j2 and Afsai = N^^at = NgSfiB 
and the expression for the Curie constant C in Eq. (23), 
Eqs. (Bll) become 



Ml2 



Bs 



h 3(/-1)mi. 3(/ + l)/i2. 



t 2(5' +1) t 2(5+1) t 



(B13) 



M2z = Bs 



h 3(/ + 1) a^u 
i 2(5+1) i 



3(/ - 1) 
2(5+1) i 



This solution is valid in both the paramagnetic and anti- 
ferromagnetic states. The reduced parallel susceptibility 
per spin is 



X(t) = lim^. 
h^O h 



(B15) 



In the AF phase, the order parameter is the staggered 
ordered moment 



mI = — Y- 



(B16) 



which is one-half the difference between the z- 
components of the ordered moments of the two sublat- 
ticcs. The term "ordered moment", when used in the 
context of a coUinear antiferromagnet, is the staggered 
moment. 

The susceptibility is isotropic at Tn (at reduced tem- 
perature t = 1). Therefore setting t = 1 and jliz = Ji^z 
in Eqs. (B13) gives 



'^-5TT^^' 



where i = 1,2. Using the expansion Bs{y) = (5 + l)y/3 
for y <C! 1 and solving for /ij^ gives 



5'+l - 

~^-^' = ^^ = W^) 

Then Eq. (B15) gives 

XII (r) XII (t) 3(/ + l) 



lim 



Hz{t,h) 



X||(Tn) X||(i = l) 5+1 h^o h 



(B17) 



(B18) 



Appendix C: Ordered Moment versus Temperature 
below Tn 



In the AF state, setting the applied magnetic field h = 
and the ordered moment fl2z{t) = —flizit) in the first 
of Eqs. (B13) gives the simple result 



4(i) = Bs 



5 + 1 



t 



(Cl) 



Thus the reduced exchange field as in Eqs. (28) and (29) 
is 



^exch(^) — 



5 + 1 



(C2) 



For specified 5, /, h and t, one can solve these two simul- 
taneous equations numerically for p,iz{t, h) and /22z(i, h). 
The average reduced magnetization per spin is 

Mz(i, h) = \[i^iz{t, h) + /X2.(i, h)]. (B14) 



Appendix D: Zero-Field Magnetic Heat Capacity 



below Tn 



In the presence of the staggered exchange field with z- 
components ffiexch = — -f^2exch and nonzero /ziz = —ijl2zi 



46 



the energy of a coUinear G-type AF system in zero ap- 
plied field is 

N 1 

Eave = —-^IJ'lzHlexch = — -^MizHi exch, (Dl) 

where the factor of 1/2 is included to avoid counting ciach 
magnetic moment twice (once in /n^ or Mi^ and again 
in -ffiexch)- The exchange field seen by sublattice 1 can 
be written using Eqs. (29) and (42) as 



Hi 



exch 



(A, - Ad)Mi 



.g2^|5(5+l 



(D2) 



Then using the expression for the saturation moment 
A*sat = gSi^B, one can rewrite this as 



^ 1 exch — I r» , -1 I 9 • 
\ S+1 J Msat 



Inserting Eq. (D3) into (Dl) gives 



' 2(5+1) 



Miz(T) 

Msat 



(D3) 



(D4) 



Using the dimensionless reduced variables introduced in 
Eqs. (B12) and the definition (B16) of the staggered mo- 
ment, we obtain 



Eave{t) 



SNkuT^S 

' 2(5 + 1) 



(D5) 



By setting equal to Avogadro's number A'a and us- 
ing the definition of the molar gas constant R = A^aAjb 
one obtains from Eq. (D5) the molar magnetic energy 



-E'ave (t) 

The molar magnetic heat capacity is then 



Cmag(0 _ 1 dEavejt) 
R R Tkr dt 



35 .,..d-^l{t) 
5 + 1^^^*^ dt ■ 



(D6) 



(56) 



Appendix E: Low- Temperature Heat Capacity of 
Spin Waves 

At low temperatures, only the lowest energy spin waves 
contribute to the heat capacity, so we can use generic 
Eq. (7) for the dispersion relation. To evaluate the in- 
tegrals in Eq. (69), we change variables in the integrals 
from wave vector q to the vector e with dimensions of 
energy and with components 



(El) 



Now the dispersion relation (7) is written symmetrically 
as 



ex 




Co 






= hVyQy, 




= hvyir/b 








— hv^ir/c. 



(E2) 



and Eq. (69) becomes 



NVs, 



(E3) 



de. 



de. 



Ep 



in which the anisotropy in the dispersion relation (7) has 
been moved to anisotropies in the limits of integration 
and in the pref actor. 

At low temperatures, only the lowest energy spin wave 
states are populated, so we can take the limits of each 
integral to be — oo to oo, which also eliminates the 
anisotropy between the limits of integration of the three 
integrals. We can then convert the integrals over the 
three Cartesian coordinates to an integral over radius in 
spherical coordinates according to e and 



/oo /"OO /"OO /■ 

d^x / dey / d^z ^ 47r / de e^. 
-oo J —oo J — CO Jo 



Now we will integrate only about the F point, so we must 
multiply by two to take into account the low-energy spin 
wave branches at the corners of the Brillouin zone of the 
primitive tetragonal direct lattice as discussed in the text. 
We then obtain 



E.: 



■ de. (low T) (E4) 



TT^K'vxVyV, Jo e^/'^B^ _ 1 
Changing variables in the integral to a; = e/k^T gives 

TT^TfvxVyVz Jo - 1 
The integral is 7r^/15, yielding 

Ibn VxVyVz 

Then setting TV = Na (Avogadro's number), the mag- 
netic heat capacity due to the spin waves per mole of 
spins is 

^mag 1 dEs,yc ( 47r^fc|V"spin ^ „3 (rp _^ rp \ 

"^ = :R^ = li5nwJ ' ^"^^"^^^^"'^ 

where R =- Npji-Q is the molar gas constant. 

We now calculate the two-dimensional spin wave 
theory prediction of Cmag to check consistency with 
Eq. (101) that was derived using chiral perturbation the- 
ory. In two dimensions (i.e., Jc = 0), the area of the 
sample is ^ = Na^ , where N is the number of spins and 
ah = is the area of the square unit cell which contains 
one spin in its basis. Then Eq. (69) is replaced by 



E — 



1 r^'' huj 

ffa'i J -Hi a J-ir/b c 
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Following the same steps as for the three-dimensional 
case above, converting the two-dimensional integrals to 
polar coordinates according to 



J — oo J — c 



f-QO 

2-K / de e, 
Jo 



and multiplying by two to take into account the spin wave 
excitations at the corners of the Brillouin zone, gives 



mag 

R 



6C(3) 
^{hv/a)'- 



(E8) 



where v is the isotropic spin wave velocity in the a6-plane 
and we have used -^§^dx = 2^(3). Equation (E8) is 
identical to Eq. (101) obtained for the isotropic Heisen- 
berg square lattice from chiral perturbation theory.^^ 
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